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ABSTRACT 

We present SpitzerlRAC (~2 deg^) and MIPS (~8 deg^) observations of the Cepheus 
Flare which is associated with the Gould Belt, at an approximate distance of ~300 pc. 
Around 6500 sources are detected in all four IRAC bands, of which ~ 900 have MIPS 
24 lira detections. We identify 133 YSO candidates using color-magnitude diagram 
techniques, a large number of the YSO candidates are associated with the NGC 7023 
reflection nebula. Cross identifications were made with the Guide Star Catalog II and 
the IRAS Faint Source Catalog, and spectral energy distributions (SED) were con- 
structed. SED modeling was conducted to estimate the degree of infrared excess. It 
was found that a large majority of disks were optically thick accreting disks, suggesting 
that there has been little disk evolution in these sources. Nearest-neighbor clustering 
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analysis identified four small protostellar groups (L1228, L1228N, L1251A, and L1251B) 
with 5-15 members each and the larger NGC 7023 association with 32 YSO members. 
The star formation efficiency for cores with clusters of protostars and for those without 
clusters was found to be ~ 8% and ~ 1% respectively. The cores L1155, L1241, and 
L1247 are confirmed to be starless down to our luminosity limit of Lboi = 0.06 Lq. 

Subject headings: infrared: general — ISM: individual (Cepheus Flare) — stars: for- 
mation 



Introduction 



The Gould Belt is a band of stars and mole c ular clouds t hat encircle the sky at an inclination 
■^20° to the Galactic Plane (jHerschell 



of ~20° to the Galactic Plane (jHerschell 118471 : iGouldl Il879l ). It is the locus of star formation 
within 140-500 pc of the Sun and includes many, well-know n star formation regions (Serpens, 
Ophiuchus, Orion, etc). The Spitzer Gould Belt Survey (SGBS: lAllen et al.ll2009l ) is a Spitzer GQ- 



4 leg a cy pr oject designed to extend the earlier successful Spitzer Cores to Disks (c2d; lEvans et al 



2003 



20081 ) legacy project and complete a census of star formatio n within 500 pc. In addition 



SGBS and c2d data complement data from the upcoming JCMT (jWard-Thompson et al.l 120071 ) 
and Herschel (jAndre &: Saracend l2005l ) Gould Belt surveys. The Spitzer Space Telescope is a 
85-cm- diameter cryogenicall y cooled satellite telescope designed to operate throughout the infrared 



regime ( Werner et 



(IRAC: lFazio et al. 



al.ll2004bl ). Its two instruments used in this study are the Infrared Array Camera 
20041') which can observe at 3.6-8.0 /im, and the Multiband Imaging Photometer 



for Spitzer (MIPS: lRieke et al.ll2004l ). which can observe at 24-160 /im. In this paper we focus on 
Spitzer data from one part of the Gould Belt known as the Cepheus Flare region. 

The Cepheus Flare is a complex of nebulae th at extends 10-20° out of the plane of the Galactic 
disk at a Galactic longitude of ~ 110° (|Hubblelll934l l. Star formation towards the Flare can be broken 
down into that associated with the wall of the Local Bubble at ~160 pc, that associated with the 
sweep of the Gould Belt at ^^300 pc, and that associated with the Perseus arm of the galaxy at 



800 pc (lYonekura et al.lll997l : iKiss et al.ll2006l ) . Figure [J shows a visual extinction ( Dobashi et al 



20051 ). CO (IDame et al.ll200ll ). and IRAS 100 //m (|Miville-Deschenes fc Lagachell20n,5l ) finding chart 
towards the Cepheus Flare. It reveals five associations of dark clouds, each located at the middle 
distance, including L1148+L1152+L1155, L1172+L1174, L1228, L1241 and L1247 + L1251. In 
this paper, we present IRAC and MIPS data obtained by the SGBS program towards these dark 
cloud associations and the related dark cloud L1221 which is outside the region shown in Figured) 

The whi te contours in Figu re [D show the integrated ^'^CO emission surrounding the dark cloud 
associations (jPame et al.l 120011). Table [D li s ts the cloud positi ons along with their inferred masses 
and ^'^CO line widths from lYonekura et al.l (1993). iKunI ( 19981 ) matched sources in the IRAS point 
source catalogs with Ha data to identify a catalog of pre-main sequence stars and dense cores in 
Cepheus. Their distribution of sources showed that star formation was occurring along the cloud 
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edges. Similarly, Figure [T] shows that regions of high extinction are also distributed towards the 
edges of the clouds (the lowest CO contours). 

The Cepheus Flare is an expansive and sparsely filled region. Studies have attempted to 
define structure within the Flare by grouping the scattered objects into associations and groups. 
Yonekura et al.l (119971 ) mapped Cepheus in CO and identified several large as sociations of dark 
clouds. The clouds presented in this paper correspond to Yonekura Group A. iKiss et al.l (j2006l ) 
studied the cloud morphology in Cepheus using their own 256 square degree extinction map. They 
identified eight cloud complexes across the wider area. 

The Cepheus Flare is bounded by a series of shells and loops. Three of these are shown in 
Figure [U as dashed lines. The Cepheus Flare Shell is an expanding supernova bubble at a distance 
comparable to the clouds in this paper. Its center is located at / ^ 120°, b ^ 17° (approximately 
7° east of L1241) and it has a radius of ~ 9.5° (50 pc at 300 pc) (jCrenier et al.lll989l : loiano et al 



20061 ). It is possible that the star formatio n within L1251 has been triggered by the passage of this 



shell (jOlano et al.l 120061: iKun et al.l 120081 ). The dark cloud L1228 is coincident with the current 
radius of the shell (IKun et al.l l2008l ) . Other star formation regions associ ated with the Cepheu s 
Flare Shell include L1333 {I = 129°, b = +13°) at a distance of 180 pc (jobavashi et al.l Il998l '). 
L 1333 is in the neighb oring constellation of Cassiopeia, on the opposite side of the Cepheus Flare 
Shell (jKun et al.ll2008l ). The Ce pheus Flare She ll and the older, larger Loop III supernova shell 
appear to be concentric spheres (jKun et al.l 120081 ). 



The Cepheus Shell should not be confused w ith the Cepheus Bubble which is an e xpanding 
dust ring that surrounds the Cep 0B2 association ( Patel et al. 1998 : Abraham et al. 2000). Spitzer 
observations of three clusters associated with Cep 0B2 found that 10% of detected disks were 
"transition" objects (their spectral energy distribut ions, SEDs, are essentially photospheric except 
for an infrared excess at the longest wavelengths; ISicilia-Aguilar et al.l 120061 ). Clouds coincident 
with this Bubble are excluded from ou r analysis as they are believed to be associated with more 
distant material ( Abraham et al. 200C). Other star formation groups in the direction of Cepheus 
include S140, Cep 0B3 and Cep 0B4, These are part of the Yonekura Distant Group, which is 
at a distance of 600-800 pc (jYonekura et al.l 119971 ). The infrared loop G109+11 was identified by 



Kiss et al.l (j2004l ) who later a ssociated it with "Void #2" in their extinction survey and a bright 



rim of excess 12 fiui emission (jKiss et al.ll2006l ). 



In Section, [2] we discuss distance estimates to the Cepheus dark cloud associations. In Section, 
[3] we outline the observation strategy, the data reduction procedure, and present the assembled 
catalog and false-color images of the region. In Section HI we use several techniques to identify YSO 
candidates based on their infrared colors. We classify the candidates based on their spectral indices 
and present color-color diagrams of the resulting list. In Section [5l we analyze the photometry of the 
YSO candidates, adding additional data where available, to produce spectral energy distributions. 
From these we calculate the Cepheus Flare luminosity function and conduct basic SED modeling. 
In Section [6l we look at the relation of the YSO candidates to their surroundings and present 
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extinction and MIPS 160 fj,m maps as well as the results of clustering analysis and comparisons 
with existing surveys. In Section 7, we discuss a star formation scenario for the Cepheus Flare 
region. In Appendix [Aj we discuss each of the individual dark cloud associations and compare 
results with archival 850 //m submillimeter continuum data. 



Distance 



The Cepheus Flare contains a number of different compo nents whose distr ibutions on the sky 
overlap, but whose velocities and distances are different (see Kun et al. 20081 fo r a re view). The 
dark cloud associations shown in Figure [1] were selected from the iDobashi et al.l (j2005l ) extinction 
map because they each had peak Ay > 3 and are within a distance of -^500 pc. The selected 
clouds divide into a Galactic West Group comprising the L1172+L1174 and L1148+L1152+L1155 
associations at a Galactic longitude of 101-105° and a Galactic East Group comprising the L1228, 
L1241, and L1247+L1251 associations at a Galactic longitude of 111-115°. There is a noticeable 
lack of regions with Ay > 3 in t he span 106-110° although it does no t mean that this span is 
completely devoid of young stars (iTachihara et al.ll2005l : iKun et al.ll2008l ) or of molecular material 
as evidenced by the CO and extinction maps. However, the lack of Ay > 3 extinction means that 
this region falls outside of our selection rules for Spitzer imaging. The arc of the Cepheus Flare 
Shell would appear to divide the active Galactic East star formation group and the barren central 
molecular mass of the Cepheus Flare. Thus the Cepheus Flare Shell appears to be triggering star 
formation as its sweeps across the Flare. 

The Galactic West Group of dark cloud associations is surrounded by common low-level inte- 
grated molecular emission suggesting that they are at a similar distance (e.g., see Figured]). Many 
of the distance estimates to these clouds ultimately rely on distance estimates to the Herbig AeBe 
star HD 200775 embedded within L1174. This is the driving source of the NGC 7023 reflection 
nebula, but estiniates of its distance have been hampered by uncertainties in its spectral type. 
Whitcomb et al.l (|l98ll ) surveyed the state of the art at the time of their publication and found 
a range of distances from 350 to 600 pc bef ore adopting: a distance of 440 pc. This distari ce was 
used by many studies of molecular cores (e.g. iMvers et al.lll987l : IWard-Thompson et al.lll999l'). This 
value is in agreement with the Hipparcos distance of 430lgQ'' pc to HD 200775 ( van den Ancker et al 



1997 



Bertout et al.lll999l ). 



Estimates of the di stance to L1174 based on diagrams of c olor excess-distance modulus give 
estimates of 300 ±20 pc (|Shevchenko et al.lll989l ') and 288 ±25 pc (jstraizvs et al.lll992l ). The second 
estimate is based on a group of reddened stars i n the proximity of the nebula, but does not include 
the actual distance to HD 200775. For that star. IStraizvs et al.l (|l992l ) derived alternative distances 
of 275 pc and 417 p c and adopted the former based on the correlation with their first estimate. 
Straizvs et al.l (jl992l ) estimated a distance to the LI 148/55 molecular ri ng of 325±13 p c. Given 
the uncertainty of the HD 200775 distance ( Kun et al. 2008) we adopt the Straizvs et al.l distances 
to L1174 and L1148±L1152±L1155. 
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The Galactic West Group of dark cloud associations are separated from the G alactic Eas t 
Group of associations by ~ 10°. At the adopted distance of these clouds this is ~ 50 pc. iKunI (jl998l ) 
used Wolf Diagrams to estimate distances of 200l2o° for L1228 and 3001^° P C for L 1241 and 
L1247/L1251. The distance to L1251 matches the distance adopted bv lLee et al.l (|2006l ). Despite 
their proximity on the sky, L1228 and L1251 are actu ally on opposite surfaces of the projected 
Cepheus Flare Shell and thus have diffe rent distances (IKun et alJ 120081 ) . For L1221, we use the 
distance adopted by lYoung et alJ (120091 ) . Table [T] lists the adopted dis tances of the dark cloud 
associations. The mass of each cloud, Mdoud) from 
adopted distances and is listed in column 7. 



Yonekura et al.l (119971 ) has been adjusted to our 



3. Observations and Data Processing 

The areas with Ay > 3 identified from the Dobashi maps were cross referenced with the 
Spitzer data archive and the "isolated cores" lists from the c2d program to create two lists of 
targets. The first list comprised those regions for which public data existed that was compatible 
with the processing requirements of the c2d/SGBS analysis pipeline and the second list comprised 
those regions that had not yet been observed with Spitzer. It was this "second list" that became 
the basis for Spitzer observations undertaken specifically for the SGBS. Table [2] lists the complete 
set of Spitzer observations (Astronomical Observation Requests - or AORs) included in this paper. 
The first column gives the name of the central cloud and the second column lists the program 
identifier (PID). The majority of the archival observations come from the c2d (PID #139) isolated 
cores program with the exception of L1228 which was taken as part of the Galactic First Look 
Survey (PID #104). The new regions are listed under the SGBS PID #30574. The extents of the 
regions listed in Table [2] are shown in Figure [1] by the rectangular footprints (white for IRAC, black 
for MIPS). 



3.1. Spitzer Gould Belt data 

Two epochs separated by 5-6 hours were used to take two complete sets of observations of the 
new SGBS clouds. The redundancy in coverage allowed for the rejection of transient phenomena 
(including foreground asteroids), the recovery of coverage lost by the blanking of bad pixels, and 
full mapping at MIPS 70 /xm where half the array was not working (each epoch was offset by half 
the array width). The unique AOR numbers for each epoch's observations are listed in columns 3 
and 4 for IRAC and columns 7 and 8 for MIPS. For the archival observations, only one epoch was 
obtained. Columns 5 and 9 list the respective dates the IRAC and MIPS observations were taken. 
Likewise, columns 6 and 10 list the version of the Spitzer Science Center (SSC) data pipeline that 
had been used to process the data prior to them being downloaded from the archive. 

The total area mapped was 2.11 square degrees with IRAC and 8.33 square degrees with MIPS. 
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MIPS mapped a larger area than IRAC because its minimum scan-lengths were longer than the 
diameter of the compact regions of high extinction. The regions for which we have four-band IRAC 
fluxes are slightly smaller than the areas listed in Table [TJ This is because the IRAC detectors are 
paired into two different pointing offsets resulting in slightly different coverage patterns between 
3.6/5.8 /im and 4.5/8.0 //m. The IRAC observations have a total integration time of 48 seconds 
per point split, between the two epochs. Each epoch consisted of two dithers. The SGBS MIPS 
observations were taken with the fast-scanning mode using a 240 arcsec step size and a cross scan 
overlap. This gave a total per point integration time of 31.4 seconds at 24 and 70 fim and 6.1 
seconds at 160 fiia. 

The basic calibrated data (BCD) were downloaded from the Spitzer archive and then ingested 



into the c2d/SGB S analysis frai nework (for a full description see lEvans et al.l 120071 or separately 
Harvey et al.ll2006 for IRAC and lRebull et al.ll2007l for MIPS). In brief, the data were inspected and 
custom data masks created to identify bad pixels. The data were then corrected for instrumental 
effects including bleed over, saturation, and banding effects for IRAC and iailbar and stim flash 



latents for MIPS. The improved data were mosaiced with the MOPEX package (jMakovoz et al 



20061 ) and s ource extraction wi th the c2dphot tool was performed independently at each wavelength 
and epoch (jEvans et al.l 120071 ) . T he separate epoch/wa velength source lists were band- merged 
together with the 2MASS catalog (jSkrutskie et al.l 120061 ) and cross identifications were made at 
better than 2 arcsec accuracy. Each source was then characterized spectrally and for the quality 
of detection. The catalog was "band- filled" to produce upper limit flux estimates for sources that 
were not detected at all wavelengths (jEvans et al.l 120071 ). 



Table [3] lists the original survey/instrument, filter, nominal wavelength, flux zero-point Sopt, 
and limiting magnitude of the ten bands from the band-merged catalog (2MASS J-K^; IRAC 
3.6-8.0 //m and MIPS 24-160 /im) plus additional bands from the second HST Guide Star Catalog 
(GSC-II, B-I), the IRAS point source catalogs (12-100 ;^m), and SCUBA data archive (450- 
850 //m). Data for these additional bands were only added for YSO candidates, as described in 
Section 15. 1[ The photometric system is based on t he Vega magriitude system using the fiux zero 
points taken from the c2 d delivery docume ntation (jEvans et al.l 120071 ). These are within 1-2 per 
cent of the SSC's IRAC0 dPeach et ahlbood ) and Mipjf zero points. We use th e c2d zero points to 
maintain compatibility with the c2d delivery documentation (jEvans et al.ll2007l ). The uncertainties 
for individual fiux measurements are > 4% and > 8% for IRAC and MIPS (24 /xm and 70 //m) 
respectively. These uncertainties do not include the zero point / absolute calibratio n uncertainities 
of 1.5% for IRAC, 4% for MIPS 24 /im, 20% for MIPS 70 //m dEvans et al.ll2007l l. The absolute 
Spitzer calibration uncertainties are lower than the individual measurement uncertainties by a 
factor of 2 or more, except for the MIPS 70 /im band which has a significantly higher calibration 
uncertainty. The limiting magnitudes are identified in Section 13.41 as the turnover in the source 



1http://ssc 


spitzer . caltech. edu/irac/calib/| 




^http: //ssc 


spitzer . caltech. edu/mips/calib/| 
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count distributions. 



3.2. Archival Spitzer data 

The c2d data used in this paper comes from their final data releas^ (DR4) and uses the 
same target selection and observation strategy as the SGBS program. A significant difference was 
that the c2d cores data were taken using the MIPS small map photometry mode rather than the 
fast-scanning techniq ue. As a result, the 70 fim coverage on the c2d cores can actually be smaller 
than the IRAC area (lEvans et al.l 120071 ). No 160 nm data were taken for the c2d isolated cores 
(requesting 160 /im data imposes a scheduling limitation as it can only be taken in particular "cold" 
campaigns). The L1148-|-L1152-|-L1155 region contains 3 c2d cores in close proximity so the entire 
area was covered by a single SGBS MIPS fast-scan to obtain complete 160 /im coverage. The only 
regions for which 160 //m data were not taken were L1221 and L1251. 

At the time of writing, the c2d studies on the individual first list cores have already been 
published or are in the process of being published. L1148 was studied by lKauffmann et al.l (j2005l ) 
who discovered a new very low lumin osity (VeLLO) Class YSO. L1155C, L1152, and L1228N 
were studied bv lChapman et al.l (|2009l ) to examine the extinction law from 3.6 to 24 /xm. Although 
L1228N was observed by c2d, we use the First Look Galactic data-set as it co vers both L1228N ari d 
L1228S. An analysis of the smah YSO group i n L122 8S was published bv IPadgett et al.l ^M)- 
L1221 was extensively studi ed bv lYouiig et al.l (j2009l ) including follow-up observations with the 
VLA. L1251B was studied bv lLee et al.l (j2006l ) who compared their YSO dete ctions with cont inuum 
submillimeter data. This was then followed-up by a molecular line study (iLee et al.l 120071 ). The 
dataset was built with the aim of maximizing coverage with a (r elatively) uniform sen sitivity. We 
chose not to add deeper data (e.g., the PID #3656 data used by IChapman et al.ll2009l ) as we wish 
to have all regions covered to approximately the same depth. 



The LI 152, L1155, a nd LI 221 studies used data from the DR4 c2d data release (jChapman et al 



200S 



Young et al.ll2009l ). the same one used by this paper. The L1251 study, however, used data 



from an earlier data release (DR3) that were processed through version Sll of the SSC pipeline 



(Lee et a 



20061'). Likewise, t he First Look study of L1228 used data from version S9.1 of the SSC 



pipeline (IPadgett et al.ll2004l ). The SSC pipeline is constantly being revised to improve calibration 
and reduce instrumental artifact^. Improvements in the calibration will change the photometry 
for regions where a newer pipeline has been used than in the original study. This is most significant 
for the L1228 region where it was found that fiuxes measured from S14/16 data had increased by 
an average of 16% against the S9 version of the data. The greatest increase was at 70 fim which 
increased by almost 40%. 



• ^http: //ssc . spitzer. caltech. edu/legacy/c2dhistory. html I 

^See |http : //ssc . spitzer . caltech . edu/ archanaly/plhistory/ for a revision history. 
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3.3. Images of extended emission 



Figure [2] shows RGB color composites of 4.5 fim (blue), 8.0 fim (green), and 24 fim (red) 
towards the regions observed for the SGBS. The green haze that covers a large area of the images 
comes from 8.0 fim background which is believed to arise predominantly from 7.7 /im PAH emission 
(jFlagey et al.ll2006l ). In L1241, the haze appears along the edge of the image and is anti-correlated 
with regions of higher 160 /im emission (see Section [6]). While the mosaic processing removed a 
number of instrumental effects, it was unable to correct for the bleeding at 8 ^m that appears 
in the color composite images as green smearing upwards from the bright center of NGC 7023. 
The red point source i n the L1148 map is the heavily reddened L1148-IRS source studied by 
Kauffmann et al.l ^20od ). 



The brightest extended structure at all wavelengths in this survey is the NGC 7023 reflection 
nebula. Figure [3] shows three-color composite images covering the visible (Digitized Sky Survey-II; 
left image), infrared (IRAC; middle image) and far-IR/submillimeter wavelengths (MIPS; right 
image). The extended cold dust emission, seen in silhouette in the first panel and emission in the 
last panel, forms a faint cross with HD 200775 at its intersection. The vertical axis of the cross is 
formed by the north-south filament that includes L1172 and L1174. The east-west axis is formed 
by the reflection nebula forming in this filament. The outflow from HD 200775 is currently inactive, 
but it has generate d an asymmetric east-west bi-conical cavity that is filled with hot atomic gas 
(|Fuente et al.lll998l ). The larger western lobe of this cavity appears as a hole in the left hand panel 
through which background stars are clearly visible. The ea stern lobe is trunca ted by dense material 
which is being heated and photo-dissociated by the YSO (|Fuente et al.lll998l ). 



The reflection nebula is confined to the north by a dense filament. This same structure is 
shown in strong 3.6 fim emission on the leading edge of colder dense material revealed by 850 /im 
dust emission (see Appendix [A|) . To the south of NGC 7023 is the L1172D dark nebula that appears 
in silhouette in the visible image, but reveals three embedded sources in the mid-infrared image. 
The dust itself emits at the longest wavelengths and shows the bright point source associated with 
the L1172 SMM 1 embedded protostar. The extended 8 fj,m emission (the red haze in the central 
image) is anti-correlated, however, with the high visual extinction in the first image and the dense 
gas in the third image. 



3.4. Source Catalog 

Tables H] and [5] list the number of objects detected in various combinations of the IRAC or 
MIPS bands within the SGBS Cepheus Delivery Catalog, which contains a subset of the data 
included in this paper. Specifically it includes the five regions L1148, L1172, L1228, L1241 and 
L1247+L1251. The other regions (L1155, L1152, and L1221) are available from the SSC as c2d 
data products. Source detection for SGBS data is defined as an object with a peak intensity >3 
X the local rms in its respective band. A 2MASS detection constitutes an object seen at >10 x 
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the local rms in both the H and Kg bands (a design constraint of the 2MASS catalog itself; see 
Skrutskie et"al]l2006l ^. 



Table U] shows that a total of 71076 sources were detected in at least one IRAC band and of 
those 6518 were detected in all four bands. Table [5] lists the number of source detections for each of 
the four IRAC bands and the first two MIPS bands for several different signal-to-noise thresholds. 
It shows that the majority of these sources were predominantly detected at shorter wavelengths. 
An important distinction between Tables H] and [5] is that the former only refers to the area common 
to all Spitzer wavelengths whereas the latter lists statistics for the entire region of the Delivery 
Catalog. 

Figure m shows source count histograms of the number density of sources per magnitude interval 
per square degree for six bands from the delivery catalog. The gray curve shows all sources. The 
limiting magnitudes listed in column 5 of Table [3] are taken as the magnitudes where the black curves 
in Figure m turn over. The total source count histogram at 3.6 //m is double-peaked. The 4.5 fiui 
source count histogram shows a second peak, but the effe ct is far less pronounc ed. Double-peaks 
we re also found i n the s ource count histo g rams for Perseus (|j0rgensen et al.ll2006l ) and Chamaeleon 
II (jPorras et al.l 120071 ). |j0rgensen et al.l (120061 1 showed that the double peak was not due to the 
presence of background galaxies by comparing counts of on-cloud and off-cloud sources and that it 
must derive from something, like extinction, that is common to the entire region. The solid black 
curves on Figure [4] plot only those sources that meet our detection criteria. The peaks of the black 
curves for 3.6 ^m (IRACl) and 4.5 ^m (IRAC2) are coincident with the first of the double-peaks. 

Over-plotted on Figure H] are the sourc e count histograrns of Galactic infrared sources to- 
wards Cepheus estimated from the model of IWainscoat et al.l (|l992l ) updated by Carpenter (pri- 
vate communication) to apply to the Spitzer bands. The source counts agree with the model at low 
magnitudes before diverging at high magnitudes. The divergence from the model increases with 
increasing wavelength. An excess of sources above the model can be attributed to the detection 
of extragalactic sources that are not present in the Wainscoat model. The majority of the 24 fim 
(MIPSl) sources are in excess of the Wainscoat model and are therefore probably extragalactic in 
origin. 



4. Candidate YSOs in Cepheus 

4.1. Identification of YSO Candidates 

Young stellar objects (YSOs) and background galaxies can have similar infrared colors, but 
they can be differentiated by the fainter apparent magnitudes of the galaxies. There are several 
published schemes for identifying YSO candidates from their Spitzer colors, each relyi ng on specific 



subsets of t he Spitzer wavelength coverage and the equivalent 2MASS fluxes (e.g., iHarvey et al 



2006 



20071 : lR.ebull et al.ll2007l : lAllen et al.ll2007l : ICutermuth et al.ll2008l ). These schemes are com- 
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monly calibrated against a known galaxy catalog. SWIRE (the Spitzer Wide-area Infrared Ex- 
tragalactic Survey, another of the Spitzer legacy programs) was specifically designed to observe 
the infrared extragalactic b ackground aloiig sigh t-lines that avoided as much foreground (Galactic) 



contamination as possible ([Lonsdale et al.ll2003l ). Thus, it provides an ideal calibration dataset to 



use with any galaxy-rejecting protocol. 

IRAC observations were targeted towards regions of high visual extinction with the expectation 
that these regions will harbor the highest concentration of YSOs. A side-effect of the different 
operational modes of IRAC and MIPS is that the MIPS data cover an area four times as large as 
the IRAC data. A galaxy-rejection scheme that used IRAC-only data would be unusable across 
3/4 of our surveyed area. Therefore we adopt different schemes for the common (IRAC-I-MIPS) 
and the MIPS only areas. A significant factor in both schemes is a 24 /xm MIPS detection. Since 
that band is not as sensitive as the IRAC bands, a third scheme is adopted for sources detected 
only with IRAC. 

Table [6] lists the YSO candidates identified towards the regions in our survey. The first three 
sections of Table[6]hst the YSO candidates identified by the 5-Band, IRAC-only, and 2MASS/MIPS 
schemes respectively (described below). The first column lists the source number. The second 
column lists the Spitzer catalog position identifier of each source (a contraction of the source's seven- 
digit right ascension and declination without delimiters). All positions are given in J2000. Column 
3 lists cross-identification names from the literature. Columns 4-9 list the fiux and statistical errors 
or where appropriate the band-filled upper limit. Column 10 lists a series of fiags denoting which 
candidate identification methods identify the source. The letters F, I and M are shown for sources 
identified by the 5-Band, IRAC and 2MASS/MIPS methods respectively. A source may satisfy 
the rule sets of more than one method, but it is assigned to a group based on the order F>I>M. 
This catalog includes data previously published in a series of c2d papers examining individual dark 
clouds. Column 11 lists references to these and other studies. 



4.1.1. 5-Band Identifications 



Harvey et al.l (120071 ) used a five-band (4.5, 8.0, and 24 fim, plus 2MASS H and Kg) scheme to 



estimate the probability that a given source is a background galaxy. The unnormalized probability 
Pgai is the product of a series of probabilities that a given source is a background galaxy. Each of 
these probabilities is based on the selection criteria from a different color-magnitude diagram. The 
final value of Pgai is moderated by a series of additional factors w hich include whether the source 
is extended at 3.6//m or 4.5//m (see Table 1 of lHarvev et al.l 120071 for a full list). Additionally, the 
scheme requires that a source be detected at each IRAC band and at MIPS 24 /xm, irrespective of 
whether that band was actually used to construct Pgai. 

Under the 5-band scheme, a catalog is filtered to reject those sources that can be adequately 
modeled by a stellar photosphere. Pgai is calculated for the remaining sources and these are then 
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filtered to retain those sources with a suitably low value of -Pgai, i-e., rejec ting sources that are 
statistically likely to be background galaxies. From their study of Serpens, iHarvev et alj (|2007l ) 
found that an upper limit of log(Pgai) < —1.47 rejected the galaxies from their SWIRE control 
catalog. This method creates a catalog of YSO candidates that is largely free from background 
galaxies. The catalog will be luminosity limited partially because it requires [24]>10 (fainter objects 
at [24] are preferentially background galaxies, but this cut off will also exclude faint YSOs). 

Figure [5] shows a histogram of log(Pgai) for sources in the Cepheus catalog with detections in 
all IRAC bands and MIPS 24 /im. The peak at -5 contains the majority of the YSO candidates 
identified by this method. There is also a small tail of sources to the left of the -1.47 divide (the 
vertical dashed line). There are 12 sources in the tail with —2.5 < log(Pgai) < —1.47 and 3 (25%) 
of them are previously known YSO candidates. We ch oose to retain the canonical dividiing lin e for 



YSO candidates for consistency with the c2d studies (jHarvey et al 



2007 



Evans et al.ll2008l ) and 



because moving the division by any significant amount would exclude the three known sources. If 
half of the tail sources were background galaxies, it would represent a contamination of ~ 5% in 
the total number of YSO candidates identified by the 5-Band method. 

Figure [6l show s a sequence o f color -color plots for sources detected in all four IRAC bands and 
shows the cut-offs iHarvey et al.l ()2007l ) used to construct Pgai . Dashed- lines show "fuzzy" limits 
while solid lines show hard limits. The left panel in each pair shows the SGBS data. Black filled- 
circles and crosses are respectively point-like and extended sources identified as YSO candidates 
by Pgai, whilst dark-gray circles are sources identified as point like (filled circles) and extended 
(open circles) galaxies. Also shown as pale-gray points are sources that were identified as stellar 
photospheres via their SEDs. The right-hand panel of each pair of plots in Figure [6] shows a contour 
plot of the SWIRE Catalog after it has been processed in a manner similar to the SGBS Catalog. 
It can be seen how the limits used in the calculation of Pgai have been chosen to reject regions of 
parameter space where there are significant numbers of SWIRE galaxies. These plots show that 
the galaxies are tightly clustered and that there is no significant overlap between the different 
candidate types. The 5-Band scheme identified 98 YSO candidates and 1 possible YSO candidate 
with a band- filled (upper limit) flux at 8 ^m (shown as an open square in Figured]). Sources that 
satisfy the 5-Band Scheme are listed with an F flag in column 10 of Table [H 



4-. 1.2. IRAC-only Identifications 



For those sources which have no detection at 24 /i m, we must us e selec tion rules based on IRAC 
colors alone. For this, we use the selection rules from lHarvev et al.l (|2008l ). They define ([4.5]- [8.0] 
< 0.5 and [8.0] > 13 — ([4.5] — [8.0])) to reject galaxies and stars. The upper row of Figure [7] shows 
a plot of [4.5]-[8.0] versus [8.0] for sources without a MIPS 24 //m detection and comparative plots 
of the 5-Band and SWIRE sources. Of the YSOs identified by the 5-Band scheme, 12 would be 
misidentified by the IRAC-only scheme as stars due to their small [4.5]- [8.0] color and six bright 
galaxies would have crept into the sample. 
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Sources that satisfy the IRAC-only scheme are hsted with an I flag in column 10 of Table [6j 
A source flagged with a lowercase f has a 24 ^um upper limit flux, but would otherwise satisfy the 
5-Band scheme. A comparison of sources with I and f flags shows that the IRAC-only scheme and 
5-Band scheme are consistent for sources without an 24 fj,m detection. A total of 10 additional 
YSO candidates was identified by the IRAC-only scheme. 

An additional constraint can be imposed in color-space as exter nal galaxies are often rich 



in PAH emission compared to YSO sources. iGutermuth et al.l (j2008l ) proposed using [4.5]- [8.0] 



vs. [3.6]-[5.8] and [5.8]-[8.0] vs. [4.5]-[5.8] color-color diagrams to reject background galaxies with 
strong PAH emission. Sources that match these rules are shown with triangle rather than circle 
markers in Figure [71 No IRAC YSO candidates have evidence of PAH emission from this technique, 
but two 5-Band YSO candidates are flagged. One of these is L1148 IRS. 



4.1.3. 2MA SS/MIPS Identifications 
An inherent strength of the 5-Band scheme is that it relies on a broad range of colo r coni - 



binations, but it can only be used in regions where there is IRAC coverage. iRebull et al.l (j2007l ) 
proposed a scheme to identify YSOs using MIPS 24 //m and 2MASS Kg and the selection rules 
Kg — [24] > 2 and Kg < 14. To this we add the limit [24] < 10 for consistency with the 5-Band 
scheme. The lower row of Figure [7] shows color-magnitude plots of Kg — [24] versus Kg for sources 
without an IRAC detection and comparative plots of 5-Band and SWIRE sources. Column 10 of 
Table [6] includes an M flag for sources that would have been identified as a YSO candidate by this 
method. A significant fraction of 2MASS/MIPS YSO candidates were found within 1 magnitude 
of the selection limits. These are marked as open circles on Figure [7] and listed with a lowercase m 
in Column 10. 

The 2MASS/MIPS scheme would have misidentified 19 of the 5-Band YSO candidates as 
background galaxies and stars. From Figure [71 we can see that the majority of these appear in the 
region of background galaxies with a Kg — [24] color of 10. This is approximately equivalent to 
the spectral index of an embedded protostellar source (a Class I YSO, see following section and 
Figure [8]) . This shows a further advantage of the 5-Band scheme in that it is more sensitive to 
embedded sources as it does not rely so heavily on Kg. Based on the MIPS scheme, 4-5 galaxies 
could have crept into the YSO sample. A total of 24 additional YSO candidates were identified by 
the 2MASS/MIPS scheme that were not identified by the 5-Band and IRAC-Only schemes. 



4.1.4- The Final Catalog 

The morphologies of all YSO candidates were checked at each wavelength listed in Table [3] to 
screen for obvious galaxy candidates or possible artifacts. During the visual checking, an unusual 
extended source was found at 21^02'"21.2^+68°04'36". It is a flat spectrum source across the 2- 
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24 fxm waveband, but rises at longer wavelengths. It is detected in all three MIPS wavebands. It 
was coincident with a spur in the extended 160 emission from the NGC 7023 nebula and a small 
i ncrease in visual ex tinction. It corresponds with a 1.4 GHz (20cm) source NVSS 210221+680436 



(|Condon et al.lll998l 1. Its colors of [5.8]-[8.0] = 2.8 and [3.6]-[4.5] = 0.0 place it well away from the 
color-color clustering of the other YSO candidates (see Section . This object has been excluded 
from further consideration as a YSO as it is probably a background galaxy. 

A total of 133 YSO candidates was identified by the three schemes. The vast majority of 
these (99) were identified by the 5-Band scheme and a further 10 were identified by the IRAC 
scheme. The 2MASS/MIPS data identified a further 24 candidates. Although the 2MASS/MIPS 
data covered a larger area than the 5-Band scheme, most of it was at a lower visual extinction. 
Thus, it could be argued that we would not expect to detect the same density of YSOs. The 
position of each YSO candidate was checked against the SIMBAD database. The SIMBAD name 
and other selected literature names for with each source are listed in column 3 of Table [H A total 
of 59 YSO candidates were found to have antecedent catalog names. 

It should be noted that not all known YSOs and YSO candidates coincident with our mapped 
area were identified by the three YSO candidate identification schemes. Section [6?2] below discusses 
the completeness of our YSO sample and Table [TBI lists 15 known sources missed by our selection 
criteria. These include luminous YS Os like HP 200775 (the dr iving source of the NGC 7023 nebula) 



and PV Cep (a "run-away" YSO, iGoodman Arcd l2004l ) which were excluded because they 



saturate the Spitzer detectors. In addition sources like XMMU J223727.7+ 751725 were excluded 
since they have no detectable infrared {2MASS / Spitzer) flux. From this list 10 are coincident 
with entries in the Cepheus catalog, but which were not identified by the three Spitzer schemes. 
Photometry for these 10 sources is listed in the fourth section of Table [6l We refer to these 10 
sources as non-Spitzer identified YSO candidates. 

Hereafter, we assume that these 143 (133 Spitzer plus 10 previously identified) YSO candidates 
are actually YSOs and analyze their properties accordingly. 



4.2. YSO Classification 



YSOs can be sepa rated into a series o f four evolutionary classes depending on either their 
infrared spectral index (jWilking et al.lll989l ) or the mean frequency of their SED (jMvers Ladd 
19931 ). The spectral index classification of a source as a Class II or Class III YSO is generally 
correlated wit h its respective class i fication as a Classical T T auri (CTTS) or Weak-Line T Tauri 
(WTTS) star (jwilking et al.lll989l : lAndre k Montmerlelll994l ). Class and Class I protostars are 
younger embedded YSOs and are differentiated from each other by the amount of material remaining 
in their envelopes. The envelope appears as a colder component to the SED and peaks towards the 
submillimeter. The original definition of a Class protostar was an object whose submillimeter 
luminosity (A > 350 ^m) contributed greater than 5% to the source's total bolometric luminosity, 
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which is equivalent to a protostar-to-envelope mass ratio of less than one (jAndre et al.lll993l ). The 
second method of YSO classification uses a bolometric temperature calculated for a blackbody 
which has the same mean frequency as the source's SED. The value of the bolometric temperature 
decreases from Class III to as the SED is increasingly dominated by long wavelength emission. 
The results for this method are described in Section 15. 2i 



Table [7] lists the derived properties of the YSOs from Tabled The first column lists the index 
number of the YSO. Column 2 lists the dark cloud from Figure [J with which the YS O candidate is 



associ ated. Locations that have a visual extinction less than one magnitude in the lDobashi et al 



(j2005l ) maps or are not part of an identifiable YSO group (see Section [6. 3p are listed as "off-cloud." 
Column 3 lists the spectral index of the source, aiR, as given by 

dlog(AS'A) , . 

aiR = —71 ^ (1 

a log A 

where Sx is the monochromatic flux density at wavelength A. The index was calculated by a 
least-squares fit to all available data in the range 2MASS (2.2 fim) to MIPS 24 fim. 



The mapping betw een the YSO evolutio nary sequence and spectra l index proposed bv lWilking et al 



(|l989l ) and amended bv I Greene et al.l (| 19941 ) and I Andre k Montmerld (| 19941 ) are Class I {ax < 0.3), 
"flat" (0.3 > ax> -0.3), Class II (-0.3 > ax> -1.6) and Class III (-1.6 > ax). The decreasing 
spectral index is attributed to the decline in the amount of circumstellar dust with advancing YSO 
evolution. The associated class for each source is listed in column 4 of Table [71 "Flat" spectrum 
sources have spectra intermediate between embedded and T Tauri stages. Class YSOs cannot be 
distinguished from Class I YSOs via this method as it is not sensi tive to the long wa velength part 
of the SED where the majority of the envelope emission radiates (lEnoch et al.l 120081 ). 



The majority of the data points used to calculate aiR are in the IRAC regime. For those YSO 
candidates without an IRAC detection, ajR is only calculated from 2MASS Kg and MIPS 24 //m 
- effectively the -fCs-[24] color. To test the robustness of using just the 2MASS/MIPS color, the 
top of Figure [8] shows ^^^-[24] versus aiR for the YSOs candidates detected with IRAC. The tight 
correlation follows a linear relationship thus proving that the equivalent color can be an acceptable 
proxy for a spectral index calculated using more intermediate points. 

The solid line through the points in Figure [8] shows the theoretical relationship between Kg- 
[24.0] and ax. For two generic bands, A and B, the color and spectral index are linearly related 
such that mA - m-B = Ciax + C2 where Ci = -2.5X, C2 = -2.5 + 2.51", X = log(AA) - log(AB) 
and Y = log(S'o,A) — log(5'o,B)- The wavelength of each band is Aa and Ab, the flux zero points 
are S'o,a and S'o,b and the color is (mA — rriB). Using the data from Table O the coefficients for 
irs-[24.0] are Ci = 2.63 and C2 = 7.55. A linear regression to the data, shown by the dashed line, 
gives coefficients of Ci = 2.64 ± 0.06 and C2 = 7.73 ± 0.06. 

The color-aiR, equation can be used to recompute the class boundaries in terms of a source's 
irs-[24.0] color such that they become -ftrs-[24.0] = 8.32, 6.74, and 3.32 for the boundaries between 
Class I/"flat", "flat" /Class II, and Class II/III respectively. These values agree with those used by 
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Lee et al. 


(2006) and 


Rebull et al. 


(2007) 



the dotted lines. 

The lower panel of Figure [8] shows a histogram of the YSO spectral indices. The black line 
histogram shows sources with IRAC coverage. The range of values of ax is -2.28 to 1.63. The gray 
histogram shows all sources and includes 2MASS/MIPS YSOs. Both distributions peak around -1 
and show an increase in the number of sources from Class I to Class II. The peak o f the Cepheus 



spect ral index histograms matc hes the peak of the distributions seen in Serpens (jHarvev et al 



20071 ). Lupus (jMerm et al.ll2008l) and IC 5146 (IHarvev et al.ll2008l ). The range of extreme spectral 
indices is less than in IC 5146 (jHarvev et al.ll2008l ). 



The number of sources classified in each Class via a\ and that number as a percentage of the 
Spitzer identified YSOs is listed in column 2 of Table El The row for Class sources is left blank 
for aiR, as they can only be distinguished by their bolometeric temperature. We can compare the 
relative fraction of aiR classes for the Ce pheus Flare against the values for other regions from the 
c2d survey (Figure 5 of lEvans et al.ll2008l ). The percentage of sources in each class closely matches 
the numbers for the entire c2d survey (shown in column 4 of Table [8]). Star formation in the 
Cepheus Flare is comprised of a series of small YSO groups, isolated cores, and a single (relatively) 
large YSO group (see Section [6]). That this mixture gives the Cepheus Flare a relative number of 
sources equal to the c2d survey may just be a coincidence or it could be that the mixture of star 
formation modes in Cepheus mirrors the balance of modes across the wider c2d sample. 

Of the individu al c2d regions, Cepheus most closely matches the oir clas s profile of Serpen s 
(jHarvev et al.l 120071 ). which is also very similar to that for the full c2d s urvev (lEvans et al.l 20081 ). 
Cepheus also has the same relative number of Class I sources as Serpens (jHarvev et al.ll2007l ). The 
relative number of Class I sources in Cepheus is ~ 10 percentage points higher than that seen in 
the Chamaeleon II and Lupus clouds although it i s ~5 percentage p oints lower that that seen in 



the IC5146 and Perseus clouds (jHarvey et al. 



2008 



Evans et al.ll2008l ). 



Kun et al.l (j2008l ) include in their review of the Cepheus star formation region a list of Classical 
and Weak-Lined T Tauri stars whose pre-main sequence nature has been confirmed by spectroscopic 
observations. There are 31 Spitzer identified YSOs from our sample that appear in that list (these 
have a reference to the Kun review in the last column of Table [6]) . Twenty-five of these sources 
(80%) are classified as Class II by their aju value. This match agrees with the usual interpretation 
of a YSO with a Class II infrared spectrum as most probably being a T Tauri protostar. 



There are 4 sources on the iKun et al.l list that are among the previously known YSOs listed 
in Table [16] that were coincident with our mapped area, but which were not identified as YSO 
candidates by the Spitzer schemes. Despite not being identified as YSOs 3 of these 4 so urces were 



coincident with entries in the Spitzer Cepheus catalog. Therefore, 34 out of 35 of the iKun et al 
stars had i nfrared emission that was detectable by Spitzer. It should be noted, however, that the 
Kun et al.l T Tauri list is from a project that is still underway and is not necessarily complete. 



Our Class III to Class II ratio is 1/8, whereas WTTS dominate CTTS in spatially complete 
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surveys. For example the W TTS/CTTS ratio in Taurus is 8 (jNeuhaeuser et al.lll995l ) and 3/2 in 
IC 348 (jLuhman et al.l 120031 ) . Infrared selection schemes will naturally be less sensitive towards 
objects with smaller infrared excesses, i.e., Class III YSOs. Estimates of completeness can only 
be made if there is a previous ce n sus o f WTTS to compare against. In a comparable study of 
the Lupus III region, iMerm et al.l (|2008l ) estimated their completeness for Class III sources with 
infrared excesses against objects with no infrared excess as ^ 50%. If the low rate of T Tauri 
sources in the Cepheus cores without detectable infrared excesses is real and similar to Lupus III, 
it would further reinforce the idea that the YSO population is comparatively young. 

Figure [7] shows that no 2MASS/MIPS YSO candidates are detected with a Ks-[2A] color 
greater than ~ 5 (or 6). From Figure [HI we see that this means that the 2MASS/MIPS method 
is significantly biased towards the detection of Class II and Class III YSOs. Deeper Kg data are 
needed to make a more complete survey of YSOs in the MlPS-only regions. 



4.3. Color-Color Diagrams 

Figure [9] shows color-color diagrams of the 143 Cepheus YSOs listed in Table El A source is 
included in a given plot if it has been detected at each of the wavelengths shown in that plot. Red 
markers show Class I sources, green markers show fia t spectrum sources, blu e markers show Class 



II sources and purple markers show Class III sources. iRobitaille et al.l (|2006l ) generated a database 
of 20,000 YSO models and used them to delineate a series of regions in color-color space that corre- 
spond to three different phases of YSO evolution. They termed these evolutionary phases "Stages" 
to differentiate them from the equivalent "Classes" (a purely observational parameter which de- 
scribes a source's infrared signature and may be influenced by reddening, e.g., as a result of varying 
disk orientation) . The approximate boundaries of the regions containing t he different evolutiona ry 



states are shown on each of the three panels in Figure [9] (see Figure 23 of IRobitaille et al.ll2006l ) 



Allen et al.l (|2004l ) compared the theoretical Spitzer colors of YSOs with data from four em- 
bedded clusters. The region they identified as corresponding to the approximate domain of Class 
II sources contains almost all of our Class II candidates as shown by the thick lined box on the 
first panel of Figure [H The small six-sided region coincident with the Allen box shows the region 
where most of the Stage II Robitaille models lie. Our Class II sources cluster around this region, 
but there are a number that lie just above it. These could be Class II sources that still retain some 
degree of reddening. The Class III sources are clustered around the zero point on each axis. Two 
large regions in the first panel are divided by a line that runs through the Stage II region and Class 
II box. The region above this line contains the majority of the Robitaille Stage I models whereas 
the region below it is where any Stage may be present. Our Class I sources agree quite well with 
the Stage I area. 

The middle-panel of Figure [9] shows IRAC and MIPS colors with the expected color-space 
limits of the three Robitaille Stages (Stage III, II, and I running left-to-right). The points show 
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the Cepheus YSO candidates with the same color code as the first panel. In this color-space, there 
is also excellent agreement between the Stage spaces and the Classes. The green points show flat 
spectrum sources. These do not corre s pond to a specific Stage, but are found along the Stage I/II 
boundary, as expected. I Evans et al.l (|2008l ) found a similar result for the entire c2d survey and 
showed that extinction correction of the individual fluxes only marginally improved the agreement 
between the Classes and the Stages. The Stage regions take account of differ ent angles of inclinatio n 
and this is partially why the Class I sources scatter over such a large area (jRobitaille et al.ll2006l ). 
The angle of inclination angle of a disk system to the observer is important as it can have a strong 



effect on the source's infrared signature and lead to the source being misclassified ([Crapsi et al 
20081 ). 



The last panel of Figure [9] shows the H — Kg versus J — H 2MASS colors of the Cepheus YSO 
candidates. The YSOs are not as separated as in the other two color-spaces, but there is a trend 
for the Class I and Flat sources to scatter away from the main locus. The box around this locus 
shows the expected domain of reddened stellar photospheres from the Robitaille models while the 
enclosed region to the right shows the region where any Stage evolutionary model can be present 
dRobitaille et ahlbood ). 



5. Cepheus YSO Properties 
5.1. Additional Photometry 



5.1.1. Guide Star Catalog 



A search was made of the HST Guide Star Catalog-II (GSC-II) for visual companions within 5 
arcseconds of the Spitzer/2MASS sources. The GSC-II is an all-sky compilation of astrometric and 
photometric data from a series of different catalogs including the Palomar Sky Survey-II (POSS-II), 
the Palomar Quick-V survey, and Tycho Catalogs (jMcLean et al I&OO). Of the 133 YSOs identified 
by Spitzer, 93 were also identified in the GSC-II. All of these had either "star" or "non-star" as 
their GSC-II spectral classification; none were classified as "galaxy." 

Table [9] lists, where available, the GSC-II derived Johnson-Cousins BVRI photometry for 
sources listed in Table [6l Column 1 lists the YSO candidates index from Tabled column 2 lists the 
GSC-II identifier of the associated source, and columns 3-6 list the equivalent flux photometry in the 
BVRI bands. The majority of the GSC-II identiflcations were from the POSS-II survey andhad^Bj, 



■Rf, a nd /n photometry. Rp and In are from the John son-Cousins photometr ic system (iReid et al 



199lh . Bj was converted to the same system using the lBlair &: Gilmord (Il982l ) color transforms. V 
data from the Palomar Quick-V survey were assumed to be in the standard Johnson system. The 
bright ~7 mag Bj^ and Vr magnitudes for HD 2 00775 (YSO #136) were from the Tycho survey and 



were converted to the Cousins-Johnson System (|ESA 



to fluxes using the Cousins photometric zero points (jBesselllll979l ). These zero points are listed in 



19971). All B VRI magnitudes were converted 
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column 4 of Table [3] and the transforms are summarized in Table [TUJ The errors on the optical 
photometry were typically 30% and could be as high as 50% in some cases. 



5.1.2. IRAS 



A search was also made for coinc ident sources froin the IRAS Faint Source Catalog within 
15 arcseconds of the Spitzer positions (jMoshir et al.lll990l ). This catalog was used as it gave more 
matches thari the Point Source Catalog and may be more reliable in regions of strong nebulosity 
([Rebull et al.ll2007l ) . The lower resolution of the IRAS survey meant that it was possible to associate 
multiple Spitzer sources with a single IRAS source. In these cases, the IRAS source was assigned 
to the Spitzer source whose MIPS 24 ^um flux was nearest to the IRAS source's 25 flux. No 
color corrections were applied to the IRAS fluxes. Column 7 of Table [9] lists the associated IRAS 
identifier and columns 8-11 list the 12.5, 25, 60, and 100 /um fluxes and upper limits. 

Twenty cross identifications (15%) were made between the 133 YSO candid ates and the IRAS 
faint source catalogue. Of these 9, are coincident with sources i n thelKunI ( 19981 ) list of IRAS-based 
YSO candidates. A further 2 sources are coincident with the iKunI (119981 ) list, but are not cross- 
identified with IRAS in Table [9] because we restricted ourselves to just the Faint Source Catalogue 
and did not include the Point Source Catalogue. The Kun IRAS Sources are listed with literature 
names beginning with K98b (Kun 98 Second Table) in Column 3 of Table [6j 



5.1.3. Suhmillimeter 

While more evolved protostars may be discerned by a negative spectral slope or a high bolo- 
metric temperature, it is more difficult to determine the evolutionary state of young YSOs without 
examining emission from their envelopes at longer wavelengths. Table [TT] lists available far-infrared 
and suhmillimeter photometry for Spitzer YSO candidates. The first column lists the YSO index 
and the second column lists the dark cloud associated with the YSO. Columns 3 to 6 list the actual 
fluxes quoted to a common precision. The last column lists, where relevant, a literature reference 
for the quoted flux. The 160 ^m fluxes were measured using a 40" radius aperture and 450, 850, and 
1200 ^m fluxes were measured in a 20 or 25" radius aperture depending on the quoted reference. 

Column 3 of Table [TT] lists MIPS 160 ^m fluxes. The fluxes were measured from the unfiltered 
BCD maps. Backgrounds were subtracted using a sky annuli of 1.0 to 1.87 times the radius of the 
aperture. Several of the apertures contained more than one YSO candidate. For these sources, 
the measured flux was divided evenly between the candidates within the aperture. Several bright 
sources, including PV Cep and L1157, were found to be completely saturated in the 160 /im BCD 
images. These sources and those within regions of large scale extended emission where flux cannot 
be reasonably assigned to a compact YSO candidate have been excluded. The MIPS 160 /im 
maps shown later are constructed from BCD images and have been filtered by a five-pixel diameter 
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median-filter to remove artifacts and to replace pixels excluded due to saturation. 



Columns 4 and 5 of Table [TT] list SCUBA 450 and 850 nm fluxes. The Submillimeter Common 
User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT) was a submil- 
limeter camera that could map the sky simultaneously at 850 and 450 ^m. Fully sampled maps 
were produced by either sca nning the hexagonal bolometer layout or offsetting the pointing-center 
in a 64-point jiggle pattern ([Holland et al.lll999l ). 



Several o f the c2d studies of the individual regions have existing SCUBA data for their sources. 
Young et al.l (120061 ) undertook a targeted SCUBA campaign to observe cores in the c2d program. 
The cores they observed included L1157, L1221, L1228, L1251, and L1155C. The YSOs in L1251B 
are tightly clustered , making a ssignm ent of the lower-resolution SCUBA flux difficult. In their 
analysis of L1251B, iLee et al.l (j2006l ) deconvolved the submillimeter emission of the embedded 
YSO candidates L1251B IRS 1,2 and 4 by using the brightnes s ratio of the resolved 2 4 //m sources 



Their estimates are listed in Table [TTl Data from lYoung et al.l (j2006l ) were used by the lYoung et al 
^200d ) study of LI 221. 



An additional search was made of the SCUBA archive for unpublished data coincident with 
other regions in this survey. Scan-maps of NGC 7023 taken on 1999 Oct 16-17 and jiggle maps of PV 
Cep taken on 1997 Oct 16 were do wnloaded. The data were re duced in the normal manner using the 
SCUBA User Reduction Facility (IJenness &: Lightfootlll998l ). The scan-maps were restored using 
the Emerson 2 technique ( Emerson Il995l ). The submillimeter zenith opacity at 850 and 450 /im was 
determined using the 'skydip' method and by comparison with the 1.3-mm sky opacity. Calibration 
was performed using observations of the planet Uranus taken during each shift. We estimate that the 
absolute calibration uncertainty is ±10% at 850 fini and ±25% at 450 //m, based on the consistency 
and reproducibility of the calibration. Secondary-beam corrected photometry was preformed in a 
40"-diameter aperture at the location of each YSO candidate. 



Kauffmann et al.l (|2008l ) surveyed a selection of the c2d srnall co res with the 1.2mm MAMBO 



bolometer array on the IRAM 30-m telescope (jKrevsa et al.l 1 19991 ) . Seven of the regions they 
surveyed are coincident with the Cepheus Flare cores that contain YSO candidates. Photometry 
for these YSOs is listed in column 6 of Table [TTl 



5.2. Spectral Energy Distributions 

Figures [lOl [HI [12] and \T3\ show plots of the spectral energy distributions (SEDs) of sources 
classified as Class I, Flat, Class II and Class III respectively. The open circles show, where available, 
the photometry from GSC-II, 2MASS, IRAC, MIPS, IRAS, and SCUBA. The arrows show the 
position of flux upper limits. The top-left hand corner of each plot is labeled with the index of 
the YSO. If the source index is followed by "xO.Ol" it indicates that the SED has been scaled 
downwards by two dex in order to place it on the same grid as the other SEDs. In Figures [12] 
and [T3l the gray lines show two comparison SEDs that have been normalized near the peak of the 
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dere ddened SEP (usually th e 2MASS J band). The solid gray hue is a NEXTGEN profile for a K7 



star ( Hauschildt et al 



19991) and the dashed gray line with error bars is the median SED for a T 



Tauri star in Taurus (jHartmann et al.ll2005l ). These SEDs are discussed in more detail in Section 
ET 



For the Class I YSOs in Figure [TOl with 3 or mor e detections longwards of 65 /im, we fit simple 
greybody SEDs assuming (3 = 2 (see lKirk et al.ll2007l for details) to characterize their submillimeter 
luminosities. Table [T2] lists the derived parameters from these fits. Column 1 lists the YSO index, 
column 2 lists the fitted dust temperature and column 3 lists the submillimeter luminosity integrated 
under the fitted greybody between 350 /^m and 2000 ^m. 

An alternatiye method to the spectral index for deriving the evolutionary classification of 
a YSO is to estimate the source's bolometric temperature, Tboi, which is the temperature of a 
blackbody that has the same mean frequency as a source's observed SED (|Mvers &: Laddlll993l ). 
Tboi can classify embedded protostars more effectively than the spectr al index as it uses the entire 
spectrum and not just the infrared portion used for the spectral index ([Enoch et al.ll2008l ). Column 
5 of Table [7] lists Tboi for each protostar calculated using all available data points. A Class source 
has Tboi < 70 K, C lass I has Thni = 70 - 650 K, Class II has Tboi = 650 - 2880 K and Class III 
has Tboi > 2880 K (IChen et al.l Il995l'l. The Thn i boun daries for the equivalent of the intermediate 



flat spectrum class were not set by 



Chen et al. 



(| 19951 ) , but lEvans et al.l (|2008l ) recently suggested 
boundari es of 350 K and 950 K. Column 6 of Table [71 lists the C lass for each YSO as derived 
from the lChen et al.l ([l995l ) boundaries with the I Evans et al.l ([20081 ) modification for flat spectrum 



sources. 



The number of YSOs in each Class as classified by Tboi is listed in Column 3 of Table [HI The 
spectral index cannot classify a source as a Class source so no total is listed the ajR column 
for such objects. While there is general agreement between the relative number of each Class of 
YSOs between both schemes (i.e., a few Class III sources compared to a lot of Class II sources), 
the bolometric method tends to skew a source's classification towards an earlier class (if the class 
changes at all). The spectral index is only calculated on flux data up to 24 ^m so flux data 
longwards of 24 /xm, as would be expected to dominate the SED of an embedded protostar, will 
not factor into a source's classification. For example, YSO #68 has a spectral index of -0.51 which 
classifies it as a Class II source whereas it has a bolometric temperature of 9.3 K which places it 
firmly in the Class regime. 



Enoch et al.l ([20081 ) examined the effect on Tboi of excluding a 160 fim data-point. They found 



that Tboi will be over-estimated and the ratio of Class to Class I sources will be skewed towards 
Class I sources. For the Cepheus data, we found that the exclusion of all data longwards of 150 /im 
(i.e., the 160 /xm and submillimeter data listed in Table [TT]) did not significantly affect the numbers 
of sources classified as Flat or Class II and III. It did, however, affect the relative number of 
Class and I sources. The number of Class and I sources classified by Tboi calculated using the 
submillimeter data is 18 and 9 respectively (as shown in Table [8]). Whereas, the number of Class 
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and I sources classifie d by Thni calculated without the submillimeter data is 4 and 16 respectively 
in agreement with the lEnoch et alj (|2008l ) finding. 

We retain the submillimeter data points when calculating the value of Tboi quoted in Table [7| 
and used to calculate the statistics in Table El The bolometric method undoubtedly gives a better 
assessment on the evolutionary status of an embedded protostar and we discuss the individual 
values for potential Class protostars in Appendix [Al The number of sources with high quality 
submillimeter data, however, is smaller than our total infrared sample size. Since this could bias 
our source classifications, we retain the oir classifications as our main scheme. This is also done 
for consistency with earlier studies. 



5.3. Luminosity Function 



The bolometric luminosities, Lboh of the YSOs are estimated by integrating under all available 
data-points in each YSO's SED using a simple trapezoidal method. The distance to each YSO 
was taken as the distance towards its respective associated dark cloud, as d i scusse d in Section [2j 
The resulting luminosities are listed in column 7 of Table [71 lYonekura et al.l (|l997l ) estimated the 
IRAS luminosity function of YSO sources associated with the molecular cores in their CO survey 
of Cepheus. They parametrized the number of sources dN in the luminosity interval L^, + dL^, as 



dN 
dL^ 



© 



(2) 



where A^o is a normalization factor and p is the power-law index of the function. lYonekura et al 
(jl997l ) fit Nq = 8.9 and p = 1.40ib0.32 above their completeness limit of IL© for their close group of 



IRAS sources. Figure [TH shows in its upper panel the function for the 133 Spitzer identified YSOs. 
The slopes were fitted by a least-squares fit above the estimated break of log(L^,/L0) = —1.5. 
The slope of the lumin osity function fitted to all sources is p = 1.61 zt 0.13. This value agrees 
with upper range of the lYonekura et al.l (jl997l ) error bars. The individual regions show a range of 
luminosities plus a single source that is approximately one dex more luminous that the rest of the 
YSO population. 

The lower panel of Figure [TH shows the same histogram data as in the upper panel but plotted 
with dN rather than log{dN / d{L / Lq)) on the abscissa. The luminosity distribution peaks at the 
point where the upper plot diverges from the single-power law, i.e., at a luminosity of 0.06 Lq. 
This Spitzer compl eteness limit is similar to that found for other c2d and G ould Belt regions (e.g ., 



Harvey et al.ll2007l ) and is 1.5 dex lower than the IRAS completeness limits (jYonekura et al.lll997l ) 



A lot of sources will have differing photometric coverage depending on DSS, SCUBA, and 
IRAS detections. Therefore, we also calculate the more homogeneous -Ljr, the luminosity integrated 
between 1-30 /xm for all sources with a 3.6 /Lim detection. The distribution of Lin is shown in Figure 
as the dashed histogram and the individual values are listed in column 8 of Table [71 There is 
good agreement in the peak position and width of the two histograms, but it is noticeable that 
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the LiR lacks -Lboi's higher luminosity tail. When compared to the equivalent Ljr distribution for 
Serpens and Lupus, the Cepheus Ljr distribution peaks in the same place, but the breadth of the 
peak is wider. All show the same peak just bef ore 0.1 Lq followed by a sharp drop to zero around 
0.01 Lq (jHarvev et ahlboOTl : Im^iu et ahlbood l. 



5.4. SED Modeling 



A dusty circumstellar disk and envelope at a temperature lower than the central YSO will 
contribute to the combined SED as excess emission at infrared and millimeter wavelengths. The 
magnitude of this infrared excess can be recovered if the SED of the star is subtracted from the 
combined spectra. The spectral type of the YSO needs to be known so that a template SED can 
be used, however. In general t he spectral type w ill not be known without a targeted campaign 
of spectrographic observations. iMerm et al.l (|2008l ) fitted NEXTGEN stellar models to continuum 



Harvev et al. 


( 


2007 


) and 


Harvev et al. 



underlying spectral type of their YSO population was a low-mass K7 T Tauri star or in the case 
of more luminous stars an AO type. We will use the same K7 assumption to calculate the infrared 
excesses for the Cepheus YSOs. 

For the Class II and Class III sources, the effect of reddening by line-of-sight extinction was 
removed from the SEDs by using the visual extinction towards each source calculated from its 
2MASS J — Kg color un der the assumptions of a K 7 underlying spectral type and a i? = 5.5 
interstellar extinction law (jWeingartner &: Draindl200ll ) . The dereddened data are shown in Figures 
12] and [13 as solid markers. A NEXTGEN K7 profile (IHauschildt et al.lll99d ) normalized against 
the dereddened J-Band flux is plotted for comparison. For a subset of sources, the dereddened 
optical photometry was significantly higher than the K7 profile that had been normalized to the 
dereddened J-Band flux. These sources include the known variable stars FT Cep, FU Cep, EH 
Cep and FV Cep. For these sources, it was necessary to modulate the estimated Ay down by a 
factor of 2-4 to make the K7 profile and the dereddened fluxes coincident. 

The dashed line in Figures [12] and [13] is the median SED for a T Tauri star in Taurus 
(jHartmann et al.l l2005l ) that has been normalized to the dereddened 2MASS J flux. The pro- 
file represents a prototypical optically thick accretion disk surrounding a T Tauri star. The SEDs 
of some YSO candidates (e.g. YSO #3) follow this profile quite closely, but others are closer to the 
K7 profile and only show a small amount of infrared excess at 24 //m (e.g. YSO #22). Approxi- 
mately 30 YSOs (25% of the disk population) have infrared excesses well below the median SED 
of the T Tauri stars in Taurus. Implying that these sources have evolved or settled disks. This 
low fraction also implies that the remaining disks in Cepheus show characteristics of being actively 
accreting and optically thick. 



Tables [13] and [H] list the SED modeling parameters and results for the Class II and Class HI 
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YSO candidates. Column 1 lists the YSO index. Column 2 lists the Ay that was used to deredden 
the SED and column 3 lists the waveband that was used to normalize the stellar profile. The 
luminosity of the YSO, Lgtar; was estimated by integrating the normalized stellar profile after it 
had been interpolated to the observed wavelengths. The luminosity of the circumstellar disk, Ldisk; 
was estimated by integrating the difference between the interpolated normalized stellar SED and 
the dereddened SED. Lgtar and log(Ldisk/-^star) are listed in columns 4 and 5 of Tables [T3l and [Til 

A histogram of log(Ldisk/-^star) is shown in Figure [15l This r atio can be used to charac- 



terize a disk as either accretion, passive reprocessing, or debris-like ([Kenvon &: Hartmann 



Hillenbrand et al.ll2008l ). The dividing- line between debris and passive disks su ggested by 



1987 



Currie 



(12008 ) and used by earlier c2d papers was log(Ldisk/-^^star) = —1-7, however the iHillenbrand et al 
(j2008l ) study of young debris disks s howed that their luminos ities were well below the -1.7 value. 
We use a division of -3 based on the IHillenbrand et al.l (j2008l ) result. These divisions are marked 
in Figure [I5l 

The dashed histogram shown in Figure [15] only includes sources with complete IRAC spec- 
troscopy (i.e., it excludes the 2MASS/MIPS identified sources). The difference between the two 
histograms shows that the majority of the sources identified as passive disks have been identified 
using the 2MASS/MIPS scheme. This pattern could be explained if Ljisk was underestimated due 
to the lack of IRAC photometry. Assuming that the estimate of Lgtar did not change, the addition of 
IRAC photometry could cause Ldisk to increase and would preferentially shift a source towards the 
accreting disk region. Thus, the observed frequency of passive disks could be lower than actually 
shown. 

It is clear that the majority of the circumstellar disks modeled by this method are accretion 
disks with a peak in the distribu tion of log(Lfiisk/L star) = 0.4. This is in agreement with the disk 
fractional luminosity for Lupus (jMerm et al.l 120081 ) and Chamaeleon (lAlcala et al.l 120081'). We d o 
not, however, find as many passive and debris disks as was found in Serpens (jHarvey et al.ll2007l ). 
Differences in stellar ages and masses between clouds could account for different distributions, but 
a detailed analysis of this possibility is beyond the scope of this paper. The large percentage of 
accreting disks is consistent with our finding that a relatively small number of sources have infrared 
excesses below the median T Tauri SED in Taurus. Together these points suggest that the disk 
average evolutionary status in Cepheus is very close to primordial and that there is little evidence 
for disk evolution. 



Cieza et al.l ((20071) introduced the second order parameters Aturnoflf and Oexccss to characterize 



the disk infrared excess. The first parameter, Aturnoflf) is the longest wavelength without significant 
infrared excess. We compute this as the band shortward of the last band where the ratio of disk flux 
to observed flux is greater than 80%. When no band has a ratio greater than 80% we set Aturnoff 
equal to the longest wavelength band. The second parameter, Oexcess, is the spectral index of data 
points longwards and inclusive of Aturnoflf computed in the same manner as oir. The calculated 
values of Aturnoflf and aexcess for the Cepheus YSOs are listed in columns 6 and 7 of Tables [13] and 
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[TH Figure [16] shows a plot of Aturnoff versus Oexcess • 
Class III YSOs are shown as asterisks. 



Class II YSOs are shown as open circles and 



Cieza et al.l (j2007l ) found an evolutionary sequence along Aexcess with Classical T Tauri stars 



having a Aexcess shortward of 2]VIASS Ks- They also showed that the majority of o^excess 

values 

cluster around -1 irrespective of Aturnoflf, but that the spread in aexcess increased with Aturnoff- 
Figure [T6l shows the same trend for clustering around aexcess ~ —1; but as most stars have had 
their stellar profiles normalized to the 2MASS J we would not expect to see values of Aturnoff equal 
to it. Surprisingly, there are no objects with Aturnoff = 4.5 /im. Xhe calculation of aexcess 

and 

Aturnoff is affected by the normalization of the stellar profile and to the assumed spectral type. 
Thus, the empty 4.5 /im band and the scatter of excesses at shorter wavelengths could possibly be 
due to imperfect knowledge of each source's spectral type. 



Harvev et al.l (|2007l ) plotted a similar diagram for Serpens and showed that the Serpens Class 
III sources predominately had a Aturnoff ^ 5.8 yum. We find a similar result for the Class III sources 
in Cepheus as shown by the values of Aturnoff in Column 6 of Table [U] and the few sources that 
appear at IRAC 3 and 4 in Figure [16l We are unable to calculate aexcess for the three sources 
with Aturnoff = 24 fj,m as none of these sources have longer wavelength data points against which to 
calculate an index. All three were Class III sources, however, and would have followed the trend 
for those sources appearing to the right in Figure [TBI 



Merm et al.l (|2008l ) showed that an object in the top right of the diagram was likely to have 
a transitional disk - an optically thick disk with a central cavity larger than several AU. These 
objects appear as a photospheric SED with an infrared excess only at the longest wavelengths. 
The object in the top right of Figure [16] is YSO #83, a Class II YSO in L1251B associated with 
IRAS F22367+7448. It has been detected in x-rays and was clas sified as a Classical T Tauri based 



on its large Ha equivalent width and spectroscopic follow-up (jSimonl l2006l : iKun &: Prusti 



1993 



Kun et al.ll2008l ). A single transitional disk source may also be indicative that most of the disks in 
this region are at an early evolutionary stage. By comparison the Spitzer st udy of three clusters 



around Cep 0B2 found that 10% of detected disks were transitional in nature (ISicilia-Aguilar et al 
20061 ). Several sources have strong 70 u m excesses, e . g. Y SO #16, 35, 122. Similar objects are 



discussed in the Spitzer studies of Lupus (jMerm et al.l 120081) and Chamaeleon (lAlcala et al.l 120081 ) 
and may represent a very young population of debris disks ([Hillenbrand et al.l 120081 ). 



In Cepheus, there is a population of isolated T Tauri stars (Z ~ 117 — 122°) which are unas- 
sociated with molecular material and are outside of the regions surveyed in this paper {I < 115°). 
The isolation of these objects compared to their evolutionary status suggests that they formed in 
situ and that their na tal clouds have been removed by interaction with the Cepheus Flare Shell 
([Tachihara et al.|[2005l ). If cloud disruption is as effective as argued bv [Tachihara et al.l ([20051 ) then 
it could explain why we are preferentially seeing young disks in our YSO sample. Our observations 
were specifically targeted towards dense clouds and would therefore be missing older YSOs whose 
clouds have already been dispersed. In Appendix IA.4[ comparison of the YSO distribution to the 
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pattern of extinction and 160 /im emission shows that the L1228 South YSO group is on the very 
edge of the L1228 cloud. It is perhaps reasonable to suggest that the L1228 cloud is currently 
being disrupted and that we are witnessing the unveiling of a cluster of T Tauri stars similar to the 
isolated group observed by Tachihara. 



6. Extended Structure and YSO Distribution 

Figure [T7] shows the distribution of YSOs and visual extinction towards the regions in the 
Cepheus survey. Figure [T2] shows the pattern of MIPS 160 //m emission towards exactly the same 
regions. All regions are shown at the same angular scale. 



6.1. Comparison of Ay and 160 maps 



Figure [T7] shows Ay calculated from the Spitzer catalog. The line-of-sight extinction towards 
sources classified as reddened stellar photospheres between wavelengths of 2 and 24 /um was calcu- 
lated assuming an extinction law of R = 5.5. These points were then filtered based on the equiv- 
alent 2MASS extinction to correct for contamination from isolated clumps of extinction. These 
irregul arly spaced, filtere d-data were then averaged onto a regularly spaced grid using a Gaussian 
kernel (jEvans et al.l 120071 ). The resolution of the maps depends on the surface density of available 
sources. It was found that the highest resolution consistent between all regions in Cepheus was 150 
arcseconds. 

The approximate areas where we have IRAC and MIPS overlap are shown by the irregular 
boxes on Figure [TTl These are the regions where there is enough photometric coverage to calculate 
the Spitzer Ay. The greyscale within the overlap regions show the local range of the 2.5 arcminute 
resolution Spitzer Ay on a linear stretch. Outside of the overlap area, the greys cale shows the 



6 arc minute resolution DSS extinction on a linear stretch between 0.5 and 3 mag (jPobashi et al 



20051 ). Dashed black contours outline the region where the DSS extinction is greater than Ay = 1 
and the Spitzer extinction is greater than Ay = 5. The Spitzer extinction maps show higher values 
of Ay than the Dobashi maps because the infrared stars are visible at higher column densities. 

Figure [18] shows 160 fim emission observed by MIPS towards towards the five regions in Figure 
[TTl The greyscale is shown on a log stretch between the local minimum (approximately 20 MJy/sr) 
and 316 MJy/sr for NGC 7023 or 100 MJy/sr for the rest. The Ay contours and overlap box from 
Figure [17] are repeated for reference. There is a strong correlation between the distribution of the 
visual extinction and the 160 //m emission. This pat tern of the 160 //in emission following th e Ay 
map has also been seen in many other regions (e.g. iRebull et al.ll2007l : [Chapman et al.ll2007l ). 



In general, the 160 fim emission above ~ 40 MJy/sr is confined within the Ay = 1 DSS 
contour. Some regions, particularly L1241 and L1247, also show a correlation at smaller angular 
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scales. There are also subtle differences, however. For example, in L1228S the extinction peaks 
approximately 4 arcminutes further north than the 160 //m emission. The center of NGC 7023 
saturates the MIPS detectors and is left blank in the map. 

Table [TSl lists parameters for the clouds and YSO clusters/groups shown in Figure [T71 We 
identify dense cores based on a comparison of the the extinction and 160 /im maps. The Ay = 5 
contour was chosen to define the cores as it was the level where most of them separated from each 
other. These are labeled on Figures [T71 and [T8l Due to the partial mapping of L1147 and small size 
of L1155E these regions are merged with L1148 and L1155C when calculating derived properties. 
L1174A and L1174B with the unlabeled L1174C are hsted together as NGC 7023. For each of the 
remaining fourteen cores, we measured the peak and mean within the Ay = 5 contour, the area 
co ntained within that contour, and the equivalent mass calculated from the Ay/Nn relationship 
of lBohlin et al.l (jl978l ). Th e c2d and Gould B elt Ay values are calculated assuming an extinction 
law of i? = 5.5. To use the lBohlin et al.l (jl978l ) relation of 



Ay_ 



1.1 X 10 



-21 —2 

cm mag, 



(3) 



which was calculated for a value of = 3.1, we include a correction of 0.716. These values are listed 
in columns 11-14 of Table [T5l The equivalent values for the YSO clusters/groups are calculated 
using a YSO density contour that defines the edge of the group (see Section . 



6.2. YSO distribution and completeness 

The distribution of YSO candidates is over-plotted on the five regions in Figure [T71 The 
markers are color coded to the class of the candidate - red for Class O/I, green for Flat spectrum, 
blue for Class II, and purple for Class III. Columns 2-6 of Table [TSl list the number and type of YSO 
candidates that were detected towards each core within the Ay = 1 contour or YSO group/cluster 
boundary. Column 7 of Table [TS] lists the ratio of Class I YSO (column 2) to the total number of 
YSOs (column 6). In general, the protostellar cores had such ratios in excess of ~ 20% while the 
cores containing clusters had values less than ~ 20%. 

A YSO was considered to be "off-cloud" if it was outside of the Ay = 1 contour and was not 
a member of a formal YSO cluster or group (see Section [613]) . Results for these YSO candidates 
are listed as a separate group at the bottom of Table [151 The isolated population of Weak-Line T 
Tauri YSOs in Ce pheus have higher YSO -to-cloud separation ('^ 10 pc) than found in other star 



formation regions (iTachihara et al.[ [20051 ). While our observations do not get that far "off-cloud", 
we can use the number of off-cloud YSOs, most of which are Class II and Class III sources, to 
estimate the density of off-cloud YSOs that are serendipitously located close to dense regions. The 
area at an extinction of Ay > 1 is 2.14 square degrees and the total mapped area is 8.33 square 
degrees. We detect 19 Spitzer identified YSOs at an extinction of Ay < 1. Assuming that these 
are not background galaxies, we estimate an off-cloud YSO density of 19/(8.33 — 2.14) = 3.1 YSOs 
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per square degree. Conversely, we detect 133 — 19 = 114 Spitzer identified YSOs at an extinction 
of Ay > 1. This gives an on-cloud YSO surface density of 114/2.14 = 53 YSOs per square degree. 

The SIMBAD object types (the "otype" field) of the sources listed in the first three sections 
of Table [6] include YSOs and YSO candidates (Y*0 and Y*?), emission-line stars (EM*), variable 
star of various classes (V*, Pu*, Mi*, Or*), reflection nebulae (RNe), and T Tauri stars (TT*). 
A number of sources have multiple object types. A search was made for additional sources with 
similar object types in a region coincident with our mapped area. The result included 14 variable 
stars located around NGC 7023 and towards the L114 8-I-L1155 r ir ig. O nly HD 200775, however, 
had a match against a source in the Spitzer catalog. iKun et al.l (|2008l ) reviewed surveys of Hq 
emission line stars towards the Cepheus Flare and produced a list of emission line stars that had 
been confirmed spectroscopically as T Tauri stars. There are 19 emission-line stars in the list of 
additional SIMBAD sources. The emissi on-line stars were pre ferentially found to wards LI 228 as thi s 
area was the subject of a Ha survey by lOgura Satd (jl99d ). Thirty-one of the iKun et al.l (120081 ) 
T Tauri stars are coincident with Spitzer YSO candidates, but none of the additional SIMBAD 
emission line stars are confirmed T Tauri star. 

The SIMBAD search was narrowed to only include objects previously identified as YSOs and 
pre-main sequence stars or candidates for either. A total of 47 YSOs were found in the SIM- 
BAD database of which 36 were cross-identified with Spitz er YSO candidat es. Half of these cross- 
identified YSOs are in the L1251B region and are from the lLee et all \20()4 ) c2d paper. Excluding 
these, 18 out of 30 SIMBAD YSOs were cross-identified with the Spitzer YSO candidat es. The 12 
SIMB AD YSOs that were not identified by the 5pifeer schemes plus 3 T Tauri stars from lxun et al.l 
([200a) that were not identified by the Spitzer schemes are li sted in Table 1161 Column 4 lists the 
catalog that the object was taken from (S for SIMBAD, K for lKun et al.l 120081 ). Column 5 lists the 
core or region coincident with the object's position. Column 6 lists whether the object position was 
coincident with an entry in the Cepheus Spitzer catalog even if that object was not classified as a 
YSO candidate or formally detected. Ten of the additional YSOs were coincident with entries in 
the SGBS Cepheus catalog (PV Cep had only 2MASS data). Each of these non-Spitzer identified 
YSO candidates has been assigned a YSO Id (as listed in column 7) and their photometry is listed 
in the fourth section of Table [6l 

The embedded Class protostar L1157 was outside of our IRAC area and has no 2MASS 
detection. L1157 was studied with Spitzer bv iLoonev et al.l (|2007l ). They presented striking 8 /im 
images of LI 157 that showed its CO outflow in emission and a flattened circumstellar envelope in 
absorption against the bright 8 ^m background. Its IRAC fluxes have been added to our catalog 
based on 4.5-arcsec radius aperture-photometry conducted on post-BCD data downloaded from the 
Spitzer data archive. Spitzer data for PV Cep exists in the data archiv^, but it is saturated at all 
four IRAC wavelengths. The 2MASS/MIPS YSO #118 is also in the PV Cep IRAC field so we 
add its IRAC fluxes to the catalog in the same manner as for L1157. 



^AORNo. 18955008 
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HD 200775 would have been picked up by the IRAC YSO scheme if we had reduced the S/N 
cut to 2. It was detected in the filters used for the 5-Band scheme, but it was excluded because 
the 3.6 /im flux was below 2 a. The source is very luminous (Lboi = 100 Lq) and is saturated in 
all MIPS bands. It is illuminating the surrounding gas gi ving it a. large r flux uncertainty than an 
isolated star wo uld have. T he position of NGC 7023 S T dSellgrenlllQSa l was identified with IRAS 
F21023+6754 bv lXunl dlQQsl ^. Kun also identified the sources IRAS F21025+6801 and F20597+6800 
as IRAS YSO candidates within the NGC 7023 region. Of these latter three sources, only IRAS 
F21025+6801 is coincid ent with an entry in th e Cepheus catalog and then it only has a 24 ^um 
upper limit. None of the iTachihara et al.l (|2005l ) x-ray T Tauri candidates were coincident with the 
area we surveyed. 

A total of 143 YSOs were identified from the Spitzer catalog and the SIMBAD database. Of 
these 133 (93%) were identified by the color-color schemes discussed in Section 14.11 and a further 
10 (7%) were coincident with entries in the catalogue. The number of previously known YSOs that 
were not identified by Spitzer, but had entries in the Cepheus catalogue is listed in Column 8 of 
Table [T5l The total number of known YSOs that were not identified by Spitzer is listed in Column 
9 of Table [151 The number of Spitzer identified YSO candidates within Ay = 5 contours discussed 
in the previous section is Column 10 of Table [T5j 



6.3. Pre-stellar cores, protostellar cores, and cores with YSO groups/clusters 

The dense cores seen in extinction in Cepheus can be split into three groups. The first division 
are pre-stellar cores. These are cores that are believed to be taking part in the star formation 
process but have not yet formed an embedded YSO (i.e., Nyso = 0). We then define the cores 
that have A'ygo > within the Ay = 5 contour as protostellar cores. Concentrations of YSOs 
i n protostellar cores with at least 5 and 35 members are termed groups and clusters respectively 



(|Lada fc Ladall2003l ) 



For pre-stellar cores, the YSO count within the Ay = 5 contour is used because the area 
subtended by the Ay = 1 contour can be large enough to give a significant chance of including a 
field YSO. For example, the area of L1241 at > 1 is approximately 0.6 square degrees. Based 
on the average off-cloud YSO surface density, we would expect 1-2 YSOs to be coincident with this 
region. There is one Class III YSO just inside the Ay = 1 border of L1241. Given its class and 
the likelihood of finding a coincident YSO within that border, we conclude that L1241 is starless. 
The parameters for the Cepheus pre-stellar cores are listed in the first section of Table [T5l The 
respective masses and areas are calculated from the Ay = 5 contour. 

The cores that have A'ygo > are termed protostellar cores. The distribution of YSOs 
in Figure [17] shows a number of protostellar cores. There have been multiple metrics devised 
for classifying the composition and concentration of YSOs in protostellar cores. For example. 



Cartwright &: WhitworthI (|2004l ) used minimal spanning trees to derive a quantity Q that could 
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differentiate between a small, large-scale density gradient and multi-scale sub-clustering. Alter- 
natively, the spatial density of sources can be estimated by calculating the angular distance tn 
to the Nth nearest neighbor for each poin t on a regularly spaced grid ([Christopher et al.l 119981 : 
Gladwin et allll999l : ICutermuth et al.lboosi '). 



For consistency with previous c2d/SGB S papers, we follow the nearest-neighbor clustering 
scheme described bv lGutermuth et alj (|2005l ). When calculating tn, the appropriate correction for 
spherical coordinates was used (this is vital for Cepheus due to its high declination). This angular 
distance is then used to calcula te a volume nuniber d ensity assui ning a spherically symm etric 
distribution of cluster members (jGutermuth et al.l |2005| ) . Following iGutermuth et al.l (j2005l ) , the 
volume density of stars is given by 



N -1 



(4) 



where M^, is the average stellar mass and N is the chosen index for the neighbor. We use N = 5 and 
assume an I MF with an average stellar mass of 0.5 M fD for consistency with other c2d and SGBS 
papers (e.g. |j0rgensen et al.l 120081 : iHarvev et al.ll2008l ). The value of /O* was calculated across the 
regions in our survey. A limitation to this technique is that it can only calculate for groups with 
5 or more members. 

The value of above which a clust e r is stable against disruption from passing interstellar 
clouds is > 1.0 M0 pc~^ (ISpitzeii 119581 : iLada &: Ladal 120031 ). Cores with peak values of p* less 
than the disruption density or with fewer than N members are listed as purely protostellar and are 
listed in the middle section of Table [T5l 

The cores with peak values of p^, greater than the disruption density are listed in the third 
section of Table [T5l Five regions had stellar mass volume densities in excess of the 1.0 pc~^ 
disruption density. These are NGC 7023/L1172, L1251A, L1251B, L1228S, and L1228N. For this 
analysis, we have only included the YSO candidates towards NGC 7023 and not the main-sequence 
members listed in the literature. The pattern of p* towards each of these groups is shown by the 
colored contours on the upper panels in Figure [T71 The blue contour traces the disruption density 
while green contours trace values of 0.125, 0.25, 0.5, 2.0, and 4.0 times the disruption density. 

Figure [17] shows that the disruption contours are approximatel y elliptical, likely becau se these 
relatively small and isolated YSO groups are themselves elongated. Ij0rgensen et al.l (120081 ) defined 
an empirical limit of 25 times the disruption density as a level which appears to more adequately 
trace the local YSO distribution. This level is shown as a yellow contour in Figure [T71 The red 
contour over NGC 7023 shows the equivalent of the yellow contour after known non-YSO cluster 
members are included. The position of the 21 non-YSO cluster members is shown by yellow crosses. 

Each group can be classified by the number of members within a particular contour. YSOs 
within a minimum contour equal to the disruption contour (blue contour) are described as a "loose" 
association and YSOs within a minimum contour equal to 25 times the disruption density (yellow 
contour) are described as a "tight" association (jjcirgensen et al.ll2008l ). A cluster is defined as 
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an association with at least 35 members (jLada &: Ladal l2003l ). Based on these definitions, no 



association in Cepheus can be considered a true cluster although NGC 7023 is close. It is a loose 
association of 32 members which contains two tight groups, A and B, with 20 and 4 members each. 
The other four cores contain individual YSO groups that have tight and loose members, although 
the tight regions do not fragment as in NGC 7023. 

We have defined four different contours that can be used to define the extent of the cores. 
These are the Ay = 1 and Ay = 5 contours used for the cores themselves and the /)* = 1 and 
25Mq pc~^ used for the YSO groups /clusters. The third section of Table [15] lists the YSO statistics 
and physical parameters for the cores with YSO groups. Each core is listed three times, first for the 
core itself as defined by the Ay = 1 contour, then for the loose and tight YSO groups as defined by 
the blue and green contours in Figure [TTl The relevant contour is used for the calculation of the 
group/cluster statistics. The tight groups for NGC 7023 A and B are listed separately. 

In general, the statistics for the entire core and for the loose groups are approximately equal, 
but there are a few noticeable differences. The statistics calculated from the "loose" contour 
characterize the cores with groups better than the Ay = 5 contour. The NGC 7023 loose cluster 
subtends an area that is significantly larger than the area of the core as measured by its Ay = 5 
contour. The L1228 YSO group straddles both Ay contours so using those contours to count the 
YSOs gives an artificially low number. 

The final column on Table [15] lists the star formation efficiency (SFE) for each core as given 

by, 

SFE = ^''""'l^ (5) 

^cloud + Alstars 



where Mdoud is the mass of the gas and dust in the core (as listed in column 14 of Table [TS]) 
calculated from the mean Ay within the core's boundary contour. The mass of the stars, Mstars) is 
given by A'yso-^* where M^, is the same average stellar mass as assumed for the calculation of p*. 

There is no clear trend between the masses of the cores and the number of YSOs they contain. 
Thus, a larger stellar content automatically translates into a higher star formation efficiency. The 
protostellar cores without YSO groups/clusters have SFEs ~ 1% while the overall SFE for the 
cores with YSO groups/clusters is ~ 8%. The SFE for the tight component of the groups/clusters 
is > 3x higher than for the loose component. This can be interpreted as either an increase in the 
in situ star formation within the dense parts of the core or a drift of YSOs into the central part 
of the core from other parts of the cloud. In the former, the SFE is real and, in the latter, it is 
over-estimated because the mass of gas that those stars formed from is under-estimated. 

The bottom line of Table [15] shows the aggregated statistics for the entire mapped region. 
Columns 2-5 reiterate the YSO class statistics from col umn 2 of Tab l e [8] T able [TTl lists additional 



star formation statistics in the same format as table 3 of lEvans et al.l (|2008l ). Columns 1-4 of Table 
[T7]list the boundary used for the analysis, the number of YSOs within that boundary AVso , the 
solid angle of the region Q, and the integrated area of the region A assuming an average distance 
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of 300 pc. Columns 5 and 6 respectively list the number of YSOs per square degree and square 
parsec. Column 7 lists the star formation rate (SFR= A'YSO-^*/'^sf )• Where we have used the 
aforementioned assumption that = 0.5 M f?^ (the same used i n the clustering analysis) and that 
the timescale for star formation is Tgf = 2 Myr (jEvans et al.ll2008l ) . Column 8 lists the star formation 
rate per square parsec (SFR/^). Column 9 and 10 list the total mass of the cloud (Mdoud) a-iid 
the surface density of the cloud (Mdoud/^)- 



7. Discussion: Star Formation in the Cepheus Flare 



Star formation in Cepheus is mostly concentrated in small, isolated groups of YSOs. A total 
of 78 out of 133 S^piteer-identified YSOs, i.e., 59%, were determined to be members of loose YSO 
groups, leaving 41% of the YSOs distribute d elsewhere. By co mparison, the average distributed 
YSO fraction in c2d clouds was only ~ 10% (jEvans et al.l 120081 ). Figure [T9l shows the distribution 
of YSO and YSO candidates towa rds the Cepheus Flare plotted over the distribution of Ha (log 
scaled greyscale; lFinkbeineiil2003l) and integrated C O emission (black contours; the same data as 



the white contours from Figure [H I Dame et al.l l200ll ) . The black foot prints show the limit of our 
survey and the red crosses show the positions of Spitzer YSO candidates from this p aper. Th e 



posit ions of spectroscopically confirmed T Tauri stars are shown by open blue squares (IKuii et al 



20081 ) and the position of X-ray WTTS are shown by filled blue squares (jTachihara et al 



200 



i 



The green circles show an unvetted list of objects described in the SIMBAD database as YSOs 
irrespective of their survey of origin. 

The masses of cores with YSO groups (A'yso ^ 5) and cores with only a few YSOs (A'^yso < 5) 
were found to be broadly similar, which suggests a factor other than core mass determines how 
many stars form. If we examine the location of the YSO groups (L1251A, L1251B, L1228S, and 
L1228N) in Figure [T9l we see that they are all located in the region nearest the Cepheus Flare Shell 
(CFS), whereas the majority of the cores with less than 5 YSOs (L1148, L1157, etc) are located on 
the opposite site of the region. There appears to be two different modes of star formation occurring 
here. The first mode of groups of YSOs has a SFE of ~ 8% and has produced a YSO population 
that is dominated by Class II sources. The second unperturbed, almost quiescent mode has a SFE 
of ~ 1% and has produced YSOs with a fairly uniform spread of infrared classes. 

There is mounting ev idence that star for i nation in the Cepheu s Flare is trigger e d /influenced 
by the passage of the CFS (|Grenier et al.lll98d : [Tachihara et ahlbood i lKun et aklboosl ). IkuuI (Il998l ) 
noted that the surface distribution of IRAS identified YSOs in the Cepheus Flare peaked towards 
the edges of clouds and inferred that their formation was triggered by external shocks. Likewise, 
Figure [19] shows that the locations of the small YSO groups in L1228 and L1251 are on the edge 
of the CO contours facing the Cepheus Void. The cometary shape of C^^O emission in L1251 and 
an unusual velocity gradient across the cloud appears to have been caused by interaction with a 
passing wind/shock wave, but the failure to de tect SiO einission from the cloud may indicate that 
the shock phase, at least in L1251, has passed (ISato et al.lll994l ). It is the cometary head of L1251 



- 32 - 



that contains the L1251A and B YSO groups. 



The L1228S YSO group not only Ues on the edge of the CO contours but, as shown in Figure 
[T71 it also lies on the edge of the DSS Ay = 1 and Spitzer Ay = 5 extinction contours. The peak 
of the L1228S Spitzer Ay map is offset from the peak of the 160 emission. The difference in 
position could be due to heating of the dust caused as the core is dispersed by interaction with the 
CFS. There are small collections of T Tauri stars (open blue boxes) on the Void side of both the 
L1228 and L1251 groups, suggesting that they have become exposed as material has been stripped 
away from t he lea ding edge of the molecular clouds. For the group of T Tauri stars near L1251, 



Toth et al.l (|l995l ) estimated that it would have taken 0.1 Myr after the passage of the shock for 
an offset of 1° to form, assuming an expansion velocity of 10 kms~^. 

There is a collection of SIMBAD YSOs and T Tauri stars close to the end of the CO filament at 
/ = 112°, b = +14°. It is outside of our selection criteria (it has peak Ay < 3), but it is coincident 
with the L1235 dark nebula. Given the position of the YSOs, the CO contour, and its orientation to 
the Cepheus Void it is possible that L1235 is another small YSO group, but one that is intermediate 
in evolution between L1228S and the Cepheus Void WTTS Group. Two other concentrations of 
YSOs are shown in FigurellQl L 1199 at a distance of 500 pc and a group associated with L1217 at 
a dist ance of 400 pc ( Kiiij[l998|). Both L1217 and L1235 are coincident with the G109+11 infrared 
loop (jKiss et al.ll2006l ). The small, isolated protostellar core LI 221 was included in our survey and 
is also coincident with the interior of the infrared loop, but its assumed distance of 250 pc places 
it foreground of L1217 and L1235. 



The age of the supernova at the center of the Cepheus Flare Shell is 40,000 yr (jGrenier et al 



19891 ). but the bubble that the shell encloses appears to be significantly older, 7 Myr, and was prob- 
ably a pre-existing wind driven cavity created by a star that subsequently exploded as a supernova 



(jOlano et al.l l2006l . and discussion therein) . Ionization pressure from H II regions can enhance 



gravi tational collapse in surrounding material, e.g., the Horsehead Nebula (IWard-Thompson et al 



20061 ). and trigger star fo rmation in the heads of pillar-like filaments, e.g., M16's "Pillars of Cre- 
ation" (jHester et al.lll996l ). We conjecture that this mode of star formation may have been respon- 
sible for some of the WTTS curren tly currently situa ted inside the Cepheus Void. The Cepheus 
Flare Shell is an expanding H I ring (|01ano et al.ll2006l ). however, and shows negligible Ha emission 
when compared with the Ceph 0B2 H II region (see Figure [T9|) . The lack of Ha emission means 
that any hot gas associated with the CFS has had sufficient time to cool or escape. 



Following the lOlano et al.l (|2006l ) analysis of the dynamics of the CFS we estimate that its 
radius was 30 pc at 4 Myr ago and 40 pc 2 Myr ago. These two radii are shown by the concentric, 
dotted circles in Figure [191 The current radius of 50 pc is shown by the solid circle. These circles 
are only sketches of the Shell's extent and ignore differences in the local density of material the Shell 
was propagating into and the effects of a 2 dimensional shape imposed on the projected distribution 
of a 3 dimensional structure. Nevertheless, the isolated Void WTTS are clearly located within the 
central circle and most of them are within a radius of 4° (the CFS radius after only 1-2 Myr) . The 
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Void WTTS have ages comparable to the age of the bubble (jTachihara et al.ll2005l ) 



The ratio of more evolved YSOs to younger YSOs (Class I+F/Class II+III), i.e., the "class 
ratio", in protostellar cores (A'yso < 5) and loose groups (A'yso ^ 5) is 11/9=1.22 and 28/53=0.52 
respectively. The class ratio for the entire region is 35/98=0.36. Taken at face value, this implies 
that the protostellar cores represent a younger YSO population than the loose groups and that the 
loose groups have a relatively narrow range of YSO ages. The class ra tio for distributed YSOs and 
for loose YSO groups averaged across the c2d clouds is 0.24 and 0.89 (jEvans et al.ll2008l ). The c2d 
average for distributed YSOs is similar to the class ratio for the entire Cepheus Flare and further 
supports the idea that relatively distributed star formation is the dominant mode of star formation 
in this region. 

Our SED analysis has shown that the disk population in Cepheus is skewed towards young, 
accreting disks. The age of the disks would be expected if the YSOs in the groups formed simul- 
taneously due to a fairly recent triggering event, such as the passage of a shock wave. The current 
best estimate for the lifetime of low mass star formation prior to the Class III phase is 2-3 Myr (see 
Evans et al.ll2008l . and discussion therein). The four loose YSO groups are in the zone between the 
current position of the CFS and its estimated position 4 Myr ago. Give n that L1251 and L1228 
are believed to be on opposite sides of the Cepheus Void (IKun et al.ll2008l ). projection effects could 
mean that L1251 is closer to the Shell wall than it actually appears. Therefore, the general empir- 
ical ages of Class II dominated groups are consistent with their formation at a time when the CFS 
was coincident with their position. 



The Void WTTS (jTachihara et al.ll2005l ) occur at the same galactic latitude (6 ^ 16°) as L1241, 
but they are distributed towards the projected interior of the Cepheus Void at 115 < / < 123°. 
The region on the other side of L1251 and L1228, in the longitude range 107 < / < 110°, appears 
in Figure [19] to hold a reservoir of material that has not yet started forming stars in the same 
manner. L1241 also represents a large mass of material that has not yet started forming stars. 
It has usually been placed at the same distance as L1251 (see Section [2]), but the lack of star 
formation may indicate that it is actually at a different distance and, like the material in the range 
107 < I < 110°, has yet to encounter the CFS. Alternativ ely, the CO emission to wards L1241 may 
be the superposition of two or more low density clouds, lYonekura et al.l (jl997l ) list two different 
components to Vlsr towards L1241. 

Simulations of the impact of ionizing radiation on a turbul ent ISM show tha t it can enhance pre- 



exist ing density contrasts and the efficiency of star formation (jPale et al 



2007 



Gritschneder et al 



20091 ). This enhancement could explain wh y we see more efficient s tar formation in the cores 
towards the CFS. The surface density of the iGritschneder et al.l ([20091) simulation at t = 250 kyr 
(their Figure 1) resembles the distribution of CO seen in the Cepheus Flare (our Figure [T9]) . 
However, the time scale and physical dimensions of their simulation are shorter and smaller than 
inferred for the CFS (see above). 



Moving further west there is a zone clear of CO emission at / ~ 106° and then we reach 
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the locations of L1172+L1174 at / ~ 104° and L1148+L1152+L1155 at / ~ 102°. Even here, 
furthest from the CFS, it appears that some degree of triggering is at work. The NGC 7023 
cluster in L1172 is being shaped by the powerful illumination of HD 200775 and holds the largest 
concentration of YSOs in the region we have surveyed. We presume that the formation of these 
YSOs has been triggered by compression of the material around the reflection nebula caused by HD 
200775. L1148+L1152+L1155 are comprised of a series of pre-stellar cores and protostellar cores 
that contain only a few YSOs each. It would appear that this final r egion is the rnost u nperturbed, 
but again there is some evidence of outside influences on the cores (jNutter et al.ll2009l ). 



The general scenario for star formation in the Cepheus Flare has the initial conditions of a 
turbulent ISM where distributed, quiescent, low-efficiency star formation is the norm. Approxi- 
mately 7 Myr ago, a progenitor high-m ass star with a st rong circumstellar wind formed at the 
center of what would become the CFS flOlano et al.ll2006l ). The wind or radiation from the star 
compressed the surrounding gas into pillars/filaments in which a series of new YSOs formed. The 
influence of the progenitor star eventually stripped the natal r nolecular gas from thi s generation 
of YSOs and left them behind as the Cepheus Flare WTTSs (|Tachihara et al.l hoQfh . The wind 



from the progenitor continued building up an expan ding shell of materia l (|01ano et al.ll2006l ) until 
it went supernova approximately 40 kyrs years ago (jGrenier et al.lll989l ). A second generation of 
star formation occurred on surfaces facing the interior of the wind bubble. The passage of the 
CFS enhanced a core's SFE relative to the original quiescent mode of star formation and produced 
a series of small YSO groups. The material of the clouds that surrounded the YSO groups was 
stripped away before the co-evolving group, on average, reached the Class III phase. 



8. Summary 

• We have presented SpitzerlRAC and MIPS observations of the L1148-t-L1152-hL1155, L1172-FL1174, 
L1221, L1228, L1241, and L1247-I-L1251 dark cloud associations in the Cepheus Flare region 

{D ~ 300 pc) as part of the Spitzer Gould Belt Legacy Survey (SGBS) of star formation 
within 500 pc. 

• The SGBS delivery catalog (which is available from the SSC as a legacy data product) contains 
over 71 thousand sources within the area common to all wavelengths and detected in at least 
one IRAC band. Of these, 6.5 thousand sources have detections at all four IRAC wavelengths 
and 888 also have detections at MIPS 24 fim. Across the entire area 392 sources were detected 
at MIPS 70 fim. 

• Three color-magnitude schemes based on 2MASS-MRAC-FMIPS, 2MASS-hMIPS, and IRAC- 
only photometry were used to reject background galaxies from the catalog. A total of 133 
candidate YSOs were identified in this manner. Two-thirds of these were classified by the 
infrared spectral index as Class II YSOs (generally classical T Tauri YSOs) and one-fifth 
were classified as Class I YSOs (embedded protostars). Cross- identifications were made for 
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93 candidate YSOs in the GSC-II and 20 candidate YSOs in the IRAS Faint Source Catalog. 
Ten additional previously known YSOs that were not identified by Spitzer were found to be 
coincident with entries in the catalog. 

The Cepheus Flare Spitzer luminosity function peaks at log(L/LQ) ~ —1.5. Above this 
value, the luminosity function has a power law index of 1.6 in agreement with that found for 
/i?j45'-only YSO candidates. 



SED modeling was conducted, following iHarvey et al.l (j2008l ) and lCieza et al.l (120071 ). to esti- 



mate the degree of infrared excess for Class II and Class III YSOs. The majority of the YSOs 
was found to have accretion-style disks. The values of Ay estimated from the 2MASS J-K^ 
color were over-estimated for known variable stars. We suggest that YSO #83 is a possible 
transitional disk candidate. 

Comparison of 2.5' resolution extinction maps to 160 ^m emission maps showed 14-18 dense 
cores (depending on consolidation) split between the six dark cloud associations. Three of 
these, L1155, L1241, and L1247, are confirmed as starless above the Ay = 5 contour. It 
was found that the morphology of the 160 /im and higher resolution Spitzer Ay maps agreed 
particularly well for quiescent pre-stellar cores, but diverged in less quiescent regions (e.g., 
the heart of the NGC 7023 reflection nebula). 

Five YSO associations with peak stellar mass volume densities greater than a theoretical 
cluster disruption density were found. L1228N, L1228S, L1251B, and L1251A are small 
groups with 5-15 members each. The larger NGC 7023 cluster contains 32 YSOs and 21 non- 
YSO members. The star formation rate for dense cores with and without formal associations 
of YSOs was found to be ~ 8% and ~ 1% respectively. No difference in the mass of pre-stellar 
and protostellar cores was found. 



We are grateful to Maria Kun and her co-authors for a preprint of their chapter, "Star Forming 
Regions in Cepheus" , from the Handbook of Star Forming Regions, Vol 1: The Northern Sky, (2008, 
ed. B. Reipurth, ASP Monograph Series). 

JMK thanks UK STFC for PDRA funding through the Cardiff Astronomy Rolling Grant. 

This work has made use of data and resources from the Spitzer Space Telescope, the Two Micron 
All Sky Survey, the Digitized Sky Survey II, the Guide Star Catalog II, the James Clerk Maxwell 
Telescope, the Canadian Astronomy Data Centre SCUBA Archive, and the SIMBAD database. 
The Spitzer Space Telescope is operated by the Jet Propulsion Laboratory, California Institute of 
Technology under a contract with NASA. The Two Micron All Sky Survey is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space Administration and the National 
Science Foundation. The Digitized Sky Survey-II is based on photographic data obtained using 



- 36 - 



The UK Schmidt Telescope. The UK Schmidt Telescope was operated by the Royal Observatory 
Edinburgh, with funding from the UK Science and Engineering Research Council, until 1988 June, 
and thereafter by the Anglo- Australian Observatory. Original plate material is copyright (c) of the 
Royal Observatory Edinburgh and the Anglo- Australian Observatory. The plates were processed 
into the present compressed digital form with their permission. The Digitized Sky Survey was 
produced at the Space Telescope Science Institute under US Government grant NAG W-2166. The 
Guide Star Catalog-II is a joint project of the Space Telescope Science Institute and the Osser- 
vatorio Astronomico di Torino. Space Telescope Science Institute is operated by the Association 
of Universities for Research in Astronomy, for the National Aeronautics and Space Administra- 
tion under contract NAS5-26555. The participation of the Osservatorio Astronomico di Torino is 
supported by the Italian Council for Research in Astronomy. Additional support is provided by Eu- 
ropean Southern Observatory, Space Telescope European Coordinating Facility, the International 
GEMINI project and the European Space Agency Astrophysics Division. The James Clerk Maxwell 
Telescope is operated by The Joint Astronomy Centre on behalf of the Science and Technology Fa- 
cilities Council of the United Kingdom, the Netherlands Organisation for Scientific Research, and 
the National Research Council of Canada. SCUBA was built at the Royal Observatory, Edinburgh. 
The Canadian Astronomy Data Centre is operated by the National Research Council of Canada 
with the support of the Canadian Space Agency. The SIMBAD database is operated at CDS, 
Strasbourg, France. Marvelous fellows, all of them. 

Facilities: Spitzer, CTI0:2MASS (), JCMT 



A. Discussion of Individual Regions 



The SCUBA Legacy Catalog (IDi Francesco et al.l l2008l ) was used to identify regions of the 



Cepheus Flare that have been observed with SCUBA. The majority of the cores in the Cepheus 
Flare have been observed with the SCUBA scan-map method at one time or another and with 
varying levels of signal-to-noise. The exceptions are L1251W and L1228S which were not observed, 
and PV Cep which has only been observed with the jiggle-map method. The SCUBA Legacy 
maps have been filt ered to remove large sca le angular variations which may or may not be real 



for any given map (|Di Francesco et al.ll2008l ). The Cepheus Flare scan-maps were rereduced by 
us to compare possible extended structure against the MIPS 160 /im emission. The entire set of 
raw data for the Cepheus Flare scan maps and the jiggle maps for PV Cep and LI 172 (which is 
higher signal-to-noise than the L1172 scan-map) were downloaded from the archive and reduced as 
described in Section 15.11 A median baseline was used during the reduction of the scan-maps. The 
data have been smoothed to a FWHM resolution of 18" to improve the image quality. 

Figure [20l shows IRAC 3.6 fim images towards twelve Cepheus Flare cores. The IRAC images 
are shown with a log stretch to emphasize faint nebulosity. The white contours show the SCUBA 
850 fim data. The contours are at intervals of 2 cr (3 cr for L1228N and 10 a for PV Cep) and 
start at 3 a. The local 1 a rms was measured as the point-to-point variation towards off-source 
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regions and was typically 10-30 mJy per 18" beam. The dashed contours show Ay from Figure [T71 
at levels of 5, 9, 13, 17, and 21 mag. The positions of YSO candidates are shown by star markers 
and each of them is labeled with their index number. Additional objects of interest are shown by 
white crosses and are discussed in the following sections. Selected outflows are shown by dashed 
lines. Table [THl lists the driving sources, Herbig-Haro numbers (where relevant), and references for 
the outflows shown in Figure [20l The relation of the cores to each other is shown in Figure [H 



A.l. The L1148+L1152+L1155+L1157 Dark Cloud Association 

The Ay and 160 ^m maps both show a ring- or loop- like structure (see Figures [T71 and [TH]) that 
includes the dense cores L1148, L1155, L1152 and L1157. This ring is remarkable because three 
out of these four cores contain relatively isolated Class protostars suggesting that star formation 
is still ongoing within this structure. Two bright point sources are visible at 160 /im, the Class 
protostar L1157 and the Herbig AeBe star PV Cep. 



A. 1.1. PV Cep 



PV Cep is a bright Herbig Ae/Be star (iLi et al.lll994l ). Figure [20k shows the region around PV 
Cep. The SCUBA 850 //m emission is point-like, but there is a slight vein of low-level emission that 
poir its to the south. T his could be CO contamination from the outflow in the 850 fim waveband (e.g, 



see iChini et al.l 120011 ). The MIPS 160 fim emission is also point-like, but is saturated in the BCD 
images. The IRAC image shown in Figure [20k is a BCD image from the SSC archive and still has 
some of the artifacts that were removed from other images by the c2d/SGBS processing. Inspection 
of the SSC images shows that PV Cep is also saturated at all four IRAC bands. Nevertheless, the 
breadth of additional photometry does allow for a fairly complete reconstruction of its SED as 
shown in Figure [TOl The SED has had to be scaled downwards to fit on the same axes as the other 
SEDs as it is so bright. 

The white crosses in Figure [20k show the positions of the HH 315 outflow. The closest two HH 
objects to PV Cep have been connected to show the orientation of the outflow. Lines connecting 
paired HH objects on eit her side of PV Cep do not cross in the same location and were taken by 
Goodman &: Arcd (|2004l ) as evidence that the source was moving to the west. Taken with other 
evidence, they concluded that PV Cep had an unusually high local motion and has escaped from 
the NGC 7023 cluster. 



Kauffmann et al 



A. 1.2. L1148 

(j2005l ) studied L1148 and reported the discovery of a very low luminosity 
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object they called L1148-IRS. To the south of LI 148 in Figure [T71 the Ay = 5 contour breaks in 
two between a northern component containing L1148-IRS and a southern component that runs off 
the bottom of the high resolution Ay box. This southern core is L1147. 

Figure [20b shows the region around L1148 at IRAC 3.6 fxm. Care must be taken in interpreting 
the SCUBA 850 /im data as much of the brighter emission is towards the edge of the map and could 
easily be Fourier art ifacts introduced by the SCU BA scan-map restoration. Comparison with the 



MAMBO filaments (Kauffmann et al 



2005 



20081 ). the Ay contours, the MIPS 160 /xm emission 



and Akari far-IR emission ([Nutter et al.ll2009l ). however, show common features that supports the 
probability that SCUBA is tracing extended structure. Taken together, they appear to show a 
loop or square of i naterial /emission with L 1148-IRS on the northeastern side. The two MAMBO 
dust filaments that iKauffmann et al.l (120051 ) observed towards LI 148 are also part of this loop. The 
filament that contains L1148-IRS (YSO ^5 in Figure [20b) follows the line of the extinction, but the 
brighter MAMBO filament is 3' to the east and is coincident with a 160 /xm feature. The SCUBA 
contours also show these two filaments aligned NE-SW with L1148 IRS at the very top of one of 
them. 

The extinction peaks around L1148-IRS and to the north-west. The 160 fiui emission also 
peaks to the north-west, but is relatively weak around L1148-IRS and is stronger towards the 
south. iNutter et ahl (|2009l ) observed the L1148-t-L1155 filament at 90, 140, and 160 /um with the 
Akari far-infrared satellite. They showed that the extended far-infrared emission towards L1148 
and L1155 followed particular edges of regions of high visual extinction and attributed the difference 
in the two distributions in LI 155 to the effect of an exterior source heating the dust along one side 
of the core. They also detect the source L1148-IRS strongly at 90 //m. There is good agreement 
between the 160 fim emission as observed by MIPS and Akari. 



A. 1.3. 11152 



Chapman et al.l (j2009l ) observed L1152 with Spitzer as part of their study of the extinction law 
in four dense cores. They found three YSO candidates in this field - the two Class II sources YSO 
#2 and #3 and the Class source YSO #1. The IRAS sources F20353+6742 and F20358+6746 
are associated with YSOs #1 and #3. Figure [20b shows the region aroun d L1152 at IRAC 3. 6 ^m. 
An alternative reduction of this SCUBA data was originally presented in lYoung et al.l ( 2006 ). 



The Herbig Haro (HH) 376 jet cuts throug h this region - its orien tation is shown by the dashed 
line and the position of HH 376A is indicated (jReipurth et al.l 119971 ). It is immediately noticeable 
that the HH 376 outflow is aligned tightly with the objects of interest in this field. The northern 
end of the outflow is coincident with the c qmetary nebula GM 3-1 2 (RNO 124 = YSO #3), this 
is a conical nebula with two helical arms (jMovsessian et al.ll2004l ). The HH-obj ect 376A has a 
bow-shock structure su ggest ing that it is moving to the S-W away from YSO #3 (jReipurth et al 



1997; Movsessian et al. 



20041 ) . Given the proximity of the isolated group of 850 fim contours to HH 
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376A, it is possible that we are seeing the bow-shaped HH376A (jMovsessian et al.ll2004l ) plowing 
into the clump of material that is ahead of it. 

YSO 7^1 is associated with the peak of the SCUBA emission in the map. It has a long bipo- 



lar in olecular (CO) outflow that is offset from the direction of the HH 376 outflow (iBontemps et al 



19961 . shown as the shorter dotted-line i n Figure [20b). The HH 376 outflow could be powered by 
YSO #1 (e.g., see iReipurth et al.l 119971 ). but given the offset between it and YSO this would ap- 
pear unlikely. The projected line of the HH outflow, however, is close enough to the dense material 
seen at 850 ^m that the outflow could possibly have influenced the evolution of YSO #1 or even 
triggered its formation. There is a region of 3.6 fj,m nebulosity that stretches away from YSO #1 
parallel to the directio n of the HH 376 outflow. YSO #1 does s ho w faint nebulosit y in JH K which 



is clipped to the NE (jHeyer et al 



1990 



Connellev et al.l 120071 ). IChapman et al.l ([20091) identify 



the nebulosity as an outflow that is visible from 2MASS Kg to IRAC 8.0 /im that is altering the 
dust grain properties within surrounding material. This nebulosity is coincident with the 850 /zm 
contours. 



ClarkI (|l99ll ) suggested that YSO #1 is a very young candidate protostar. The spectral index 



and Tboi estimates of the classification of YSO #3 and #1 agree that the former is more evolved 
than the latter. YSO #1 has a value of Tboi = 33 K suggesting that it is a Class source. It also 
has -Z^submm/^boi = 5.0% which is the canonical value dividing Class and Class I sources. YSO 
#3 has neither a conti nuum detection at m illimeter wavelenghts (e.g., a disk; iTerebey et al.lll993l ) 
nor a maser detection (jFuruya et al.ll2003l ). 



A. 1.4. L1155 



The 160 /im emission in LI 155 is fragmented within the Ay = 1 contour. There is a clear 
east-west break in the 160 (im. emission between the smaller eastern core and the western core. 
The larger western core has t wo peaks which correspond to the position of L1155C (north-peak) 
and L1155H (sout hern peak) (IKirk et al.ll2005l ). The smaller eastern core has been listed as either 



L1155E or LI 158 (Kun et al 



SE of L1158 (jKirkeialJl2005|) 



20081 ). The L1155D core is the faint peak of emission slightly to the 



There is general agreement between th e distribution of 160 ^m and 850 /^m emission towards 
LI 155, as seen with other pre-stellar cores (jKirk et al.l 120071 ) . although the SCUBA emission does 
peak slightly further to the north. The strengths of the L1155C and L1155H 160 /im peaks are 
approximately equal . Both the 160 /xm and 8 5 /zm emission run along the side of the low-resolution 
Dobashi et al.l (120051 ) Ay map. iNutter et al.l (120091 ) used a comparison of Akari, SCUBA, and ISO 
data to show that there was a monotonic spatial shift of the emission peak with wavelength that 
was caused by an edge-to-center negative temperature gradient of 2 K created by the external 
illumination o f L115 5 by the nearby A star BD+671263. This star was one of those used by 
Straizys et al.l (119921 ) to derive his distance estimate to these cores (see Section [2]). 
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A.1.5. 11157 



The L1157-MM embedded protostar is usually given the sa me name as the dark cloud. It was 
not observed by our campaign, but was observed with Spitzer bv iLoonev et al.l (|2007l ) who detected 
a flattened pseudo-disk in absorption at 8 /im that matched a similar structure see in N2H^ and 
DCO^. The orientation of the disk was perpendicular to the orientation of the YSO's outflow. It 
was added to our catalog as YSO #134. It has a bolometric te mperature of 37 K, a -Lsubmm/^boi 
ratio of 2.8%, and is listed as a Class source in the literature (jKun et al.l 120081 ). 



Figure [20b shows the region around L1157 at 3.6 fim and SCUBA 850 /im. Both tracers show 
the same basic pattern, a cross-like structure formed by the N-S outflow and EW nebulosity. The 
outflow from L11 57 has been exte ns ively studied and i s considered the "prototype for chemically 
active outflows" (jKun et al.l 120081 ). IChini et al.l (j200lh compared 850 /im emission to CO (1-0) 
emission towards LI 157. They saw a similar alignment of the outflow to the 850 /xm emission and 
concluded that it must contain significant CO line contamination. 

The 160 /um emission across L1157 is dominated by a bright point-source centered upon the 
YSO. Fainter extended emission appears to the SE in the approximate position of the contours in 
the bottom-left corner of Figure [20j There is no such emission to the north of the YSO seen at 
either 160 /xm or 850 /Lim and the YSO sits on the northern edge of the core as mapped by the 
Dobashi Ay extinction. This suggests that the northern portion of the outflow is emerging into a 
lower density environment than the southern outflow. 



A. 2. The L1172-KL1174 (NGC 7023) Dark Cloud Association 

The L1172-I-L1174 (NGC 7023) complex is a clustered region of star formation. The extinction 
map in Figure [T71 shows the same T-shaped cloud as Figure[3]with regions of high extinction towards 
the extremities of the T and a void near the center. The NGC 7023 reflection nebula, otherwise 
known as the Iris Nebula, sits at the center of the T and is powered by the seventh magnitude Herbig 
AeBe star HD 200775. The 160 /im emission shows the same basic T-shape, but is dominated by 
bright emission from the center of the nebula that completely saturates the MIPS detectors. The 
dark nebula associated with the dense material east and west of NGC 7023 is called L1174. We use 
the designations L1174A, L1174B, and L1174C to refer to the material to the east, to the west, and 
in the center of NGC 7023 respectively, based on decreasing order of peak Ay. The dark nebula 
that forms the T-stem is LI 172. 

Figure [20f shows a close up of the region around the NGC 7023 nebula and cluster. The 
asymmetric bi-conical cavity to either side of HD 200775 is clearly visible. The position of HD 
200775 is shown by the central source YSO #138. There is excellent agreement between the 
lower SCUBA contour and the Ay contours. The SCUBA and extinction contours show the dense 
material that forms L1174, crossbar of the T, with the reflection nebula immediately to the south. 
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The eastern lobe of the nebula appears in bright 3.6 /im emission while only the partial rim on 
the western lobe appears. The eastern lobe is constrained by a greater level of extinction than the 
western lobe. The 3.6 /im emission traces the southern and western part of the western lobe, but 
the northern part is much fainter. This same pattern can be seen in the 160 //m map where the 
south-west part of the lobe appears in strong emission. Where the broken-western lobe encounters 
the edge of L1174B at 2l''00'"20^ -F68°13'00", there is a smah group of YSO candidates. This is the 
NGC 7023 Tight B group as identified in the clustering analysis. An enlargement of these YSOs 
is shown in the smaller image in Figure [20^ where they appear to form a small ring-like nebula 
visible at 3.6 fim. 

L1174C, the central extinction peak just to the north of HD 200775, appears as a short filament 
with three separate 850 fiui peaks. The middle peak is coincident with the embedded YSO #34. 
It has insufficient data points to calculate Lsubmm) but it has a bolometric temperature of just 
14 K making it a probable Class source. The southern edge of the filament follows the line of a 
3.6 fim filament. This filament has been mapped with multiple wavelengths, includin g earlier Spitzer 
observations, and traces th e dissociation front at the end of the dense molecular gas (lAn &: Sellgren 



2003; 



Werner et al.l l2004al ). Conditions along the NGC 7023 diss ociation front ha ve made it an 



ideal location to test dust-a ssociated photo-lum i nescence mode l s (IWitt et al.ll2006l ) and emission 
properties of PAH features (jFlagev et al.ll2006l ). I Werner et al.l (|2004al ) reported the existence of 
a small ring of IRAC 4.5 //m emission betw een HD 200775 and the bright ridge just north of it, 
also seen at 3.3 /im bv lAn &: SellgrenI (j2003l ). The feature is also present in our maps, but is not 
highlighted by the stretch we have used to display the images. The top of the bright emission from 



the n ebula forms a sharp dark line that is coincident with a strong linear HI filament (jFuente et al 



19981 ). The strongest part of the filament is also coincident with the 850 /im emission from L1174C. 

Figure [20f and the extinction map show that the eastern lobe is truncated at 2l'*02™'10* by 
the edge of the extinction from L1174A. It is a large block of extinction that runs eastwards for 
about 10'. It is noticeable that star formation is not occurring within this extinction and that all 
YSO candidate markers are positioned towards its edge or at lower Ay. 

The molecular gas in NGC 7023 has been extens iv ely stud i ed, bu t the associated cluster 
has received relatively little attention (|Kun et al.ll2008l ). ISellgrenl (jl983l ) identified thirty cluster 
members of various spectral types. We identify 9 of Sellgren's cluster members as YSO candidates 
(these are labeled as such in Table [6]). Non-YSO cluster members are shown on Figure [20f as white 
crosses and in Figure [T7] as yellow crosses. Sellgren identified seven probable pre-main sequence 
stars based on their variability, infrared excess, and hydrogen-line emission. Of those seven, we 
identify all but the variable stars SX Cep and HZ Cep as YSO candidates. 

The majority of the YSO candidates in this region are part of the Tight A Group and the 
majority of these are situated northwards of the two lobes of the nebula in the region coincident 
with the area of the highest SCUBA emission. The smaller eastern lobe is coincident with several 
YSOs and stellar cluster members, but the western lobe is noticeably less populated. The majority 
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of the stellar cluster members which do not have YSO candidate cross-identifications are situated 
to the south and east of HD 200775. The clustering analysis was repeated to include the non-YSO 
cluster members. The revised 25Mq pc~^ contour (see equation [5]) is shown in Figure [T7] by the 
red contour. The inclusion of the non-YSO members extends the Tight A group to cover most of 
the eastern lobe. 

Figure [20h shows the L1172D region which is at the northern tip of the southern extinction 
maximum in Figure[T7l I t contains a t r iplet o f YSOs in its densest region. The region was originally 



mapped with SCUBA bv lVisser et al.l (|2002l ) who labeled the peaks associated with YSO #49 and 
#50 as L1172-SMM1 and L1172-SMM2. They believed L1172-SMM2 to be starless, but the high 
sensitivity of our IRAC observations has revealed the presence of a faint YSO at that location. 
The third submillimeter source, L1172-SMM3, is 3' further to the south and is confirmed starless. 
SMMl and SMM2 both have very low values of Tboi, 42 K and 24 K respecti vely, indicating tha t 



they are Class sources. SMMl has a broad CO outflow that was detected by lMyers et all ()l988l ) 

and its i>submm/-^^boi ratio is greater than 5%. 

The Class II YSO #120 (not shown) is coincident with the listed position of L1171 and a small 
clump of Ay = 1 extinction. 



A.3. L1221 



The protostellar core L1221 lies towards the southern edge of the G109+11 infrared loop (the 
top of which is shown in Figure [1]), closer to the Galactic plane than the other regions in our 
sample. Figure [20l shows two SCU BA sources towards L1221. The northern submillimeter source 
contains a pair of YSO candidates. I Young et al.l (|2009l ) first published the Spitzer data for LI 221 
and detected three infrared sources embedded within the SCUBA emission which they labeled IRS 
1, 2, and 3. These correspond to the YSO candidates we index as #63, #105, and #64 respectively. 
YSOs #63 and #64 have values of Tboi of 66 K and 21 K and -^^submm/-^^boi ratios of 4 and 9% 
respectively. Based on these values YSO #64 is a definite Class source while YSO #6 3 is a 
Class I source. The approximate orientation of a broad molecular outflow (|Lee &: Hdl2005l ) that 
has been detected towards L1221 IRS 1 is shown by a dashed line in Figure [20l. The properties of 
the SCUBA envelopes around IRSl and 3 are very similar and the relative position of the infrared 
sources to the env elopes are virtually identical, but the protostars themselves are not as similar 
(|Young et ahlbood ). 



A. 4. The L1228 Dark Cloud Association 

The extinction map towards L1228 shows two distinct peaks. A northern peak at a declination 
of +77°35' and a lower Ay southern peak at +77°08'. We caU these two cores L1228N a nd L1228S 



respectively, although we note that L1228N has also been called L1228A in the literature (|Kun et al 
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20081 ). L1228N and L1228S each contain small groups of YSO candidates with sufficient stellar 



densities to be picked up by the clustering analysis. The SIMBAD database shows 10 YSOs in the 
entire L1228 region of which 4 are just outside of the region mapped with Spitzer (all scattered to 
the west). L1228N has received more attention than L1228S and is the site of several Herbig Haro 
outflows. 



A.J.I. L1228N 



The extinction map towards L1228N shows a core that is elongated approximately E-W while 
the distribution of YSO candidates is approximately N-S. The 160 ^m. map shows a similar mor- 
phology, but is dominated by a central bright point source. The source is coincident with YSO 
#9 and the position of IRAS 20582+7724. Figure [20) shows IRAC and SCUBA emission towards 
the L1228N region. It also shows two families of HH objects as white crosses and the approximate 
orientation of the HH 199 outflow which em anates from YSO #9 and the HH 200 outflow which 
emanate s from YSO #7 (IBallv et al.lll995l ). Additionally, YSO #9 has an east-west molecular 
outflow (lArce &: Sargentll2006l ) which appears as faint nebu losity in the IRAC ima ge. An enhanced 
image generated from the same IRAC data was shown by IVelusamy et al.l (120081 ) as part of their 
demonstration of the application of the HiRes image deconvolution technique to Spitzer data. 



Chapman et al.l (|2009l ) examined the mid-infrared extinction in L1228N and concluded that 



the outflows were altering the extinction law by destroying larger dust gr ains within the outflow 
cavities. The four VLA sources were detected in this region at 3.6 cm by iReipurth et al.l (|200J). 
VLA 1 and VLA 4 match the positions of YSO candidates ^,^9 and ^^7 from our survey, but no 
YSO matches are found for VLA 2 and VLA 3 though they may have IRAC 3.6 /nm detections. 
The positions of VLA 2 and VLA 3 are labeled in Figure [20j . 



The SCUBA emission towards L1228 is dominated by a strong point source coincident with 
the position of YSO #9 and a similar point source seen at 160 //m. Surrounding this is low-level 
emission that is elongated in a N-S direction, perpendicular to the direction of the HH 199 outflow. 



The east-west 3.6 ^um nebulosity surrounding YSO #9 i s co-incid e nt wi th a 2.2 /xm jet (jBally et al 



19951 ) . The same emission was seen in the K band by iHodappI (119941 ) who classified the sources 
surrounding it as a small cluster - albeit at the limits of their definition of a cluster. YSO #9 
has a broad SED, meaning that it is spectrally classified as a flat spectrum source rather than 
a borderline O/I source as would expected by its bolometric temperature (79 K). YSO #7 is not 
detected by SCUBA, but has a bolometric temperature (54K) below the Class cut-off. 



A.4.2. L1228S 



Figure [20k shows L1228S. There are no SCUBA data for this core. The 5|?iteer observations of 
L1228S were first presented bv [Padgett et al.l (|2004l ) who showed a cluster of 9 sources with Class 
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II or III SEDs. L1228S is situated at the furthest southern extent of the Dobashi Ay = 1 contour. 
As noted previously, the L1228S extinction appears to be separated from the 160 /xm peak. The 
proximity of L1228S to the edge of the cloud and the proximity of LI 228 to the Cepheus Flare 
Shell leads to the possibility that the cloud around t hese YSQs is bein g removed by the passage of 
the Shell, leaving the YSOs free of their natal cloud (|Kun et al.ll2008l ). 



A.5. L1241 

L1241 is a relatively massive starless core between the active regions of L1251 and L1228. 
L1241 is not shown in Figure [20] as it has not been observed with SCUBA and is not associated 
with any sources in our YSO catalog. It is associated with a large mass of CO on the west side of 
the Cepheus Flare (see Figure [T] and Table [1]) , but it appears to show little evidence for current star 
formation. Only a single YSO candidate, a Class III YSO at 21^56'^13.3^+76°58'14.2", was found 
towards this core, but it is well away from the Ay peak and this core should still be considered 
starless. No YSO candidates were founded in the SIMBAD database. The nearest SIMBAD YSO 



is 2° away ([K98c EM* 61]) and the nearest iTachihara et al.l (120051 ) x-ray selected T Tauri is 3° 



away. Why such a large mass of gas is not undergoing star formation is not immediately obvious as 
the regions L1251 and L1228 seem to be under going star formation triggered by the passage of the 
Cepheus Flare Shell. If L1241 is not undergoing star formation then it may be because it has not 
yet encountered the Cepheus Flare Shell and is at a marginally different distance than originally 
thought. 



A.6. L1247+L1251 



The map of Ay towards the L1247-I-L1251 dark cloud association shows a chain of cores 
proceeding east-west with the peak Ay decreasing towards the east. The cores L1251A, B, and 
W all contain YSO candidates (see Figure [20l-n). but L1247 appears starless (not shown). The 
YSO candidates in this region form th ree distinct grou ps - the L1251A and L1251B groups and 
a small quart e t in L 1251W. We follow iLee et al.l (j2006l ) in naming these cores A and B after the 
Sato &: Fukuil (119891 ) outflow sources and we adopt the name L1251W to refer to the elongated 
western core. 

The L1247 core has no associated YSOs and is confirmed starless within the sensitivity of our 
survey. MIPS 160 /xm was not taken for the c2d cores region so there is no 160 //m data available 
for L1251. Data were taken for L1247 however, it shows excellent correlation between the 160 i-im 
emission and the visual extinction. The 160 /xm emission and the extinction map peak in the same 
position and show a narrow extension/filament to the north of the peak. 

In L1247-I-L1251 a total of 39 YSO candidates was identified from the Spitzer photometry. Of 
these, 16 are not in the SIMBAD database. The majority of these, however, are from the L1251A 
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core and will be discussed in the paper bv lLee et al.l ((20091). There is a one-to-one correspondence 
between the YSQ candid ates #66, 67, and 68 in L1251A and the ammonia cores T3, T2, and Tl 
(jToth fc WalmslevI 119961 ) ■ The outflow HH 189 (jEiroa et al.l 11994 ) is shown by the white crosses 
and dashed line in Figure [20hi. It emanates from the tight YSO group associated with the small 
knot of SCUBA emission, but it is not entirely clear which YSO is driving it. The approximate 
orientation of the HH 149 is shown by the dashed line in Figure [20t i. The driving source outflow is 



coincident with YSO #143 (ISato &: Fukui 



198S 



Balazs et al.lll992l ;). 



The L1251B cluster was studied bv lLee et al.l (|2006l ) using the same data presented here. They 
found tha t a tig ht cluster of YSO candidates centered around the bright source L1251B IRS 1. See 
Lee et al.l (j2009l ) for further discussion of this region. 
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Fig. 1. — : A finding chart for the -^300 pc Cepheus Flare region showing the positions of the regions 
discussed in this paper. The areas mapped with Spitzer are shown as black (24 //m MI PS) and whit e 
(3.6 IRAC) footprints. The labels are the names of selected Lynds dark nebulae (lLyndslll962l ). 
Also shown are the reflection nebula NGC 7023 and th e high velocity YSO PV Cep that is reported 
to have been ejected from it (iGoodman &: Arcd 12004 ). The underlying gray-scale map show s IRIS 



(Improved Reprocessed IRAS Survey) 100 ^m emission ( Miville-Deschenes Lagache 20051). The 



white contours show integrated CO emission at levels of 5.5, 12, and 18.5 K kms ^ ( Dame et al 



200ll ). The black contours show visual extinction towards the Cepheus region , generated from the 



Digitized Sky Survey with levels at Av=l, 2 and 3 mag (jDobashi et al.ll2005l ). The dashed-black 
lines denote large scale loops and shells (see text for details). The L1221 region is not shown, but 
see Table [1] for its relative position. 
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Fig. 2. — : IRAC and MIPS composite RGB images 4.5 fim (blue), 8.0 fim (green), and 24 fim (red) 
emission towards the regions in this survey. The nebulosity in LI 172 is the NGC 7023 reflectio n 
nebula. The red source at the center of LI 148 is the source discovered by lKauffmann et al.l (|2005l ). 
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Fig. 3. — : Three composite RGB maps showing how the pattern of emission towards the NGC 7023 
region changes between the optical, infrared, and far-infrared regimes. Left: POSS-II Bj (blue), 
Rf (green), and In (red); Middle: IRAC 3.6 fim (blue), 5.8 fim (green), and 8.0 fim (red); Right: 
MIPS 24 fxm (blue), 70 nm (green), and 160 /xm (red). The dense gas appears in extinction in 
the optical image and emission in the far-infrared MIPS image. Embedded YSOs are hidden by 
the extinction in the optical image, but they appear as bright point-sources in the IRAC infrared 
image. 
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Fig. 4. — : Source counts per square degree per magnitude for six bands from the Cepheus Delivery 
Catalog. The gray line shows all sources whilst the black line shows sources with a S/N>3. The 
turn over in each black line is taken as the limiting magnitude for that band. The 2MASS catalog 
only includes sources with S/N>10 so no gray line plot is shown for the Kg band. The dashed black 



line is the expected background source count towards Cepheus calculated from the IWainscoat et al 
(|l992l ) model of Galactic infrared source counts with a visual extinction of Ay = 1- 
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Fig. 5. — : A histogram of log(Pgai) for the sources in the Cepheus catalog that were detected in each 
IRAC band and MI PS 24 /um. The dashed hne shows the log(Pgai) = —1.47 criterion estabhshed by 
Harvey et al.l (j2007l ) as the divide between YSO candidates (YSOc) and galaxy candidates (GALc). 
The majority of the galaxy candidates appear as a distribution at —1 < log(Pgai) < 0, whereas the 
majority of the YSO candidates appear as a separate peak at log(Pgai) = —5. 
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Fig. 6. — : Color- magnitude diagrams after iHarvey et alJ (|2007l ) showing the separation of YSO 
and galaxy candidates and the dividing lines that were used to construct -Pgai- Left panels: SGBS 
Cepheus data. The markers show point-line/extended YSOs (black filled-circles/crosses) and point- 
like/extended background galaxy candidates (dark gray filled/open-circles) that have been sepa- 
rated via their Pgai value. Stellar sources (light gray filled circles) with the same S/N requirements 
as the YSOs and galaxies are shown for comparison. The open square is a source with only an 8 fim 
flux upper-limit, but otherwise satisfies the 5-band YSO selection criterion. Right panels: Contour 
plots of the number density of sources identified as neither YSOs or stars from the SWIRE catalog. 
Dashed black lines show the "fuzzy" borders used to calculate Pgai while solid lines show the hard 
[24.0] < 10 mag cut-off. 
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Fig. 7. — : Color-magnitude diagrams for IRAC and 2MASS/MIPS detections showing the dividing 
hnes used to separate YSO candidates from other sources. Top row: IRAC-only detections. The 
left panel shows the IRAC-only data, the middle and right panels the 5-Band and SWIRE galaxies 
as per Fig ure El The triangle m arkers show PAH-strong sources as identified by the selection 
rules from iGutermuth et al.l (|2008l ). Lower row: 2MASS/MIPS only detections. Left panel shows 
2MASS/MIPS only detections, middle and right panels are 5-Band detections and SWIRE galaxies 
as per Figure [H The solid lines show the selection cut-offs. Dashed lines show a 1 mag offset from 
the selection cut-offs and open markers show YSO candidates within that zone. Non-YSOs in some 
panels have been shown as a density plot for clarity. 
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Fig. 8. — : Upper: Plot showing the equivalence of ax and the Kg- [24.0] color for the 5-Band 
YSOs. The solid line shows the theoretical relationship between a\ and ii's-[24.0]. The dashed 
line is a linear regression to the data. Lower: Histogram of spectral indices for YSO candidates. 
The black-line shows those YSOs detected by IRAC. The gray-line shows all YSOs including the 
2MASS/MIPS YSOs. The dotted lines show the divisions of the classical YSO spectral classes and 
are extended to show the equivalent Kg- [24.0] colors in the top diagram. The majority of YSOs 
have indices in the Class II regime. Not all of the YSOs included in the solid histogram are plotted 
in the upper graph as not all IRAC YSOs have 24 /xm detections. 
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Fig. 9. — : Color-color diagrams showing the correspondence between the colors of the YSO can- 
didates and various color-regions predicted by models of YSO spectra. The points are color coded 
to the YSOs ajR classification - red are Class I sources, green are flat spectrum sources, blue are 
Class II sources, and purple are Class III sources. Approximate regions that contain the Robitaille 
models are shown on each panel by thin lines (jRobitaille et al.l 120061 ) . Th e box with the th ick 
outline plotted on the [5.8]-[8.0] vs. [3.6]-[4.5] plot is the region identified by I Allen et al.l (|200J) as 
being the approximate domain of Class II sources. The three Robitaille regions on this plot show 
the region where any Stage model can be found (bottom right), where the majority of Stage II 
models are found (small six-sided region), and where the majority of Stage I models can be found 
(top). The regions plotted on the [8.0]-[24] vs [3.6]-[5.8] plot approximate regions occupied by (from 
the left) Stage III, II, and I models. The box enclosing the majority of the sources in the H — Kg 
vs. J — H plot shows the region occupied by reddened stellar photospheres. The region to the right 
of this denotes the area where any evolution stage can occur. 
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Fig. 10. — : SEDs for Class I YSO candidates towards the Cepheus Flare. The open circles with 
error bars show the data points from GSC-II, 2MASS, IRAC, MIPS, IRAS, and SCUBA where 
available. Upper-limits arc also shown where available. The solid curve shows a simple grcybody fit 
to the long wavelength data (see text for details). The YSO Id is shown in the top left hand corner 
of each box. If this is followed by "xO.Ol" it indicates that the SED has been scaled downwards by 
two dex in order to place it on the same grid as the other SEDs. 
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Fig. 11. — : SEDs for flat spectral YSO candidates towards the Cepheus Flare. Details as described 
for Figure [TOl 
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Fig. 12. — Class II candidate SEDs. Details as described for Figure [TOl Additionally, the filled 
circles show the dereddened data. The gray lines show two comparison SEDs that have been 
normalized near the peak of the dered dened SED (usually the 2MASS J band). The solid gray line 



is a NEXTGEN profile for a K7 star fjHauschildt et a 



bars is the median SED for a T Tauri star in Taurus (iHartmann et al.ll2005l ) 



19991) and the dashe d gray line with error 
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Fig. 12. — : (...continued) 




Fig. 12. — : (...continued) 
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Fig. 13. — : Class III candidate SEDs. Details as described for Figure [T2l 
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Fig. 14. — : Top: Bolometric luminosity functions for the 133 Spitzer identified YSOs. Each 
luminosity function is offset by 2.0 dcx from the preceding one. The index of the luminosity 
function is shown for each region. No lower limit is shown for bins with only one source. Lower: 
Histogram of bolometric luminosity (Lboi) for all sources (solid line) and infrared luminosity (^ir) 
for sources detected with IRAC (dashed line). Both plots show a break at our limiting luminosity 
of L ~ 0.06 Lq. 
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Fig. 15. — : Histogram of the ratio of disk to star luminosity (^disk/^star) as derived from the SED 
modehng for Class II and Class III candidate sources. The solid line is for all sources and the 
dashed line excludes sources identified with the 2MASS/MIPS scheme. The vertical dashed lines 
show the expected regions for accreting, passive and debris disks. The histogram shows that the 
majority of the disks are of the accreting type. 
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Fig. 16. — : Plot of Aturnoff) the wavelength where the dereddened SED diverges from the normalized 
stellar photosphere, against Oexcess the spectral index of SED points longwards of Aturnoflf- Each 
marker represents a single Class II or Class III candidate as shown by the key. The wavelength 
bands are labeled with their instrument band. The Class III candidates only show an excess at the 

longest values of Aturnoflf- 
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Fig. 17. — : Distribution of YSO candidates, clustering results, and extended structure towards the 
Cepheus dark clouds associations. The greyscale shows the distribution of visual ext inction derived 



via tw o different datasets. The extinction as derived from the Digitized Sky Survey ([Dobashi et al 



20051 ) is shown across the entire map. The higher-resolution Spitzer extinction maps are superim- 



posed over this in the regions where there is IRAC/MIPS overlap (as shown by the boxes). Two 
dashed contours, Dobashi Av=l and Spitzer Av=5, are shown. The colored markers show the 
location and spectral type of the YSO candidates. The colors are the same as used in Figure [9] 
(red/green/blue/purple for Class I/Flat /II/III). The colored contours, green unless noted, show 
from the clustering analysis in units of 0.125, 0.25, 0.5, 1.0 (blue), 2.0, 4.0, 25.0 (yellow) x M0 
pc~^. 
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Fig. 18. — : Distribution of MIPS 160 /im emission towards the Cepheus dark cloud associations. 
The greyscale shows a log stretch between the local minimum (~ 20 MJy/sr) and 316 MJy/sr (NGC 
7023) or 100 MJy/sr (all other regions). The Dobashi Ay=\ and Spitzer Av=5 contours and the 
overlap boxes from Figure [T7] are shown for reference. No MIPS 160 data exists for L1251 so 



an IS QPHOT 200 /xm map scaled to the 160 /im intensity in L1247 is shown instead (jLemke et al 



19961 ). The line separating the 200 fim and 160 fim maps is shown by a dashed line. 
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Fig. 19. — : The distribution of YSOs and YSO candidates towards the Cepheus Flare. A key 
to the different YSO symbol is shown top-right. The areas mapped with Spitzer at 24 ^um ar e 
shown as black footprints. The labels are th e names of selecte d Lynds dark nebulae (lLyndall962l ). 
The greyscale shows log scaled Ha emissi on (lFinkbeineiil2003l ). The contours shows integrated CO 
emission at 5, 10, 15, 20 and 25 K kms~^ ( Dame et al.ll200ll ). The soli d black line shows the present 
extent of the Cepheus Flare Shell that surrounds the Cepheus Void (joiano et al. 2006 ). The two 
dotted circles show the estimated extent of the Shell 2 Myr and 4 Myr ago. The dashed-black lines 
denote the G109+11 and Cepheus 0B2 Bubble. 
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Fig. 20. — : Images of IRAC 3.6 /xm emission towards the centers of the Cepheus Flare cores. The 
sohd contours show SCUBA 850 /um emission in intervals of 2 a starting at 3 (intervals of 3 o" and 
10 a are used for L1228N and PV Cep respectively). The edges of the mapped SCUBA areas are 
shown by the gray contours. The dashed contours show Spitzer Ay in intervals of 4 mag starting 
at 5 mag. The straight dashed lines show the approximate orientation of selected outflows. The 
positions of YSO candidates are shown by the star markers. These are labeled with the YSO index. 
The positions of other objects discussed in the main text are shown by the white crosses. 
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Fig. 20. — : (...continued) 
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Fig. 20.—: 



(...continued) 
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Fig. 20. — : (...continued) 



Table 1. Dark Cloud Association Properties 



Association I b R.A. Dec. Distance Mdoud AV(^^CO) IRAC/MIPS Area 

[°] [°] (2000) (2000) [pc] [Mo] [kms-^] [Sqr Deg] 



L1148+L1152+L1155 


102.4 


15.5 


20^40"^ 


67° 40' 


325±13 


830 


1.0 


0.31/2.46 


L1172+L1174 


104.0 


14.2 


21^*02™ 


68° 10' 


288±25 


730 


1.6 


0.45/1.93 


L1221 


110.7 


9.73 


22'^28" 


69°07' 


250±50 


40 


1.5 


0.05/0.40 


L1228 


111.7 


20.2 


20''57™ 


77°35' 


200+^°° 


230 


1.6 


0.44/0.77 


L1241 


113.4 


16.9 


22'^00™ 


76° 30' 


3oo;f° 


1500 


2.2 


0.30/1.67 


L1247+L1251 


114.1 


14.8 


22'*30'" 


75° 10' 


30015° 


790 


1.9 


0.56/1.10 



Table 2. AORs in the SGBS Cepheus Catalog 



Dark Cloud PID IRAC AORs MIPS AORs 



Association 




Epoch 1 


Epoch 2 


Date 


Pipeline 


Epoch 1 


Epoch 2 


Date 




Pipeline 


L1148 


30574 


0019963136 


0019963392 


2006-11-26 


S15.0.5 


0019962624 


0019962880 


2007-02- 


■26 


S15.3.0 




30574 


0019995136 


0019995648 


2006-11-26 


S15.0.5 












L1152 


139 


0005156352 




2004-07-03 


S13.2.0 


0009438976 




2004-09- 


■22 


S13.2.0 


L1155C 


139 


0005158656 




2004-07-03 


S13.2.0 


0009408768 




2004-12- 


■02 


S13.2.0 


L1155E 


139 


0005159168 




2004-07-03 


S13.2.0 


0009441792 




2004-09- 


■22 


S13.2.0 


L1172 


30574 


0019965184 


0019965696 


2006-11-26 


S15.0.5 


0019994624 


0019994880 


2007-02- 


■26 


S15.3.0 




30574 


0019996928 


0019997184 


2006-11-26 


S15.0.5 














30574 


0025037312 


0025037568 


2008-03-12 


S17.0.4 












L1221 


139 


0005165312 




2004-08-19 


S13.2.0 


0009428736 




2004-09- 


■24 


S13.2.0 


L1228 


104 


0006577664 




2003-12-06 


S14.0.0 


0006577408 




2003-12- 


■08 


S16.1.0 


L1241 


30574 


0019967232 


0019967488 


2006-11-26 


S15.0.5 


0019964416 


0019964928 


2007-02- 


■23 


S15.3.0 




30574 


0020000256 


0020000512 


2006-11-26 


S15.0.5 












L1247 


30574 


0019969024 


0019969280 


2006-11-26 


S15.0.5 


0019996416 


0019996672 


2007-02- 


■23 


S15.3.0 


L1251A 


139 


0005165824 




2004-04-04 


S13.2.0 


0009425664 




2004-11- 


■30 


S13.2.0 


L1251E 


139 


0005167360 




2004-10-29 


S13.2.0 


0009424384 




2004-09- 


■24 


S13.2.0 


L1251C 


139 


0005166592 




2004-05-18 


S13.2.0 


0009424896 




2004-09- 


•24 


S13.2.0 
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Table 3. Spectral Bands In This Survey 



Origin 


Filter 


Wavelength 


Sopt 


Limiting 


Ref 






[Mm] 


[Jy] 


Magnitude^ 




GSC-II'' 


B 


0.44 


4260. 




1 


GSC-II'^ 


V 


0.55 


3640. 




1 


GSC-II'' 


R 


0.64 


3080. 




1 


GSC-II*^ 


J 


0.79 


2550. 




I 


2MASS 


J 


1.24 


1594. 


16.3 


2 


2MASS 


H 


1.66 


1024. 


15.5 


2 


2MASS 


Ks 


2.16 


666.7 


15.3 


2 


IRAC 


IRACl 


3.6 


280.9 


17.2 


3 


IRAC 


IRAC2 


4.5 


179.7 


16.7 


3 


IRAC 


IRAC3 


5.8 


115.0 


15.5 


3 


IRAC 


IRAC4 


8.0 


64.13 


14.3 


3 


IRAS 


IRASl 


12. 






4 


MIPS 


MIPSl 


24. 


7.14 


10.2 


3 


IRAS 


IRAS2 


25. 






4 


IRAS 


IRAS3 


60. 






4 


MIPS 


MIPS2 


70. 






3 


IRAS 


IRAS4 


100. 






4 


MIPS 


MIPS3 


160. 






3 


SCUBA 


SHORT 


450. 






5 


SCUBA 


LONG 


850. 






5 



''Taken as the turnover in the S/N limited source count distribu- 
tion. 

''The GSC-II data have been converted into the Johnson-Cousins 
photometric system from the natural systems used by its component 
surveys 



References. — (1) Bessell 19791: (2)k 


jkrutskie et al. 2006: (3) 


Evans et al. 2007: f4) Neugebauer et al.l 


1984 (5) Holland et al. 


1999 
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Table 4. Detection (S/N > 3) statistics across the delivered area 



Detected with/in. . . Detections 



. . . at least one IRAC band 71076 

. . . all four IRAC bands 6518 

. . . three IRAC bands 6167 

. . . two IRAC bands 34841 

. . . only one IRAC band 23558 

. . . MIPS 1 2140 

. . . all four IRAC bands and MIPS 1 888 

. . . all four IRAC bands and with 2MASS'' 3839 

. . . all four IRAC bands, but not with 2MASS^ 2679 

. . . 2MASS alone'' 1 

. . . IRAC, but not in 2MASS'' 67080 



Excluding Extended Sources 



. . . all four IRAC bands 6183 
. . . all four IRAC bands and with 2MASS'' 3821 
. . . IRACl+2 and 2MASS'' 3971 



''A 2MASS detection is taken as one with S/N > 10 in both 

H and Kg 
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Table 5. Detection of sources above S/N Thresholds for the entire mapped area at each 

wavelength. 



Wavelength (/xm) 
Detection 3.6 4.5 5.8 8.0 24 70 



S/N > 3 80010 61001 

S/N > 5 63148 44093 

S/N > 7 48854 32061 

S/N > 10 32037 19894 

S/N > 15 14672 9379 



14672 10221 9338 392 

8391 6172 2559 349 

6402 4589 1021 180 

4270 3079 389 2 

2082 1626 1 1 



Table 6. Cepheus YSO Candidates 



Id. Name/Position Literature Name'^ 3.6//m 4.5//m 5.8/nm 8.0/im 24.0/im 70.0/im Flags'' Ref. 

(SSTgbsJ) [mJy] [mjy] [mjy] [mjy] [mjy] [mjy] 



1. 5-Band YSO Candidates 



1 


2035463+6753020 


L1152 1; [K98b] 4 




3.53 


± 


0.19 


6.30 


± 


0.33 


5.45 


± 


0.27 


3.31 


± 


0.16 


273. 


± 


25. 


3670. 


± 341. 


F 


1; 


8 


2 


2036116+6757093 


L1152 2 




41.3 


± 


2.0 


47.9 


± 


2.3 


56.2 


± 


2.7 


74.5 


± 


3.6 


173. 


± 


16. 


347. 


± 33. 


FIM 


1 




3 


2036198+6756316 


L1152 3; [K98b] 
[K98c] EM*6; GM 3-12 


6; 


228. 


± 


12. 


243. 


± 


16. 


215. 


± 


10. 


262. 


± 


13. 


634. 


± 


59. 


1080. 


± 102. 


FIM 


1; 


8; 10 


4 


2039430+6708301 






2.39 


± 


0.12 


2.93 


± 


0.14 


3.46 


± 


0.18 


4.39 


± 


0.21 


9.20 


± 


0.87 






FI 






5 


2040567+6723050 


L1148-IRS; [K98b] 8 




1.32 


± 


0.07 


1.31 


± 


0.06 


0.797 


± 


0.049 


1.16 




0.06 


75.3 


± 


7.0 


331. 


± 33. 


FI 


2; 


8 


6 


2055371+7738196 


L1228 4 




5.96 


± 


0.29 


5.13 


± 


0.25 


4.25 


± 


0.22 


4.21 


± 


0.20 


4.05 


± 


0.39 






FIM 


1 




7 


2057067+7736561 


L1228 5; HH 200 IRS 




6.87 


± 


0.38 


9.72 


± 


0.48 


12.5 


± 


0.6 


16.3 


± 


0.8 


269. 


± 


25. 


2060. 


± 236. 


FIM 


1; 


15 


8 


2057079+7736597 


L1228 6 




3.03 


± 


0.15 


3.35 


± 


0.17 


3.59 


± 


0.18 


8.53 


± 


0.41 


153. 


± 


14. 






Flm 


1 




9 


2057130+7735437 


L1228 7; [K98b] 


12; 


260. 


± 


13. 


406. 


± 


23. 


513. 


± 


24. 


742. 


± 


36. 


1860. 


± 


181. 


9000. 


± 983. 


FIM 


1; 


8; 11 






[RR96] HH 199 1 










































00 
to 


10 


2057155+7734236 


L1228 8 




2.39 


± 


0.12 


2.52 


± 


0.12 


2.62 


± 


0.13 


2.83 


± 


0.14 


4.70 


± 


0.45 






Flm 


1 


11 


2057170+7736586 


L1228 9 




39.3 


± 


1.9 


42.5 


± 


2.1 


39.1 


± 


1.9 


41.1 


± 


1.9 


77.1 


± 


7.1 






FIM 


1 


1 


12 


2058400+7727454 


L1228 10 




20.2 


± 


1.0 


23.5 


± 


1.1 


23.6 


± 


1.1 


24.3 


± 


1.2 


41.2 


± 


3.8 






FIM 


1 




13 


2059154+7733229 






25.1 


± 


1.2 


19.5 


± 


0.9 


17.1 


± 


0.8 


20.0 


± 


0.9 


22.8 


± 


2.1 






FIM 






14 


2059175+7728579 






19.2 


± 


0.9 


17.0 


± 


0.8 


14.3 


± 


0.7 


13.3 


± 


0.6 


13.7 


± 


1.3 






FIM 






15 


2059228+6814438 


V*' FT Cep; [K98b] 
[K98c] EM*26 


14; 


348. 


± 


18. 


368. 


± 


19. 


365. 


± 


18. 


500. 


± 


24. 


919. 


± 


86. 


600. 


± 62. 


FIM 


8; 


15 


16 


2100203+6808270 


NGC 7023 S 41 




41.4 


± 


2.0 


44.0 


± 


2.1 


44.2 


± 


2.1 


47.1 


± 


2.2 


119. 


± 


11. 


713. 


± 143. 


FIM 


7 




17 


2100207+6813169 






3.52 


± 


0.19 


13.1 


± 


0.7 


21.3 


± 


1.0 


25.5 


± 


1.2 


602. 


± 


56. 


3280. 


± 391. 


FI 






18 


2100221+6812585 






9.24 


± 


0.53 


31.5 


± 


1.6 


47.4 


± 


2.2 


59.8 


± 


2.8 


142. 


± 


13. 






FI 






19 


2100289+6813023 






21.9 


± 


1.1 


22.5 


± 


1.1 


21.5 


± 


1.0 


23.6 


± 


1.1 


23.1 


± 


2.2 






FIM 






20 


2100321+6812475 






125. 


± 


6. 


108. 


± 


5. 


96.7 


± 


4.6 


89.1 


± 


4.3 


69.0 


± 


6.4 


293. 


± 72. 


FIM 






21 


2100366+6803446 






19.6 


± 


1.0 


12.5 


± 


0.6 


8.58 


± 


0.41 


5.38 


± 


0.25 


1.80 


± 


0.29 






F 






22 


2100380+7706598 






3.13 


± 


0.15 


2.24 


± 


0.11 


1.60 


± 


0.09 


0.961 


± 


0.052 


3.68 


± 


0.38 






F M 






23 


2100555+6811273 






8.94 


± 


0.43 


7.80 


± 


0.37 


6.97 


± 


0.34 


8.37 


± 


0.40 


9.70 


± 


0.95 






FIM 






24 


2100581+6809383 


NGC 7023 S 36 




17.1 


± 


0.8 


14.6 


± 


0.7 


13.8 


± 


0.7 


17.1 




0.8 


20.9 


± 


2.0 






FIM 


7 




25 


2101037+6813093 






24.0 


± 


1.2 


22.4 


± 


1.1 


19.4 


± 


0.9 


20.6 




1.0 


47.6 


± 


4.5 






FIM 






26 


2101087+6812526 






3.74 


± 


0.18 


3.54 


± 


0.17 


3.42 




0.17 


3.80 


± 


0.24 


5.25 


± 


0.61 






FIM 






27 


2101264+6810385 


NGC 7023 S J 




59.1 


± 


3.2 


90.1 


± 


4.4 


81.5 


± 


4.0 


146. 


± 


10. 


221. 


± 


28. 






FIM 


7; 


15 


28 


2101264+7707028 


L1228 Sl(a) 




4.46 


± 


0.22 


4.92 


± 


0.23 


4.84 


± 


0.24 


5.51 


± 


0.26 


21.5 


± 


2.0 






Flm 


3;15 


29 


2101273+6811384 






16.5 


± 


0.8 


16.5 


± 


0.8 


18.1 


± 


0.9 


21.9 


± 


1.0 


31.4 


± 


3.0 






FIM 







Table 6 — Continued 



Id. Name/Position Literature Name'' 3.6/im 4.5/im 5.8//m S.O/xm 24.0/xm 70.0/im Flags'" Ref. 

(SSTgbsJ) [mjy] [mjy] [mjy] [mjy] [mjy] [mjy] 



30 


21 ni 2Q3-)-77n2373 


7,1228 S6 


167. 


-|- 


9. 


115. 


-|- 


7, 


92.0 


-|- 


5.0 


97.5 


it 


5.0 


175. 


-j- 


16. 


145. it 17. 


FIM 




15 


31 


21 (11 297-1-681 3468 




18.7 


-|- 


0.9 


17.6 


-|- 


0.9 


17.9 


-|- 


0.8 


18.4 


it 


0.9 


45.6 


-|- 


4.2 




FIM 






32 


21 ni 31 (1-1-7701 


1,1228 S7- fK98rl F,M*40 


18.8 


-|- 


0.9 


13.2 


-|- 


0.6 


8.97 


-|- 


0.43 


5.81 


it 


0.28 


53.6 


-|- 


5.0 


65.5 it 7.7 


F M 


3- 




33 


2101 326-1-7701 175 


L1228 S8 


52.0 


-I- 


2.9 


53.3 


-I- 


2.8 


50.1 


-I- 


2.5 


63.8 


it 


3.0 


78.2 


+ 


7.3 


74.3 it 9.8 


FIM 


3- 


15 


34 


21 01 32S-I-6S1 1 202 




3.07 


-|- 


0.16 


8.78 


-j- 


0.46 


14.8 




0.7 


23.3 


it 


1.3 


180. 


+ 


19. 




FI 






35 


2101351 -1-7703567 


1,1228 S4 


22.0 


-|- 


1.1 


22.6 


-|- 


1.1 


21.0 


-|- 


1.0 


22.1 


it 


1.0 


91.3 




8.5 


406. it 48. 


FIM 


3- 


15 


36 


2101358-1-6813257 




62.1 


-I- 


3.1 


66.5 


-(- 


3.4 


44.3 


-I- 


2.1 


61.9 


it 


3.0 


69.3 


± 


6.5 




FIM 






37 




T,1 C19 


4.77 


-[- 


0.23 


4.72 


-[- 


0.23 


4.38 


-|- 


0.22 


5.05 


it 


0.24 


22.5 


-j- 


2.1 


fiQ 7 -1- c « 
\jo. I HZ o.u 


FIM 


9. 

•Ji 


15 


38 


21 ni 42^-l-fi81 2'i7l 


NfOr! 7023 S N 


35.3 


-)- 


1.8 


29.4 


-1- 


1.6 


26.6 


-|- 


1.3 


28.3 


it 


1.4 


52.4 


-|- 


4.9 




FIM 


7- 


15 


39 


21 01 437+675026Q 




21.3 


±- 


1.0 


18.1 


-|- 


0.9 


14.4 


-|- 


0.7 


14.3 


it 


0.7 


20.7 


-|- 


1.9 




FIM 






40 


2101439+6814033 




31.6 


-I- 


1.5 


43.3 


+ 


2.2 


51.3 


-I- 


2.4 


46.9 


it 


2.2 


286. 


-I- 


27. 


266 + 36 


FIM 






41 


2101454+6804213 




13.5 


-I- 


0.7 


12.8 


+ 


0.7 


9.49 


-I- 


0.49 


17.9 


± 


0.9 


10.3 


-j- 


1.0 




FIM 






42 


2101467+6808453 


V* FU Oen- NGC 7023 S O 


42.1 


-I- 


2.0 


36.0 


-I- 


1.7 


25.3 


-I- 


1.2 


20.1 


it 


1.0 


35.4 


-I- 


6.0 




F M 


7- 


15 


43 


2101 496+7705479 


T,1228 S3 


18.3 


-|- 


0.9 


13.8 


-|- 


0.7 


12.4 


-|- 


0.6 


21.5 


it 


1.0 


133. 


-j- 


12. 


QCQ + S5 


FIM 




1 

'-•bo 


44 


2101527+6809519 


NGC 7023 S 44 


10.6 


-£ 


0.5 


13.4 


-£ 


0.7 


14.6 


-£ 


0.8 


20.4 


it 


1.2 


25.1 


-£ 


3.9 




FIM 


7 




45 


2102025+6811587 




1.81 


± 


0.09 


2.73 


± 


0.13 


2.95 


± 


0.21 


<6.79 


32.2 


± 


5.7 




F M 




1 


46 


2102050+7657185 


L1228 S9; [K98c] EM*43 


52.0 


± 


2.5 


40.1 


± 


1.9 


31.4 


± 


1.5 


39.8 


it 


1.9 


219. 




20. 


294. it 29. 


FIM 


3- 


8; 15 


47 


2102140+6813595 




12.6 




0.6 


11.2 


± 


0.6 


9.75 


-I- 


0.49 


10.2 


it 


0.5 


12.4 


± 


1.2 




FIM 






48 


91 021 41 +Rj^n7'?n^ 




7.76 


-|- 


0.40 


12.0 


J- 


0.6 


15.4 


-|- 


0.8 


27.1 


it 


1.6 


107. 


-)- 


10. 




Flm 






49 


2102212-1-6754203 


L1172 SMM 1 


13.8 


± 


0.7 


32.3 


it 


1.6 


43.1 


it 


2.1 


52.0 


it 


2.5 


252. 


it 


23. 


705. it 72. 


FI 


6 




50 


2102273+6754186 


L1172 SMM 2 


1.03 


± 


0.06 


1.15 


it 


0.07 


0.659 


it 


0.049 


0.356 


it 


0.033 


5.63 


it 


0.64 




F 


6 




51 


2102280+6808477 




4.91 


± 


0.24 


4.87 


± 


0.24 


4.84 


it 


0.24 


6.94 


± 


0.33 


21.4 


± 


2.0 




FIM 






52 


2102282+6803298 




51.1 


it 


2.6 


45.1 


± 


2.2 


44.7 


it 


2.7 


64.4 


it 


3.2 


92.7 


± 


8.6 


88.0 it 12.3 


FI 


15 




53 


2102299+6754083 




100.0 


± 


5.0 


120. 


it 


6. 


127. 


it 


6. 


130. 


it 


7. 


166. 


± 


15. 


249. it 35. 


FIM 






54 


2102324-1-7655100 




3.55 


± 


0.17 


2.96 


it 


0.14 


2.48 


it 


0.13 


2.14 


it 


0.11 


1.28 


± 


0.19 




Flm 






55 


2102330-1-6807292 


V* FW Cep 


40.1 


± 


2.0 


34.8 


it 


1.7 


30.6 


it 


1.5 


29.4 


it 


1.4 


14.6 


± 


1.4 




FIM 


15 




56 


2102549-1-6806210 




14.9 


± 


0.7 


14.4 


it 


0.7 


12.1 


± 


0.6 


11.9 


it 


0.6 


31.5 


± 


2.9 




FIM 






57 


2102595+6806322 


[K98b] 30; GN 21.02.4.02 


13.8 


± 


0.7 


18.5 


it 


0.9 


27.1 


it 


1.3 


78.6 


it 


3.9 


239. 


± 


22. 


87.5 it 19.0 


FIM 


8; 


15 


58 


2102596+6808119 




14.1 


± 


0.7 


10.1 


± 


0.5 


6.87 


it 


0.33 


4.24 


it 


0.22 


1.74 


± 


0.24 




F m 






59 


2103076-1-6808340 




27.2 


it 


1.3 


23.9 


it 


1.1 


20.5 


it 


1.0 


22.6 


it 


1.1 


32.4 


± 


3.0 




FIM 






60 


2103244-1-6759067 


V* EH Ccp; [K98b] 33; 
[K98c] EM*42 


310. 


it 


15. 


288. 


it 


20. 


280. 


it 


14. 


364. 


it 


18. 


435. 


± 


40. 


97.4 it 11.2 


FIM 


8; 


15 


61 


2103594H-6749296 




15.9 


it 


0.8 


14.6 


it 


0.7 


12.5 


it 


0.6 


12.4 


it 


0.6 


18.5 


± 


1.7 




FIM 






62 


2104156-1-6742464 




1020. 


it 


167. 


2850. 


it 


157. 


2780. 


it 


182. 


1780. 


it 


100. 


366. 


± 


34. 




FI 
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Id. 


Name/Position 


Literature Name^ 


3.6/im 


4.5/im 


5.8/im 


8.0/im 


24.0Mm 


70.0(Um 


Flags'" Ref. 




(SSTgbsJ) 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 





63 


2228030+6901 166 


L1221 IRS 1 




109. 


± 


7 




149. 


± 


1 


64 




1,1991 TRS 
[RR98] HH 363 VLA 4 


■J) 


0.567 


-j- 




033 


3.10 


-|- 


Q 


65 


2229333+7513160 






5.62 


± 





28 


10.5 


± 





66 


2229594+7514037 


[K93c] 4; [TW96] T3 




15.8 


± 





9 


21.8 


± 


1 


67 


2230318+7514094 


[TW96] T2 




0.134 


± 





017 


0.421 


± 





68 


Z^OXUUUT^ 1 U J.OO ( ^ 


[KQ3rl fi' [TWQfil Tia 
ixvyot^j \jj 1 J. vv yuj i- xa 




0.430 




Q 


029 


1.07 







69 


2234405+7517442 






30.2 


+ 


1 


5 


28.4 


+ 


1 


70 


99344S3-I-7'S1 3352 






8.97 







44 


7.30 







71 


2235008+751 5364 






22.0 




I 


I 


23.4 




I 


72 


2235025+751 7584 


fRD95l 3 




35.6 


+ 


1 


g 


43.1 


± 


2 


73 


2235042+7518206 






32.4 


+ 


1 


g 


32.0 




1 


74 


2235096+751 6094 






7.96 


+ 





40 


8.40 


+ 





75 


22351 41 +751 5026 






1.13 







06 


1.50 







76 


2235167+7518470 


[K93b] 39 




75.9 


± 


4 





74.2 


± 


3 


77 


2235254+7517561 






102. 




5 




108. 




5 


78 


2235266+751 6369 






34.5 


+ 


1 


g 


53.2 


+ 


2 


79 


2235272+7518018 






37.2 


+ 


1 


9 


34.2 


+ 


1 


80 


2236056+751 8325 


fKq3hl 2- [K93rl 
fK98bl 105 




78.2 


-|- 


4 




102. 




5 


81 




ixvyot/j -i-U, irvyouj J-U / 




54.8 


-|- 


2 




54.9 




2 


82 


2236591+7521207 






2.87 


± 





19 


2.38 


± 





83 


2237496+7504063 


[K93b] 2 




18.9 


± 


1 





13.3 


± 





84 


2238116+7512144 






8.52 


± 





42 


7.31 


± 





85 


2238152+7507203 






13.9 


± 





7 


11.8 


± 





86 


2238188+7511536 






109. 


± 


9 




154. 


± 


9 


87 


2238296+7514265 






8.62 


± 





43 


6.98 


± 





88 


2238405+7508412 






7.38 


± 





44 


6.76 


± 





89 


2238428+7511367 


L1251B IRS 4 




3.91 


± 





40 


7.96 


± 





90 


2238469+7511337 


L1251B IRS 1; [K93c] 

[K98b] 111 


14; 


30.9 


± 


2 





115. 




6 


91 


2238507+7510351 


L1251B IRS 8 




4.10 


± 





21 


3.50 


± 





92 


2238530+7511235 


L1251B IRS 2 




6.43 


± 





46 


10.2 


± 





93 


2239047+7511010 


L1251B IRS 9 




1.72 


± 





09 


1.63 


± 





94 


2239133+7512155 


L1251B IRS 16 




41.1 


± 


2 





112. 


± 


6 


95 


2239147+7507161 






6.84 


± 





35 


4.65 


± 






392. ± 20. 387. ± 26. 1940. ± 187. 6940. ± 641. FIM 4 



1 s 




-1- 


9^^ 


84 
o.o^ 


-1- 


n iQ 


A7 ^ 


+ 


4 4 


^OSO + 4RQ TTT 
ouou. in ^uy. r i 


4- 1 


5 


15.6 


± 


0.8 


18.2 


± 


0.9 


57.1 


± 


5.3 


FI 




2 


25.2 


-j- 


1.3 


27.3 




1.4 


277. 




26. 


FIM 


9; 12 


042 


0.315 




0.040 


0.162 




0.033 


5.10 




0.49 


F 


12 


Uo 


U. 1 DO 


+ 




n ^79 


+ 


n n^9 


1 fi'i 


+ 


24 


... F 


Q- 1 9 


4 


27.7 


-|- 


1.4 


35.6 


-|- 


1.8 


141. 


-|- 


13. 


FIM 






o.uo 


+ 


n '^9 


Q OA 


+ 




18 2 


+ 


1 7 


FIM 








+ 


1 2 


^^9 


+ 


1 .o 


104 


+ 


1 n 

-LU. 


FIM 








-j_ 


2 2 


47 fi 


+ 


9 


fi7 1 
U i . 1 


+ 


U.O 


FIM 


14 


5 


36.0 


-£ 


1.7 


40.8 


-|- 


2.0 


83.7 




7.8 


FIM 




41 


6.19 




0.32 


5.13 


-|- 


0.24 


12.9 


-1- 


1.2 


FIM 


1 


08 


1.69 


± 


0.11 


2.56 


± 


0.12 


20.4 


± 


1.9 


FI 


00 


8 


73.0 


± 


3.6 


85.9 


± 


4.2 


109. 


± 


10. 


FIM 


15 




121 


+ 

HI 


u. 


123 


+ 


o. 


7R 9 


+ 


7 7 


FIM 




a 
u 


oo.u 


+ 






+ 


9 7 


7Q 9 


+ 


R 1 


FIM 




7 


33.0 




1.6 


35.5 


-|- 


1.7 


48.8 


-|- 


5.1 


FIM 






145. 




7 


246. 


-|- 


12. 


544. 


-|- 


51. 


FIM 


8; 9; 15 


8 


58.1 


± 


2.8 


72.6 


± 


3.7 


210. 


± 


20. 


FIM 


8; 9; 15 


14 


2.11 


± 


0.13 


2.06 


± 


0.11 


2.47 


± 


0.32 


Flm 




7 


9.61 


± 


0.48 


8.06 


± 


0.39 


94.8 


± 


8.8 


FIM 


5; 9; 15 


37 


8.07 


± 


0.39 


9.46 


± 


0.46 


28.3 


± 


2.6 


FIM 


5 


6 


10.1 


± 


0.5 


9.48 


± 


0.45 


16.9 


± 


1.6 


FIM 


5 




227. 


± 


11. 


285. 


± 


27. 


321. 


± 


30. 


FIM 


5; 15 


35 


6.73 


± 


0.34 


8.28 


± 


0.42 


10.7 


± 


1.0 


FIM 


5 


35 


6.24 


± 


0.31 


6.42 


± 


0.31 


9.46 


± 


0.88 


FIM 


5 


47 


8.94 


± 


0.65 


13.2 


± 


0.7 


391. 


± 


37. 


Flm 


5 




269. 


± 


13. 


376. 


± 


20. 


2790. 


± 


274. 


Flm 


5; 8; 9 


17 


3.05 




0.16 


2.74 


± 


0.13 


3.66 


± 


0.38 


FIM 


5 


6 


12.0 


± 


0.7 


16.0 


± 


0.8 


463. 


± 


43. 


FIM 


5 


08 


1.63 


± 


0.09 


1.89 


± 


0.09 


3.48 


± 


0.36 


FI 


5 




206. 


± 


11. 


234. 


± 


11. 


343. 


± 


32. 


FI 


5 


23 


3.22 


± 


0.17 


1.99 


± 


0.10 


0.870 


± 


0.155 


F m 
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Id. 


Name/Position 


Literature Name'' 






4.5/im 


5.8/im 


8.0/im 


24.0//m 


70.0(um 


Flags'^ 




Ref. 




(SSTgbsJ) 








[mjy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 








96 




[K93b] 45; 


[K93c] 


ID, 


OQ Q _|_ 1 Q 


32.1 ± 1.6 


29.8 ± 1.4 


33.1 + 1.6 


fi7 9 -L ft 9 




FIM 




9; 12; 15 






[TW96] Hla 
























Q7 




[K93b] 46 






1 1 n -1- Q 


162. ± 13. 


245. ± 12. 


( . in lo. 


9^n -t- 9fi 




FIM 


u, 




98 










iQ 4. _l_ n Q 


18.3 ± 0.9 


19.1 ± 0.9 


1 c -L n Q 
-LO. u in u. y 


11.2 ± 1.0 




FIM 


5" 


15 


99 


99A1 047-4-7'^ 1 n4QR 








4 on _|_ n 91 


2.79 ± 0.14 


1.98 ± 0.11 


j-.o J. HZ u.u f 


o.oo zn u.'iu 




F M 
















2. 


IRAC YSO Candidates 














100 


2100224+6813042 








2.22 ± 0.14 


4.43 ± 0.24 


5.10 ± 0.25 


5.91 ± 0.29 


<23.0 




fl 






101 


2101030+6813072 








10.9 ± 0.5 


11.9 ± 0.6 


9.88 ± 0.47 


12.3 ± 0.6 


<16.2 




fl 






102 


2101102+6811332 








1.63 ± 0.08 


1.39 ± 0.07 


1.28 ± 0.07 


1.39 ± 0.09 


<3.39 




I 






103 


2101125+6810195 


NGC 7023 S I 






34.2 ± 1.7 


27.4 ± 1.4 


26.4 ± 1.3 


33.6 ± 1.8 


<55.2 




fl 


7 




104 


2101361+6808226 


NGC 7023 S E 






24.1 ± 1.2 


31.1 ± 1.5 


29.2 ± 1.6 


55.2 ± 4.4 


<89.3 




fl 


7; 


15 


105 


2228018+6901191 


L1221 IRS 2 






83.6 ± 4.5 


137. ± 9. 


324. ± 16. 


381. ± 23. 


<368. 




n 


4 




106 


2235323+7518377 








1.84 ± 0.10 


1.71 ± 0.09 


1.51 ± 0.10 


1.66 ± 0.08 


<2.20 




n 






107 


2238424+7511454 


L1251B IRS 3 


[K93b] 


44; 


15.5 ± 0.9 


12.5 ± 0.6 


11.0 ± 0.7 


10.5 ± 0.6 


<82.7 




fi 


5; 


12; 15 






[TW96 H2b] 
























108 


2238440+7511266 


L1251B IRS 5 






9.89 ± 0.50 


8.02 ± 0.39 


5.50 ± 0.30 


5.12 ± 0.28 


<29.4 




fi 


5 




109 


2238481+7511486 


L1251B IRS 6 






5.64 ± 0.29 


5.99 ± 0.29 


6.55 ± 0.35 


6.84 ± 0.35 


<10.9 




fi 


5 
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Id. Name/Position Literature Name'^ 3.6jUm 4.5/im 5.8jUm S.OjUm 24.0/im VO.OjUm Flags'" Ref. 

(SSTgbsJ) [mjy] [mjy] [mjy] [mjy] [mjy] [mjy] 



3. 2MASS-MIPS YSO Candidates 



1 1 n 
iiU 


zUoiiUU-|-D ( UUz ( o 












i.iiD ± U.zi 


■ • • m 




111 
ill 


onQOQO'7_i_f;'7i oon? 
ZxJoZoZ ( -rO 1 izzu ( 












OA n _l_ O Q 

Z4t.u :ji z.o 


■ ■ • m 




1 1 o 

LIZ 














n one; I n i fi i 
U.yUO it U.iD4 


• • • m 




1 1 Q 














o.ZZ it U.vSO 


- ■ - m 


















A nn -1- n Ai 

^.VU It U-4± 


OZ.O IC 0.4 ivi 




lie; 

iio 


ZUv5oOZ4+D / UZ400 




1 no -1- n nfi 
i.Uz it U.UD 




U.4C51 it U.U40 




i.U4 it U. i / 


• ■ • m 




liD 


zUoyu4t/ -|-O(oiiio 






o _i_ A 1 o 

^.00 it XJ.IZ 




i.oy it U.U/ 


i.OU it \j.ZZ 


■ ■ ■ m 




117 


2044098+6720237 


V* FI Cep 










59.7 ± 5.5 


m 




118 


2047218+6748507 












3.44 ± 0.44 


m 




119 


2048103+6803019 












3.27 ± 0.38 


m 




120 


2053476+6818039 












26.4 ± 2.5 


60.3 ± 7.2 M 


1 


121 


2054530+6820374 












26.2 ± 2.4 


M 


00 


122 


2056331+6814482 












3.25 ± 0.35 


54.5 ± 6.9 m 


1 


123 


2058467+7740256 


[K98c] EM*30 










115. ± 11. 


M 


8; 15 


124 


2100529+7703149 


L1228 S5; [K98b] 24; 

[K98c] EM*35 




55.4 ± 2.9 




38.4 ± 2.5 


95.9 ± 8.9 


64.8 ± 10.4 M 


3; 8; 15 


125 


2102204+6825241 


V* FV Cep; [K98cl Em* 38 




83.4 ± 4.1 




116. ± 6. 


104. ± 10. 


M 


8; 15 


126 


2103416+6823457 






3.84 ± 0.18 




2.93 ± 0.14 


3.88 ± 0.40 


M 




127 


2104282+7703171 












34.4 ± 3.2 


51.7 ± 6.4 M 




128 


2105593+7642272 




1.32 ± 0.06 




0.787 ± 0.049 




1.72 ± 0.19 


m 




129 


2156133+7658142 












4.68 ± 0.47 


33.8 ± 7.3 m 




130 


2210100+7525255 












0.805 ± 0.266 


m 




131 


2214393+7520584 












0.881 ± 0.154 


m 




132 


2217097+7645029 












1.63 ± 0.30 


m 




133 


2236029+7523260 












7.40 ± 0.70 


M 




4. Known YSO Candidates 


134 


2039062+6802150 


L1157 MM; [K98b] 7; HH 365 


1.34 ± 0.29 


4.87 ± 0.49 


3.60 ± 0.51 


2.29 ± 0.44 


80.3 ± 7.5 


10700. ± 1160. 


8 


135 


2045539+6757387 


V* PV Cep; [K98b] 10; 
[K98c] EM* 9; HH 215 














8 


136 


2101251+7706541 


L1228 Sl(b) 


3.67 ± 0.18 


2.95 ± 0.14 


2.27 ± 0.12 


1.37 ± 0.07 


<0.179 




3; 15 



Table 6 — Continued 



Id. 


Name/Position 


Literature Name^ 


3.6/jm 


4.5/im 


5.8^m 


S.O^im 


24.0/im 


70.0^*111 Flags'' Rcf. 




(SSTgbsJ) 




[mjy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 



2360. ± 271. 10700. ± 887. 7160. ± 406. 

13.1 ± 0.7 6.65 ± 0.46 
<0.00645 



137 2101271+6810380 NGC7023 RS 2 47.7 ± 2.5 43.2 ± 2.1 59.6 ± 3.2 <56.6 

138 2101369+6809479 HD 200775; V* V380 Cep; 

NGC 7023 S A 

139 2101436+6809361 NGC7023 RS S3; NGC 28.4 ± 1.4 18.8 ± 0.9 

7023 S D 

140 2101493+7657481 IRAS F21025+7645; <0.00552 <0.00546 

[K98b] 32 

141 2103121+6813012 IRAS F21025+6801; 

[K98b] 31 

142 2234119+7518099 [K93b] 1; GSC 04600-02450 ■ ■ ■ 91.3 ± 14.6 686. ± 35. 474. ± 31. 

143 2235234+7517077 [TW96] Nla; [K98b] 103; 373. ± 26. 1000. ± 114. 1860. ± 119. 

[K93c] 8 



<20.7 

<0.0264 
<0.140 
1440. ± 137. 



15 
7 

15 



9; 15 
8; 9; 12 



References. — The SIMBAD so urce prefix, where relevant , is listed after each ref ere nce in square brac kets (a lower case l ette r denotes a p articular table in that 
ref erence). 1 - IChapman et al.ll2009l. 2 -[Kauffman n et ah II2OO5 I. 3 - Padgett et al. ''2004 4 - 'Young et al.!l20 09l. 5 -iLee et ahlbood. 6 - IVisser et ahl H 

m 



8 - IKudIIi998| [K98] 



[TW96]. 13 - iRodn'guez 



9 



Kun fc Prus"tilir9 93 [K93], 10 _ 
Reipurthlll99d [RR98], 14 - ' 



■GvurB udagyan & Magakyan 1977 [ GM] , 11 - Rodriguez , 
Rosvick fc Davidgelll995l [RD95], 15 - iKun et al.il20oi 



Reipurthlliggj [RR96], 12 




''For c2d papers (refs 1-5) the literature name, where one is given, is taken as the core name followed by the index or name given to the source in the actual paper. 

''YSO identification flag. The source was identified as a YSO candidate by: F - the 5-Band method with a 24/^m detection, f - the 5-Band method without a 24/.im 
detection, I - the IRAC colour- magnitude method, M - the 2MASS/MIPS colour- magnitude method (greater than 1 mag away from the limits), m - the 2MASS/MIPS 



colour-magnitude method (less than 1 mag away from the limits) 
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Table 7. YSO Properties 



Index Region am. Qm Tboi Tboi ibol ^IR 

Class [K] Class [Lq] [Lq] 
(a) (b) (c) (d) (e) 



1 


L1152 


1.0 


I 


33. 





1.6 


0.060 


2 


L1152 


-0.31 


II 


160. 


I 


0.50 


0.34 


3 


L1152 


-0.50 


II 


270. 


I 


2.1 


1.6 


4 


L1147 


-0.089 


F 


630. 


F 


0.017 


0.016 


5 


L1148 


0.81 


I 


81. 


I 


0.44 


0.017 


6 


L1228N 


-1.2 


II 


1500. 


II 


0.023 


0.019 


7 


L1228N 


0.66 


I 


54. 





0.22 


0.037 


8 


L1228N 


0.76 


I 


35. 





0.13 


0.019 


9 


L1228N 


0.096 


F 


79. 


I 


2.3 


0.90 


10 


L1228N 


-0.61 


II 


750. 


F 


0.0052 


0.0052 


11 


L1228N 


-0.60 


II 


820. 


F 


0.093 


0.089 


12 


L1228N 


-0.36 


II 


670. 


F 


0.039 


0.038 


13 


L1228N 


-1.2 


II 


1600. 


II 


0.13 


0.099 


14 


L1228N 


-1.2 


II 


1200. 


II 


0.055 


0.051 


15 


NGC7023 


-0.49 


II 


530. 


F 


2.4 


1.8 


16 


NGC7023 


-0.55 


II 


380. 


F 


0.36 


0.26 


17 


NGC7023 B 


1.6 


I 


59. 





0.72 


0.12 


18 


NGC7023 B 


0.95 


I 


49. 





0.34 


0.12 


19 


NGC7023 B 


-0.87 


II 


840. 


F 


0.095 


0.095 


20 


NGC7023 B 


-1.3 


II 


970. 


II 


0.84 


0.72 


21 


NGC7023 


-2.3 


III 


2500. 


II 


0.34 


0.19 


22 


L1228 S 


-1.2 


II 


1600. 


II 


0.014 


0.011 


23 


NGC7023 A 


-1.1 


II 


1600. 


II 


0.099 


0.074 


24 


NGC7023 A 


-1.0 


II 


1500. 


II 


0.17 


0.14 


25 


NGC7023 A 


-0.83 


II 


1400. 


II 


0.24 


0.18 


26 


NGC7023 A 


-0.91 


II 


1200. 


II 


0.027 


0.024 


27 


NGC7023 A 


-0.54 


II 


1000. 


II 


0.73 


0.64 


28 


L1228 S 


-0.24 


F 


590. 


F 


0.011 


0.011 


29 


NGC7023 A 


-0.62 


II 


780. 


F 


0.082 


0.082 


30 


L1228 S 


-1.1 


II 


1700. 


II 


0.92 


0.59 


31 


NGC7023 A 


-0.61 


II 


870. 


F 


0.11 


0.10 


32 


L1228 S 


-0.50 


II 


730. 


F 


0.21 


0.13 


33 


L1228 S 


-0.77 


II 


1000. 


II 


0.20 


0.15 


34 


NGC7023 A 


1.0 


1 


14. 





0.74 


0.060 


35 


L1228 S 


-0.44 


II 


53. 





0.15 


0.066 


36 


NGC7023 A 


-0.91 


II 


1100. 


II 


0.35 


0.32 


37 


L1228 S 


-0.31 


II 


340. 


1 


0.018 


0.013 


38 


NGC7023 A 


-0.91 


II 


1400. 


11 


0.32 


0.24 


39 


L1172 


-1.1 


II 


1200. 


II 


0.13 


0.12 


40 


NGC7023 A 


0.37 


I 


310. 


I 


0.27 


0.16 


41 


NGC7023 


-1.2 


II 


1300. 


II 


0.11 


0.10 


42 


NGC7023 A 


-1.2 


II 


1500. 


II 


0.34 


0.27 


43 


L1228 S 


-0.29 


F 


560. 


F 


0.15 


0.083 
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Table 7 — Continued 



Index Region Tboi Ttoi i/bol -t-iR 

Class [K] Class [Lq] [Lq] 
(a) (b) (c) (d) (e) 



44 


NGC7023 A 


-0.50 


II 


830. 


F 


0.069 


0.069 


45 


NGC7023 A 


0.023 


F 


1000. 


II 


0.044 


0.027 


46 


L1228 S 


-0.46 


II 


760. 


F 


0.32 


0.19 


47 


NGC7023 


-1.0 


II 


1200. 


II 


0.077 


0.071 


48 


NGC7023 


0.53 


I 


390. 


F 


0.067 


0.065 


49 


L1172 


1.1 


I 


42. 





0.58 


0.13 


50 


L1172 


-0.29 


F 


24. 





0.24 


0.0039 


51 


NGC7023 


-0.36 


II 


730. 


F 


0.032 


0.031 


52 


NGC7023 


-0.62 


II 


920. 


F 


0.76 


0.13 


53 


L1172 


-0.55 


II 


210. 


I 


0.53 


0.43 


54 


L1228 S 


-1.5 


II 


1600. 


II 


0.015 


0.011 


55 


NGC7023 


-1.5 


II 


2000. 


II 


0.44 


0.27 


56 


NGC7023 


-0.53 


II 


730. 


F 


0.057 


0.057 


57 


NGC7023 


0.51 


I 


500. 


F 


0.29 


0.17 


58 


NGC7023 


-2.1 


III 


2000. 


II 


0.16 


0.11 


59 


NGC7023 


-1.00 


II 


1300. 


II 


0.21 


0.17 


60 


L1172 


-0.81 


II 


1200. 


II 


3.3 


2.1 


61 


Off-cloud 


-0.97 


II 


1500. 


II 


0.13 


0.100 


62 


L1172 


-2.0 


III 


1800. 


II 


45. 


37. 


63 


L1221 


0.84 


I 


66. 





3.0 


0.73 


64 


L1221 


1.0 


I 


21. 





1.4 


0.0097 


65 


L1251W 


0.12 


F 


95. 


I 


0.54 


0.042 


66 


L1251W 


0.49 


I 


130. 


I 


1.3 


0.11 


67 


L1251W 


0.71 


I 


6.5 





0.048 


0.0012 


68 


L1251W 


-0.51 


II 


9.3 





0.017 


0.0016 


69 


L1251A 


-0.29 


F 


680. 


F 


0.19 


0.18 


70 


L1251A 


-0.81 


II 


1400. 


II 


0.10 


0.081 


71 


L1251A 


-0.23 


F 


620. 


F 


0.13 


0.13 


72 


L1251A 


-0.54 


II 


720. 


F 


0.18 


0.18 


73 


L1251A 


-0.48 


II 


770. 


F 


0.19 


0.18 


74 


L1251A 


-0.85 


II 


930. 


F 


0.041 


0.041 


75 


L1251A 


0.39 


I 


330. 


I 


0.0087 


0.0087 


76 


L1251A 


-0.83 


II 


1200. 


II 


0.60 


0.49 


77 


L1251A 


-1.1 


II 


1100. 


II 


0.77 


0.71 


78 


L1251A 


-0.28 


F 


610. 


F 


0.17 


0.17 


79 


L1251A 


-0.91 


II 


1200. 


II 


0.30 


0.25 


80 


L1251A 


0.033 


F 


470. 


F 


0.76 


0.66 


81 


L1251A 


-0.46 


II 


910. 


F 


0.48 


0.43 


82 


L1251A 


-1.1 


II 


1200. 


II 


0.016 


0.016 


83 


L1251B 


-0.56 


II 


1200. 


II 


0.35 


0.20 


84 


L1251B 


-0.44 


II 


810. 


F 


0.055 


0.052 


85 


L1251B 


-0.90 


II 


1200. 


II 


0.092 


0.077 


86 


L1251B 


-0.69 


II 


880. 


F 


1.1 


1.0 
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Table 7 — Continued 



Index Region aiR aiR Ttoi Ttoi ibol ^IR 

Class [K] Class [Lq] [Lq] 
(a) (b) (c) (d) (e) 



87 


L1251B 


-0.99 


11 


1400. 


11 


0.088 


0.071 


88 


L1251B 


-0.96 


II 


1400. 


II 


0.072 


0.056 


89 


L1251B 


1.2 


I 


35. 





1.3 


0.075 


90 


L1251B 


1.6 


I 


49. 





16. 


0.99 


91 


L1251B 


-1.1 


II 


1300. 


II 


0.030 


0.026 


92 


L1251B 


0.97 


I 


61. 





1.6 


0.099 


93 


L1251B 


-0.68 


11 


900. 


F 


0.010 


0.010 


94 


L1251B 


0.94 


I 


470. 


F 


0.50 


0.50 


95 


L1251B 


-2.1 


III 


2100. 


II 


0.098 


0.063 


96 


L1251B 


-0.82 


II 


1400. 


II 


0.39 


0.29 


97 


L1251B 


-0.72 


II 


840. 


F 


1.4 


1.2 


98 


L1251B 


-1.1 


II 


920. 


F 


0.087 


0.087 


99 


Off-cloud 


-1.3 


11 


2100. 


II 


0.071 


0.043 


100 


NGC7023 A 


0.14 


F 


27. 





0.16 


0.0053 


101 


NGC7023 A 


-0.95 


11 


1700. 


II 


0.086 


0.050 


102 


NGC7023 A 


-1.4 


II 


1800. 


II 


0.016 


0.011 


103 


NGC7023 A 


-1.3 


II 


1700. 


II 


0.29 


0.22 


104 


NGC7023 A 


-0.052 


F 


690. 


F 


0.039 


0.039 


105 


L1221 


1.0 


I 


640. 


F 


0.20 


0.20 


106 


L1251C 


-1.3 


11 


1500. 


11 


0.014 


0.012 


107 


L1251B 


-1.7 


III 


2200. 


II 


0.24 


0.14 


108 


L1251B 


-1.7 


III 


1600. 


II 


0.068 


0.056 


109 


L1251B 


-0.69 


II 


1000. 


II 


0.023 


0.023 


110 


Off-cloud 


-1.3 


II 


3300. 


III 


0.091 


0.069 


111 


Off-cloud 


-0.016 


F 


890. 


F 


0.070 


0.061 


112 


Off-cloud 


-1.9 


111 


4200. 


111 


0.32 


0.19 


113 


Off-cloud 


-0.75 


11 


2100. 


11 


0.054 


0.047 


114 


L1147 


-1.2 


II 


1000. 


II 


0.22 


0.17 


115 


L1147 


-1.2 


II 


1600. 


II 


0.012 


0.0100 


116 


Off-cloud 


-1.7 


III 


3500. 


III 


0.17 


0.068 


117 


Off-cloud 


-2.1 


III 


2000. 


II 


12. 


11. 


118 


Off-cloud 


-0.74 


11 


1400. 


11 


0.049 


0.049 


119 


Off-cloud 


-1.8 


111 


2500. 


11 


0.70 


0.61 


120 


Off-cloud 


-1.1 


11 


1800. 


11 


0.88 


0.75 


121 


Off-cloud 


-1.0 


II 


2100. 


II 


0.62 


0.55 


122 


Off-cloud 


-1.0 


II 


360. 


F 


0.070 


0.067 


123 


L1228N 


-0.89 


II 


1600. 


II 


0.89 


0.85 


124 


L1228 S 


-0.93 


II 


1100. 


II 


0.31 


0.22 


125 


NGC7023 


-0.82 


11 


1100. 


11 


0.68 


0.51 


126 


Off-cloud 


-1.2 


II 


1800. 


11 


0.047 


0.036 


127 


Off-cloud 


-0.94 


II 


1500. 


II 


0.34 


0.30 


128 


Off-cloud 


-1.0 


II 


1800. 


II 


0.0070 


0.0053 


129 


L1241 


-2.1 


III 


3000. 


III 


2.1 


1.6 
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Table 7 — Continued 



Index 


Region 




aiR 






^bol 


^IR 








Class 


[K] 


Class 


[^0] 


[Lq] 






(a) 




(b) 


(c) 


(d) 


(e) 


130 


Off-cloud 


-1.4 


II 


2100. 


II 


0.048 


0.048 


131 


Off-cloud 


-1.6 


II 


3300. 


III 


0.11 


0.082 


132 


Off-cloud 


-1.0 


II 


1800. 


II 


0.039 


0.039 


133 


L1251C 


-0.91 


II 


1600. 


II 


0.14 


0.13 


134 


L1157 


0.94 


I 


37. 





8.9 


0.021 


135 


Off-cloud 


0.66 


I 


100. 


I 


46. 


23. 


136 


L1228 S 


-1.7 


III 


1200. 


II 


0.0056 


0.0056 


137 


NGC7023 A 


-0.80 


II 


1400. 


II 


0.23 


0.20 


138 


NGC7023 A 


-1.1 


II 


2400. 


II 


100. 


43. 


139 


NGC7023 A 


-2.6 


III 


1900. 


II 


0.27 


0.27 


140 


L1228 S 














141 


NGC7023 














142 


L1251A 


-0.29 


F 


500. 


F 


4.1 


3.1 


143 


L1251A 


2.4 


I 


79. 


I 


33. 


8.3 



^The infrared spectral index, ajR, calculated by a least-squares fit between 
2MASS Ks and 24 fim. Where IRAC data were not available qir was calculated 
from the ii'a-[24] colour as described in Section 14.21 

''The bolometric temperature, T^oi, calculated using all available data from 
B-band to 1200 fim. 

'^The equivalent class calculated from T^oi (see Section 15.21 1. 

''The bolometric luminosity, Lboli calculated using all available data from B- 
band to 1200 fim. 

''The infrared luminosity, Lir, calculated using all available data from 1 to 
30 fim. 
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16, 12, 60, and 12 
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Table 8. Class statistics of Spitzer identified YSOs 



Class 


No. 


of YSOs 






aiR 




c2d^ 







17 ( 13%) 




I 


21 ( 16%) 


10 ( 8%) 


16% 


"Flat" 


14 ( 11%) 


37 ( 28%) 


12% 


II 


87 ( 65%) 


64 ( 48%) 


60% 


III 


11 ( 8%) 


5 ( 4%) 


12% 



^To tal YSQ sta.tistics for all c2d re- 
gions (jEvans et alj|2008l ). 
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Table 9. Additional Fluxes 



Id. 


GSC-II Name 


B V 


R 


/ 


IRAS Name 


12.5/um 


25/im 


60/um 


100;um 






[mjy] [mjy] 


[mJy] 


[mJy] 




[Jy] 


[Jy] 


[Jy] 


[Jy] 


1 










t z0o5o+d742 


J' n neon 
< 0.0522 


0.302 


4.00 


7.00 


2 






0.109 


0.641 












3 










1 20358+6746 


0.447 


1.000 


1.68 


< 5.71 


4 






n nono 


0.185 












5 










r 20404+0712 


< 0.0636 


0.126 


0.320 


< 4.39 


6 


INU Y ZUU ( zoo 




U.Uo7z 


0.715 












7 

8 


JNU Y4UUoo44 
JNU Y4UUoo4i 




1.62 
0.0215 


6.17 












9 


"NTnx/ /I nnonnT 

JNU Y4UUozy7 




n no A 

0.0848 


0.652 


r 20582+7724 


1.19 


3.23 


11.2 


17.1 


11 


ATnA/" A nno o n 

JN U Y 4UUoo4U 




n nooo 


0.324 












12 


INU Y4UlzoDi 




0.0497 


0.190 












13 


INU Y4UUozoo 




1.65 


T 1 K 

7.15 












14 


JNUY4U0oUzo 




n nc A A 

0.0d44 


0.477 












15 


NlKbulilyb 


0.278 7.13 


5.76 


51.3 


b 20587+6803 


0.558 


0.952 


0.822 


< 13.4 


16 


M1 T./'oni nr^nn 

JN lKbUlU799 


0.426 


2.45 


2.10 












17 










b 20597+6801 


< 0.136 


0.712 


< 7.33 


< 88.6 


20 


JNIK0UII77O 


n nnoo a n 1 00 

0.U0234 U.188 


2.27 


7.92 












21 


■\Ti T/'Onnm nn 

JNlKbU09199 


0.408 2.84 


14.9 


20.1 












22 


JNUY 4005594 




0.0693 


0.652 












23 


JNlKbUll4Ul 


0.217 


0.308 


2.72 












24 


"\T1 Ti"" n r- n 






4.95 












25 


M1 T/"C"ni "1 oni 

JN iKbUiloUi 


0.361 


1.10 


8.06 












26 


JN iKbUil72D 




0.0245 


0.143 












OT 

11 


iNlKoUoUoUo 




0.0384 














30 


■\Tn'\7' A nnno 00 

JNUY4UUUZOO 




52.1 


97.8 












31 


M1 Tx"" n /1 1 1 

JN iKbUio411 






0.293 












32 


Mnv /I nn /I /1 Q 
IN U Y 4UU44oo 




5.16 


14.7 


b 21023+ /60O 


0.171 


0.225 


< 1.09 


< 25.9 


33 


NUY4uU4o7d 




5.45 


14.0 












35 


N0Y4004980 




0.325 


1.39 












36 


N1KS011820 


0.198 


1.08 


3.18 












37 


N0Y4013403 






0.145 












38 


N1KS026404 


1.25 




10.2 












39 


N1KS004876 




0.0760 


0.763 












40 


N1KS011945 


0.198 


1.79 


1.09 












41 


N1KS009227 


0.130 


0.344 


0.901 












42 


N1KS026141 






14.9 












43 


N0Y4005414 




3.02 


16.5 












45 


N1KS026373 






5.09 












46 


NOY4000358 




15.0 


28.7 


F21028+7645 


0.0847 


0.198 


< 0.283 


< 4.63 


47 


N1KS017685 




0.0462 


0.272 












48 


N1KS026025 






0.183 












49 










F21017+6742 


< 0.136 


0.296 


1.31 


4.76 


51 


N1KS036074 




0.0350 














52 


N1KS008923 


0.295 1.81 


12.4 


28.5 
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Table 9 — Continued 



Id. 


GSC-II Name 


B 


V 


R 


/ 


IRAS Name 


12.5//m 


25/im 


60;im 


100/im 






[mjy] 


[mJy] 


[mJy] 


[mJy] 




[Jy] 


[Jy] 


[Jy] 


[Jy] 


54 








0.126 


0.821 












55 


M "1 jy c ri r\ ri o o 

IN IKbUUUzzo 


0.307 


4.67 


15.7 


22.0 












57 


M 1 jy o n r» o o c 

IN IKbUUyooo 


0.313 


3.00 


5.60 


4.81 












58 


iVTi TyOAi r\Ar\ci 
In iJxoUiU4UD 


0.0124 


0.259 


1.63 


4.81 












59 


JN IKbU 10491 


0.0276 


0.388 


2.17 


3.69 












60 


JN IKbUUUozo 


5o.4 


48.0 


108. 


49.0 


t 21U2y+d747 


0.562 


0.690 


< 0.520 


< 4.98 


61 


N1KS004279 


0.140 


0.338 


2.45 


3.55 












62 


M "1 jy o ri ri ri r\ 

IN IKbUUDbUU 


1.16 


12.2 


44.1 


533. 


r 210o5H-o7oU 


1.35 


0.562 


< 0.315 


< 7.01 


66 












r zzzyUH- ( 4oo 


< u.uy7z 


0.333 


0.774 


< 10.6 


69 


JN0XD010892 






0.0234 


0.239 












70 


JNUXDUlUU42 


r\ r\noi A 
U.UU0I4 


0.131 


0.147 


1.70 












76 


AT n "V" T~\ m 0'7'7 Q 






0.0271 


0.176 












To 


IN UAUUl ilUo 


0.332 


1.86 


1.45 


15.8 












77 


MnvT^m new 


n r\'7'i T 

U.UYoY 


0.300 


0.291 


1.92 












79 


iNuADUiUoTD 


0.0162 


0.348 


0.576 


5.04 












80 


MnvT^m nm o 


0.142 


0.754 


2.81 


13.3 


b 22350+7502 


0.317 


0.598 


< 7.05 


< 71.8 


81 


TVT n "V T~\ n o n n n 


n nnono 

0.00808 


0.182 




3.65 












83 


JNuXUUUdoIy 


0.818 


5.77 


7.88 


18.1 


T^t^c^lCT 1 ''7 A AO 

r 22cSo7+y44o 


< 0.0783 


0.132 


< 0.289 


< 2.65 


84 


JNUXDUzo54o 








0.220 












85 


JN UXDUU7o4o 


0.0361 


0.195 


0.350 


1.61 












86 


iN UAJJUUyUoo 


u.U4yu 


0.726 


2.03 


8.52 


1* z2o7z+74od 


0.178 


< 0.247 


< 21.5 


< 89.7 


87 
88 


iVTn'VT^nAn'zno 
iNUAJJUUy i\)Z 

JNUADUUoOoY 


0.0116 

0.0101 


0.140 
0.162 


0.194 
0.187 


1.67 
1.41 












90 












1 22376+7455 


0.693 


5.47 


31.1 


67.6 


91 


IN UAUUUodU / 






0.0205 


0.285 












92 


IN U ADU23537 








0.277 












95 


iNUA_DU0746D 


0.0476 


0.381 


0.712 


4.68 












96 


1nUA_Du0o510 


0.380 


2.36 


3.54 


11.2 












97 


N UA.Du09o47 


0.0598 


1.68 


5.76 


19.3 


b 22385+7457 


0.309 


0.330 


< 1.03 


< 29.3 


99 


N0XD0O8377 


0.0445 


0.627 


0.535 


3.86 












101 


N1KS011801 




0.361 


1.10 


8.06 












102 


N1KS024073 






0.0484 


0.742 












103 


N1KS026819 








13.8 












106 


N0XD023601 








0.0917 












107 


N0XD009004 


0.343 


2.21 


4.88 


14.0 












108 


N0XD008908 






0.157 


0.978 












110 


N1JU019152 


0.242 


1.01 


1.60 


1.81 












111 


N1KP005099 


0.383 


0.512 


0.525 


0.696 












112 


N1KP000227 


3.09 


6.62 


9.74 


10.4 












113 


N1KR012754 






0.930 


0.735 












114 


N1KP016685 


0.424 


1.16 


2.49 


4.56 












115 


N1KP030229 






0.0252 


0.139 












116 


N1KP010819 


1.61 


3.88 


6.80 


10.2 












117 


N1KP012225 




0.300 


7.25 




F20436+6709 


0.212 


0.0681 


< 0.221 


< 4.05 
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Table 9 — Continued 



Id. 


GSC-II Name 


B 


V 


R 




IRAS Name 


12.5//m 


25/um 


60/xm 


lOO/xm 






[mJy] 


[mJy] 


[mJy] 


[mJy] 




[Jy] 


[Jy] 


[Jy] 


[Jy] 


119 


N1KQ003034 






5.71 














120 


N1KQ008280 


0.187 


0.706 


9.83 


17.6 












121 


N1KQ009201 


0.0680 


0.842 


2.91 


7.92 












123 


N0Y4008425 






2.52 


7.92 












124 


NOY4000284 






8.56 


19.9 


F21016+7651 


< 0.0966 


0.149 


< 0.669 


< 19.2 


125 


N1KS019292 


0.682 


2.23 


19.8 


19.9 


F21017+6813 


0.119 


0.140 


< 0.357 


< 7.68 


126 


N1KS019057 






0.598 


1.29 












127 


N0Y4004915 






3.23 


9.69 












128 


N0Y7009266 






0.0598 


0.390 












129 


N0XW000215 






46.6 


73.5 












131 


N0X6013952 


0.196 


1.48 


1.79 


2.85 












133 


N0XD012039 






0.0425 


0.419 












134 












F20386+6751 


< 0.0656 


0.226 


9.97 


42.0 


135 


N1KR002065 


0.0871 


0.150 


10.7 


4.51 


F20453+6746 


12.9 


30.6 


46.4 


50.0 


137 


N1KS030305 






0.0384 














138 


N1KS000207 


3430. 


3850. 
















142 


N0XD000257 


2.55 


10.8 


47.5 


37.9 


F22331+7502 


1.16 


1.74 


1.84 


< 21.9 


143 












F22343+7501 


5.74 


27.4 


61.1 


77.9 
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Table 10. Color-color transforms to the Cousins-Johnson System 



System 




Transforms 


Reference 


Quick 


V-- 


~^Quick 




Tycho 


V-- 


=Vt - 0.090(5t - Vt) 


The Hipparcos and Tycho 




B-- 


=0.850(Bt -Vt)-V 


Catalogs fESA 19971 


POSS-II 


B- 


={Bj - 0.28y)/0.72 


Blair & Gilmore 1982 




R- 


=Rf 






I- 


=/n 
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Table 11. MIPS 160 ^um, SCUBA 450 & 850 ^um, and MAMBO 1200 /im photometry of Cepheus 

YSO candidates 



YSO # Region 


160/im 


450^m 


850/xni 


UOOfim 


Ref.'^ 




[Jy] 


[Jy] 


[Jy] 


[Jy] 





1 


L1152 


8.8 


± 





3 


8.3 


± 


4.4 


1.2 


± 





2 






1 




2 


L1152 


1.6 


± 





1 
























3 


L1152 


1.6 


± 





1 
























5 


L1148 
























0.059 ± 


007 


4 




7 
8 


L1228N 
L1228N 


3.9 
3.9 


± 
± 






1 
1 
























9 


L1228N 


15 


± 


1 




1.1 


± 


1.4 


2.0 


± 





3 


0.89 ± 


01 


1, 


4 


17 


L1174 


4.2 


± 





1 
























18 


L1174 


4.2 


± 





1 
























34 


L1174 










9.4 


± 


0.4 


1.4 


± 





1 










35 


L1228S 


6.2 


± 





3 
























49 


L1172 


1.7 


± 





1 


3.1 


± 


1.5 


0.35 


± 





07 






3 




50 


L1172 


1.7 


± 





1 








0.32 


± 





06 


0.24 ± 


01 


3, 


4 


53 


L1172 


1.7 


± 





1 
























63 


L1221 










10 


± 


6 


2.0 


± 





2 






1 




64 


L1221 










11 


± 


6 


2.0 


± 





3 






1 




65 


L1251W 
























0.35 ± 


01 


4 




66 


L1251W 
























0.21 ± 


01 


4 




67 


L1251W 
























0.46 ± 


01 


4 




68 


L1251W 
























0.36 ± 


01 


4 




89^^ 


L1251B 










2.1 


± 


1.1 


0.58 


± 





07 






1, 


2 


90'' 


L1251B 










21 


± 


11 


5.8 


± 





7 






1, 


2 


92'' 


L1251B 










1.3 


± 


0.7 


0.34 


± 





04 






1, 


2 


100 


L1152 


4.2 


± 





1 
























134 


L1157 










13 


± 


7 


2.5 


± 





5 






1 




135 


PV Cep 










12 


± 


1.0 


1.5 


± 





1 










139 


L1251A 










16 


± 


9 


2.3 


± 





3 






1 





References 
([i002) 



1 



4 - Kauffmann et al 



(|2008D 



Young et all (|2006l ). 2 - iLee et al.l (|2006|), 3 - IVisser et al 



The 850 and 450 /x m aper ture radius was 20" for I Young et al.l (|2006r ) and 25" 



for Visser et al. 



( 2002h . The lYoung et all (|2006D flux for L1251B was portioned 
between the sources by Lee et al. ( 20061 ) based on the relative brightness of the 
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sources at 24 /im. The iKaufFmann et al.l (j2008l ) 1200 /im fluxes were measured in 
a 20"radius aperture. The 160 fim fluxes were measured in a 40" radius aperture. 

^The errors for these sources a re taken as the same as the relative errors for the 
combined source, L1251E-1, from lYoung et all |2006h . 
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Table 12. Class O/I derived submillimcter parameters 



Id. 


^dust 


-^submm 


-^submm /-^'bol 




[K] 




[%] 


1 


17. 


0.080 


5.0 


49 


15. 


0.054 


9.4 


63 


18. 


0.13 


4.3 


64 


17. 


0.13 


9.1 


90 


16. 


0.50 


3.1 


134 


18. 


0.25 


2.8 


135 


19. 


0.17 


0.37 


143 


19. 


0.23 


0.69 
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Table 13. SED modelling parameters for Class II YSOs 



Index Ay 


Anorm 




log(L(jiac/i'star) -^turnoft ^excess 








[Mm] 



2 


2.3 


J 


0.31 


0.20 


1.7 


-0.28 


3 


9.5 


J 


3.5 


-0.29 


3.6 


-0.55 


6 


0.50 


J 


0.030 


-0.76 


5.7 


-1.0 


10 


3.8 


Ks 


0.0040 


-0.045 


3.6 


-0.70 


11 


9.0 


J 


0.29 


-0.38 


2.2 


-0.87 


12 


13. 


J 


0.099 


1.1 


2.2 


-0.74 


13 


3.0 


J 


0.25 


-0.93 


5.7 


-0.86 


14 


7.7 


J 


0.21 


-0.67 


3.6 


-1.3 


15 


0.42 


J 


1.3 


0.16 


2.2 


-0.72 


16 


1.5 


J 


0.33 


-0.084 


2.2 


-0.26 


19 


19. 


■J 


0.52 


-0.58 


1.7 


-1.5 


20 


2.0 


J 


1.2 


-0.31 


2.2 


-0.97 


22 


3.2 


J 


0.032 


-1.5 


7.9 


0.16 


23 


1.5 


J 


0.16 


-1.0 


5.7 


-0.81 


24 


2.7 


J 


0.28 


-0.80 


5.7 


-0.77 


25 


3.6 


J 


0.50 


-0.44 


5.7 


-0.38 


26 


5.5 


J 


0.075 


-0.93 


3.6 


-0.89 


27 


3.9 


J 


1.3 


-0.70 


3.6 


-0.41 


29 


15. 


J 


0.28 


-0.49 


2.2 


-1.1 


30 


0.19 


J 


0.85 


-0.49 


5.7 


-0.80 


31 


10. 


J 


0.41 


-0.15 


3.6 


-0.67 


32 


1.4 


J 


0.20 


-0.34 


7.9 


-0.35 


33 


0.26 


J 


0.14 


-0.078 


3.6 


-0.87 


35 


1.8 


J 


0.087 


0.15 


2.2 


0.18 


36 


1.9 


J 


0.40 


-0.17 


2.2 


-1.0 


37 


2.2 


J 


0.013 


0.096 


1.7 


-0.32 


38 


0.93 


J 


0.40 


-0.55 


3.6 


-0.76 


39 


6.2 


J 


0.42 


-1.00 


3.6 


-1.1 


41 


1.5 


J 


0.16 


-0.43 


3.6 


-1.1 


42 


1.1 


J 


0.42 


-0.41 


3.6 


-1.1 


44 


4.8 


J 


0.068 


-0.14 


2.2 


-0.64 


46 


0.39 


J 


0.29 


-0.43 


5.7 


-0.043 


47 


8.0 


J 


0.29 


-0.57 


3.6 


-1.1 


51 


11. 


J 


0.13 


0.71 


3.6 


-0.35 


52 


2.3 


I 


2.3 


-1.2 


5.7 


-0.79 


53 


20. 


J 


1.7 


-0.24 


1.7 


-0.82 


54 


1.8 


J 


0.023 


-0.99 


5.7 


-1.5 


55 


0.48 


J 


0.46 


-0.44 


5.7 


-1.6 


56 


14. 


J 


0.17 


-0.43 


2.2 


-0.95 


59 


1.2 


J 


0.30 


-0.54 


3.6 


-0.90 


60 


0.18 


J 


2.8 


-0.22 


3.6 


-1.2 


61 


2.8 


J 


0.24 


-0.50 


3.6 


-0.94 


68 


5.1 


IRAC2 


0.0037 




7.9 
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Table 13 — Continued 



Index 




■^norm 


^star 


log(-t'disc/-t'star) 


•^turnoff 


Q^excess 






[fMUl] 


[i©] 




[/im] 




YU 


3.1 


J 
J 


0.22 


-1.1 


5.7 


-0.35 


72 


16. 


J 


0.51 


-0.31 


2.2 


-1.0 


73 


5.8 


J 


0.32 


-0.23 


1.7 


-0.77 


74 


14. 


J 


0.18 


-0.83 


3.6 


-0.99 


76 


0.38 


J 


0.46 


-0.075 


2.2 


-0.88 


77 


6.1 


J 


2.3 


-0.63 


3.6 


-1.3 


79 


4.4 


J 


0.71 


-0.47 


3.6 


-0.90 


81 


3.5 


J 


0.89 


-0.51 


3.6 


-0.33 


82 


6.2 


J 


0.031 


-0.70 


3.6 


-1.1 


83 


1.7 


J 


0.49 


-0.39 


7.9 


1.3 


84 


8.4 


J 


0.16 


-0.24 


3.6 


-0.44 


85 


1.2 


J 


0.12 


-0.48 


2.2 


-0.98 


86 


2.0 


J 


1.5 


-0.25 


3.6 


-0.61 


87 


3.3 


J 


0.19 


-0.90 


5.7 


-0.74 


88 


2.1 


J 


0.13 


-0.92 


0.7 


-0.72 


91 


4.1 


J 


0.072 


-0.98 


3.6 


-1.1 


93 


9.9 


J 
J 


0.025 


-0.60 


3.6 


n T7 
-0.77 


96 


3.1 


J 
J 


0.73 


-0.36 


0.7 


-0.45 


97 


1.9 


J 


1.9 


-0.14 


3.6 


-0.64 


98 


18. 


7 

J 


0.43 


-O.oU 


1 ^7 

1.7 


-1.6 


99 


0.66 


7 
J 


0.095 


-1.3 


7.9 


-0.027 


101 


7.2 


J 


0.22 


1.1 


3.6 


-1.1 


102 


1.6 


J 


0.025 


-1.1 


5.7 


-0.72 


103 


3.4 


J 


0.55 


-0.44 


5.7 


-0.20 


106 


3.0 


7 
J 


0.027 


-1.1 


0.7 


-0.67 


109 


11. 


7 

J 


0.079 


-0.75 


2.2 


-1.3 


110 


0.00 


J 


0.099 


-1.4 


2.2 


-1.3 


113 


0.49 


7 
J 


0.058 


-0.86 


2.2 


-0.79 


114 


0.20 


J 


0.23 


-1.2 


2.2 


-0.48 


115 


5.6 


J 


0.039 


-0.59 


5.7 


-0.52 


118 


2.2 


J 


0.068 


-2.0 


2.2 


-0.82 


120 


2.4 


J 


1.3 


-0.66 


2.2 


-0.98 


121 


1.3 


J 


0.86 


-2.3 


2.2 


-1.1 


122 


2.7 


J 


0.099 


-1.4 


2.2 


-0.41 


123 


1.4 


J 


0.94 


-0.42 


2.2 


-0.98 


124 


0.96 


J 


0.35 


-0.41 


2.2 


-0.95 


125 


1.1 


J 


0.70 


-0.059 


2.2 


-0.86 


126 


1.9 


J 


0.075 


-1.1 


4.5 


-0.99 


127 


1.5 


J 


0.45 


-1.6 


2.2 


-0.90 


128 


1.3 


J 


0.011 


-1.2 


5.7 


-0.46 


130 


2.0 


J 


0.065 


-2.7 


2.2 


-1.5 


131 


0.00 


J 


0.12 


-1.3 


2.2 


-1.6 


132 


1.2 


J 


0.049 


-2.3 


2.2 


-1.1 


133 


3.1 


J 


0.23 


-2.3 


2.2 


-1.0 
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Table 13— Continued 



Index 


Ay 


Anorm 




log(idisc/i'star) 


•^turnoff 


CKexcess 






[/im] 






[^m] 




137 


5.9 




0.72 


-1.1 


3.6 


-0.80 


138 


3.5 


J 


130. 


1.0 


4.5 


0.78 
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Table 14. SED modelling parameters for Class III YSOs 



Index Ay Anorm istar log(L;jisc/Lgtar ) Atm-noff ^excess 

[/im] [Lq] lum] 



21 


1.3 


J 





58 


2.9 


J 





62 


4.5 


J 


150 


95 


1.4 


J 





107 


1.8 


J 





108 


4.5 


J 





112 


0.00 


J 





116 


0.00 


J 





117 


3.3 


J 


25 


119 


0.39 


J 





129 


0.00 


J 


2 


136 


4.8 


J 




139 


1.4 


J 






46 


-0.57 


24. 


35 


-0.99 


7.9 




-2.1 


25. 


16 


-1.1 


7.9 


36 


-0.41 


7.9 


19 


-1.0 


5.7 


35 


-0.75 


2.2 


15 


-0.62 


7.9 




-2.2 


2.2 


82 


-3.2 


2.2 


7 


-1.5 


24. 



1.7 

.37 -2.1 7.9 
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Table 15. Statistics of star forming dark clouds in Cepheus 



Dark Cloud 




Spitzer identified YSOs 


a 


non-S. 


YSOs'' 


Nyso 


Av 


Area'^ 




SFE 




I 


F 


II 


III 


Tot. 


I/Tot. 


Inc. 


Tot. 


Ay >5 


Max 


Mean 




















[%] 








[mag] 


[mag] 


[pc-^] 


[Me] 


[%] 


Prestellar Cores 


L1155 







































L1241 











1 


1 














11.2 


6.2 


0.058 


150 





L1247 





























8.9 


6.2 


0.015 


38 





Sub Total 











1 


1 

















0.073 


188 




Protostellar Cores 


(peak f 


. < 


disruption 


density) 




















L1148+L1147 


1 


1 


2 





4 


25 








2 


8.8 


6.0 


0.051 


150 


1 


L1152 


1 





2 





3 


33 





















L1157 




















1 


1 















L1172 


1 


1 


3 


1 


6 


17 








4 


13.1 


7.0 


0.040 


107 


2 


L1221 


3 











3 


100 





















L1251W 


2 


1 


1 





4 


50 








4 


11.7 


6.7 


0.040 


110 


2 


Sub Total 


8 


3 


8 


1 


20 




1 


1 


10 






0.131 


367 


1 


Cores with YSO Groups/Clusters (peak p* > disruption 


density) 
















NGC 7023 Core 


6 


3 


26 


2 


37 


16 


4 


6 


34 


14.9 


7.0 


0.040 


107 


14 


. . . Loose 


5 


3 


23 


1 


32 


16 


4 


6 






5.7 


0.078 


169 


9 


. . . Tight A 


1 


2 


7 





10 


10 


3 


3 




8.5 


6.0 


0.004 


8 


38 


. . . Tight B 


2 


1 


2 





5 


40 










10.7 


6.3 


0.001 


2 


57 


L1228N Core 


2 


1 


7 





10 


20 








6 


12.9 


7.7 


0.059 


84 


3 


. . . Loose 


2 


1 


2 





5 


40 












7.0 


0.021 


27 


9 


. . . Tight 


2 


2 


2 





6 


33 












9.1 


0.001 


2 


65 


L1228S Core 





2 


9 





11 





2 


3 


8 


10.0 


6.5 


0.028 


34 


10 


. . . Loose 


4 


2 


7 





13 


31 


2 


3 






4.7 


0.016 


14 


31 


. . . Tight 


4 





2 





6 


67 





1 






6.2 


0.001 


1 


69 


L1251A Core 


1 


4 


11 





16 


6 


2 


2 


17 


14.6 


7.3 


0.026 


79 


10 


. . . Loose 


1 


3 


8 





12 


8 


1 


1 






6.4 


0.016 


42 


12 


. . . Tight 





1 


2 





3 





1 


1 






6.9 


0.001 


3 


32 


L1251B Core 


4 





12 


3 


19 


21 





1 


17 


22.1 


8.3 


0.042 


144 


6 


. . . Loose 


4 





9 


3 


16 


25 












7.6 


0.022 


68 


11 


. . . Tight 


3 





2 


2 


7 


43 












14.3 


0.002 


13 


22 


Sub Total Core 


13 


10 


65 


5 


93 


14 


8 


12 


82 






0.195 


447 


8 


. . . Loose 


16 


9 


49 


4 


78 


21 


7 


10 








0.153 


320 


11 


. . . Tight 


12 


6 


17 


2 


37 


32 


4 


5 








0.010 


29 


39 


Off Cloud 






























Sub Total 





1 


14 


4 


19 





1 


2 














Grand Total 


21 


14 


87 


11 


133 


16 


10 


15 


92 






0.400 


1003 


4 



''The number of Spitzer identified YSO candidates by Class (as determined by the aijs spectral index) associated with the named 
core (within the projected Ay = 1 contour) or cluster (within the relevant tight or loose contour). 

''The number of non-Spitzer identified YSOs as listed in Table [T6l The column labeled 'Inc.' lists the number of non-Spitzer 
identified YSOs that were coincident with entries in the Cepheus catalog. The column labeled 'Tot.' lists the total number of 
YSOs missed by the Spitzer YSO schemes in that region. 



'^The area of the region used to calculate the mean Ay, the cloud mass Mdoud, and the star formation efficiency (SFE). For a 



- 107- 



core, this is the area that is enclosed by the Ay = 5 contour. For a loose and tight cluster, this is the area that is enclosed by the 
p» = 1.0 and 25 Mq pc~^ contours respectively. 
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Table 16. Non-identified SIMBAD YSOs 



Source Name'' 


Right Ascension 


Declination 


Origin*^ 


Region 


Spitzer 


Spitzer 




(2000) 


(2000) 






Catalog. 


YSO Id^ 


L1157 


20:39:06.2 


-1-68:02:15 


S 


L1157 


Y 


134 


PV Cep 


20:45:53.9 


-1-67:57:39 


S 


Off-cloud 


Y 


135 


IRAS F20597+6800 


21:00:18.4 


-f68:12:02 


s 


NGC 7023 


N 




L1228 Sl(b) 


21:01:25.1 


-1-77:06:54 


K 


L1228 


Y 


136 


NGC 7023 RS 2 


21:01:27.1 


-f68:10:38 


K 


NGC 7023 


Y 


137 


IRAS F21022+7651 


21:01:29.4 


-1-77:03:45 


S 


L1228 


N 




HD 200775 


21:01:36.9 


-1-68:09:48 


s 


NGC 7023 


Y 


138 


NGC 7023 RS 3 


21:01:43.6 


-1-68:09:36 


K 


NGC 7023 


Y 


139 


IRAS F21025+7645 


21:01:49.2 


-1-76:57:51 


S 


L1228 


Y 


140 


NGC 7023 S T 


21:02:07.9 


-1-68:09:16 


s 


NGC 7023 


N 




IRAS F21025+6801 


21:03:12.0 


-f68:13:02 


s 


NGC 7023 


Y 


141 


CSC 04600-02450 


22:34:11.9 


-f 75:10:10 


S,K 


L1251A 


Y 


142 


TW96 Nla 


22:35:23.3 


-1-75:17:09 


S 


L1251A 


Y 


143 


XMMU J223727.7+751525 


22:37:27.5 


4-75:15:25 


S,K 


L1251B 


N 




KP93 3-19 


22:40:58.5 


-1-75:03:56 


S 


L1251 (Off-cloud) 


N 





''Name of the source in either the SIMBAD database or lKun et"al1 (|2008l '). 

''The o rigin of the YSO candidate. S - a YSO or pre-main sequence listing in the SIMBAD database, K - a T Tauri 
star from lKun et al.l l l2008h 

'^Does the source have an entry in the Cepheus Spitzer catalog? In this instance, the only condition for inclusion is 
to have a line entry in the Spitzer catalog irrespective of whether that source was formally detected in any particular 
band. 

'^The YSO identification number from Table|6]for sources that are found in the Spitzer Cepheus catalog. 
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Table 17. Star formation statistics for the Cepheus Flare 



Boundary 


Nyso 


CI"- 






Nyso/A 


SFR'= 


SFR/A 


Mcloud 


Mcloud M 






[deg2] 




[deg-2] 


[pc-^] 


[Mq Myr-l] 


[M0 Myr-l pc-^] 


[Me] 


[M0 pc-2] 


Av < 1 


19 


6.2 


170 


3.1 


0.11 










Ay > 1 


114 


2.1 


59 


53 


1.9 


29 


0.49 






Av > 5 


92 


0.40 


11 


230 


8.4 


58 


2.1 


1003 


91 



''The solid angle enclosed by the contour 

cross sectional area assuming an average distance of 300 pc (1 sqr deg = 27.4 pc ) 

"^The star formation rate (SFR) calculated assuming a star formation timescale of 2 Myr and an average stellar mass of 
0.5 M0 
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Table 18. Known outflows in Cepheus 



Region 


Lit. Name 


Assoc. YSO 


HH No. 


Example References 


OfT-core 


PV Cep 


135 


215, 315, 415 


^oodman & Arce 2004 


L1152 


GM 3-12 


3 


376 


Rciourth ct al. 1997; Movsessian et al. 2004 


L1152 


L1152 SMM 


1 




Bontemps et al. 1996 


L1157 


L1157 IRS 


134 


375 


Umemoto et al. 1992: Loonev et al. 2007 


L1172D 


L1172 SMM 1 


49 




Mvers et al. 1988 


L1174 


HD 200775 


136 




Puente et al. 1998 


L1221 


L1221 IRS 1 


105 




Umemoto et al. 1991: Lee & Ho 2005 


L1228N 


HH 200 IRS 


7 


200 


Ballv et al. 1995 


L1228N 


HH 199 1 


9 


199 


Ballv et al. 1995 


L1251B 




Undetermined 


189 


Eiroa et al. 1994 


L1251A 




143 


149 


Sato & Pukui 1989: Balazs et al. 1992 



